Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,473 G→T 5.4% intergenic (+55/‑33) yaiT → / → yaiV pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,4730GT 5.4% 752.1 / 18.9 209intergenic (+55/‑33)yaiT/yaiVpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (78/118);  new base T (3/10);  total (81/128)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

AGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394324‑394622
                                                                                                                                                      |                                                                                                                                                     
aGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTTGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAgg                                                                                                                                                       <  2:593648/150‑1 (MQ=255)
   aCTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGAtt                                                                                                                                                    <  2:1898942/150‑1 (MQ=255)
    cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGAttt                                                                                                                                                   >  2:1334021/1‑150 (MQ=255)
    cTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGAttt                                                                                                                                                   >  2:1695488/1‑150 (MQ=255)
     ttCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTg                                                                                                                                                  >  2:3125215/1‑150 (MQ=255)
      tCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGa                                                                                                                                                 >  1:680865/1‑150 (MQ=255)
      tCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGa                                                                                                                                                 <  1:1385607/150‑1 (MQ=255)
       cAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGAc                                                                                                                                                <  2:1919698/150‑1 (MQ=255)
         gcgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACgg                                                                                                                                              >  2:651825/1‑150 (MQ=255)
         gcgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACgg                                                                                                                                              <  2:1545983/150‑1 (MQ=255)
         gcgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACgg                                                                                                                                              <  2:2891531/150‑1 (MQ=255)
         gcgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACgg                                                                                                                                              <  2:3812610/150‑1 (MQ=255)
          cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACggg                                                                                                                                             <  2:3676819/150‑1 (MQ=255)
          cgcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAATCCTTGACggg                                                                                                                                             <  2:2324062/150‑1 (MQ=255)
           gcCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTATTCACTCCGAAGAAATACTGTTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGc                                                                                                                                            <  1:2546861/150‑1 (MQ=255)
            ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                           >  2:3769257/1‑150 (MQ=255)
            ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                           <  1:1264290/150‑1 (MQ=255)
            ccTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCa                                                                                                                                           <  1:3272898/150‑1 (MQ=255)
                 aCCGATGCTAACTACCTCGGTGGTGGTGACGGAGATCAAGACGGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCAGTACt                                                                                                                                      <  2:3885740/150‑1 (MQ=255)
                  ccGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAATCCTTGACGGGCATTattt                                                                                                                                     <  1:2438178/150‑4 (MQ=255)
                   cGATGCTACCTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAATCCTTGACGGGCATTattta                                                                                                                                    <  2:1930334/150‑5 (MQ=255)
                      tGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAATTATACCTGGTAATATTCTTCCCTCCGAAGAAATACTGCTTATTTAATCT‑AAATACTGCCCGTCAAGGATTTGACGGGCATTACTGCaaa                                                                                                                                 <  2:2333168/150‑1 (MQ=255)
                      tGCTAACCACCGCGGGGGTGGTGACGTAGAGCAAGACTGGTCCGCGAACGTGGGTGTTCAATCTACCTGGTAATCTTCTTCCCTCCGAAGAAATACTGGTAATTTAATCT‑AACTACTGCCCGTCAAGGATTTGACGGGCATTACTGCaaa                                                                                                                                 <  2:2558364/150‑1 (MQ=255)
                            cTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCATGGATTTGACGGGCATTACTGCAAAGGAcgc                                                                                                                           <  2:3649984/150‑1 (MQ=255)
                             tACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGAcgcg                                                                                                                          >  1:1151888/1‑150 (MQ=255)
                               ccTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTATTATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATACTGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCa                                                                                                                        <  2:2849188/150‑1 (MQ=255)
                                 tCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTTGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCaaa                                                                                                                      <  2:1508756/150‑1 (MQ=255)
                                 tCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCaaa                                                                                                                      >  1:1425693/1‑150 (MQ=255)
                                 tCGGTGGGGGTGACGTAGAACAAGACTGGTCCGCGAACGTGGGTTTTAAATATACGGGTTAATATTCTTCACTCCTAAGAAATACTTGTAATTTAATAT‑AAAAAATGCCCGTCATGGATTTGACGGGCATTACTGCAAAGGACGCGCaaa                                                                                                                      <  1:3350630/150‑1 (MQ=255)
                                  cGGTGGTGGTTACTTAGATCAAGACTGGTCCGCTAACTTGGGTGTTAAATATACCTGTTAATATTCTTCAATCCGAATAAATACTGTTAATTTAATAT‑AAATAATGCCCGTAAAGTTTTTGACGGGCATTACTGCAAAGGACGCGCAAAt                                                                                                                     <  1:937603/150‑1 (MQ=255)
                                        tggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATc                                                                                                               >  1:834014/1‑150 (MQ=255)
                                         ggtgACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCt                                                                                                              <  2:335094/150‑1 (MQ=255)
                                            gACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTa                                                                                                           >  2:3702887/1‑150 (MQ=255)
                                                  gATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTaaa                                                                                                     >  1:3002237/1‑150 (MQ=255)
                                                    tCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAcc                                                                                                   >  2:1983784/1‑150 (MQ=255)
                                                    tCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAcc                                                                                                   >  1:3414319/1‑150 (MQ=255)
                                                          cTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCa                                                                                             <  1:866447/150‑1 (MQ=255)
                                                          cTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCa                                                                                             <  1:2573105/150‑1 (MQ=255)
                                                             gTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGa                                                                                          >  2:1078473/1‑150 (MQ=255)
                                                              tCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGTTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGaa                                                                                         <  1:1213535/150‑1 (MQ=255)
                                                              tCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGaa                                                                                         <  1:1234884/150‑1 (MQ=255)
                                                                 gcgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttt                                                                                      <  2:210941/150‑1 (MQ=255)
                                                                 gcgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAAttt                                                                                      <  2:3452753/150‑1 (MQ=255)
                                                                  cgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTg                                                                                     <  2:2712609/150‑1 (MQ=255)
                                                                  cgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTg                                                                                     <  2:3949014/150‑1 (MQ=255)
                                                                  cgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTg                                                                                     <  1:331336/150‑1 (MQ=255)
                                                                  cgAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTg                                                                                     <  2:2269360/150‑1 (MQ=255)
                                                                                aaaTATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCGAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATaaa                                                                       >  1:2140217/1‑150 (MQ=255)
                                                                                aaaTATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATaaa                                                                       <  1:2518428/150‑1 (MQ=255)
                                                                                 aaTATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AATTAATGCCCGTCAAGGATTTGACGGGCATTTCTGCAATGGACGCGCAAATGTTATCTGTTGTTAAACCTCTTCTGGAATTTGGTAAGCTCGATaatt                                                                      >  1:1970296/1‑148 (MQ=255)
                                                                                     taCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGttt                                                                  <  2:793143/150‑1 (MQ=255)
                                                                                      aCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTg                                                                 >  1:3051180/1‑150 (MQ=255)
                                                                                      aCCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAAGGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTg                                                                 >  1:1996802/1‑150 (MQ=255)
                                                                                          ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             >  2:3206458/1‑150 (MQ=255)
                                                                                          ggTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCa                                                             >  1:1548820/1‑150 (MQ=255)
                                                                                           gTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAg                                                            <  1:3701738/150‑1 (MQ=255)
                                                                                            tAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGa                                                           >  1:3234033/1‑150 (MQ=255)
                                                                                                 ttcttcACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACgg                                                      <  1:3525346/150‑1 (MQ=255)
                                                                                                   cttcACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTa                                                    >  1:110337/1‑150 (MQ=255)
                                                                                                    ttcACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAc                                                   >  1:2962397/1‑150 (MQ=255)
                                                                                                    ttcACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAc                                                   >  1:628554/1‑150 (MQ=255)
                                                                                                    ttcACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAc                                                   >  2:1697965/1‑150 (MQ=255)
                                                                                                       aCTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTAcgcg                                                >  1:69190/1‑150 (MQ=255)
                                                                                                             aagaaATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGa                                          <  1:3166936/150‑1 (MQ=255)
                                                                                                             aagaaATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGa                                          <  1:3166915/150‑1 (MQ=255)
                                                                                                              agaaATACTGGTAATTTAATCTAAAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAAGCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGa                                          <  2:695815/150‑1 (MQ=255)
                                                                                                                aaaTACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGtt                                       <  1:3769257/150‑1 (MQ=255)
                                                                                                                aaaTACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGtt                                       <  2:3051180/150‑1 (MQ=255)
                                                                                                                  aTACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTta                                     <  1:2645185/150‑1 (MQ=255)
                                                                                                                    aCTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaat                                   <  1:1294679/150‑1 (MQ=255)
                                                                                                                     cTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaata                                  <  2:1765402/150‑1 (MQ=255)
                                                                                                                      tGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaataa                                 >  2:3499896/1‑150 (MQ=255)
                                                                                                                      tGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTtaataa                                 >  2:2630033/1‑150 (MQ=255)
                                                                                                                        gTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATg                               <  1:1601317/150‑1 (MQ=255)
                                                                                                                           aTTTAATCT‑AAATAATGCCCTTCAAGTATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGaaa                            <  2:747133/150‑1 (MQ=255)
                                                                                                                            tttAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGaaat                           >  1:1064596/1‑149 (MQ=255)
                                                                                                                                aTCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCaa                       <  2:3624903/150‑1 (MQ=255)
                                                                                                                                 tCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAg                      >  2:1928040/1‑150 (MQ=255)
                                                                                                                                  cT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGt                     <  2:3603794/150‑1 (MQ=255)
                                                                                                                                   t‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    >  1:703015/1‑150 (MQ=255)
                                                                                                                                   t‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    >  2:3119660/1‑150 (MQ=255)
                                                                                                                                   t‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGtt                    >  1:2858247/1‑150 (MQ=255)
                                                                                                                                          atGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAttt              <  2:3865538/150‑1 (MQ=255)
                                                                                                                                           tGCCCGTCAAATCCTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATAtttt             >  1:3921800/1‑150 (MQ=255)
                                                                                                                                            gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc            <  1:463029/150‑1 (MQ=255)
                                                                                                                                            gCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTc            <  1:2615305/150‑1 (MQ=255)
                                                                                                                                             cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  1:548492/1‑150 (MQ=255)
                                                                                                                                             cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  1:1681205/1‑150 (MQ=255)
                                                                                                                                             cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  1:1823751/1‑150 (MQ=255)
                                                                                                                                             cccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTcc           >  1:1823727/1‑150 (MQ=255)
                                                                                                                                              ccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCa          >  1:1929008/1‑150 (MQ=255)
                                                                                                                                              ccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCa          <  2:245519/150‑1 (MQ=255)
                                                                                                                                              ccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCa          >  1:2788625/1‑150 (MQ=255)
                                                                                                                                              ccGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCa          <  1:2934817/150‑1 (MQ=255)
                                                                                                                                               cGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAg         <  2:2353728/150‑1 (MQ=255)
                                                                                                                                                 tCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTg       <  1:3077020/150‑1 (MQ=255)
                                                                                                                                                   aatgATTTGTCGGGCATTACTTCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGAt     <  2:2495542/147‑1 (MQ=255)
                                                                                                                                                   aaGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGAt     <  1:512631/150‑1 (MQ=255)
                                                                                                                                                    aGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATg    >  2:3264915/1‑150 (MQ=255)
                                                                                                                                                      gATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTc  >  1:328196/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGTAATATTCTTCACTCCGAAGAAATACTGGTAATTTAATCT‑AAATAATGCCCGTCAAGGATTTGACGGGCATTACTGCAAAGGACGCGCAAATGTTATCTGTAGTTAAACCTCTTCAGGAATTTGGTAAGCTCGATAAATGTTTGTCCAGATACGGTACGCGCTTCGAGTTTAATAATGAAAAGCAAGTTATATTTTCCAGTGATGTC  >  NC_000913/394324‑394622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑