Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,761 T→C 5.4% intergenic (+93/‑67) rpsA → / → ihfB 30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,7610TC 5.4% 365.7 / 5.5 110intergenic (+93/‑67)rpsA/ihfB30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (45/59);  new base C (2/4);  total (47/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 110 total aligned reads displayed.

AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAA  >  NC_000913/963615‑963908
                                                                                                                                                  |                                                                                                                                                   
aaCTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTAc                                                                                                                                                   >  4:888929/1‑149 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgc                                                                                                                                             >  1:317336/1‑148 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                           >  4:121833/1‑149 (MQ=255)
         aacaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCt                                                                                                                                          <  4:777919/149‑1 (MQ=255)
          acaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCtt                                                                                                                                         <  3:857932/149‑1 (MQ=255)
           caacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                         <  2:634453/148‑1 (MQ=255)
           caacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                         <  1:302061/148‑1 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGACGGCTACGGCCGCCCTTaa                                                                                                                                       >  4:651403/1‑149 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACCTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTaa                                                                                                                                       <  3:516501/149‑1 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                   <  1:23990/148‑1 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAAt                                                                                                                                  <  3:152406/149‑1 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAAt                                                                                                                                  <  4:744060/149‑1 (MQ=255)
                       tGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcagca                                                                                                                            >  3:126573/1‑149 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                      >  4:695923/1‑149 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                      >  3:650325/1‑149 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                      >  4:634975/1‑149 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagc                                                                                                                       >  2:746520/1‑148 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCATCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAACCACTAAGGGCGGCTACGGGCGCCCTTAATCAATGCAGCAAcagc                                                                                                                       <  2:383373/148‑1 (MQ=255)
                                 aaaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGcc                                                                                                                   <  1:138655/148‑1 (MQ=255)
                                 aaaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCg                                                                                                                  <  3:163627/149‑1 (MQ=255)
                                  aaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCg                                                                                                                  >  2:870430/1‑148 (MQ=255)
                                   agcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGc                                                                                                                 >  2:808863/1‑148 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCtt                                                                                                               >  3:893082/1‑149 (MQ=255)
                                       gcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTaa                                                                                                             <  2:196287/148‑1 (MQ=255)
                                         tAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAtt                                                                                                           <  1:241035/148‑1 (MQ=255)
                                          aaaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAttt                                                                                                          <  1:938753/148‑1 (MQ=255)
                                          aaaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAtt                                                                                                           <  2:796874/147‑1 (MQ=255)
                                          aaaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAtt                                                                                                           >  1:796874/1‑147 (MQ=255)
                                           aaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTg                                                                                                         <  2:522734/148‑1 (MQ=255)
                                             ggCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCt                                                                                                      <  4:838770/149‑1 (MQ=255)
                                               cGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCAATAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCttt                                                                                                    >  3:661597/1‑149 (MQ=255)
                                                 aGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTaa                                                                                                  <  3:506909/149‑1 (MQ=255)
                                                 aGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTa                                                                                                   >  1:876574/1‑148 (MQ=255)
                                                 aGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTa                                                                                                   <  1:133925/148‑1 (MQ=255)
                                                    aaTTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTCGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAgg                                                                                                >  2:631652/1‑148 (MQ=255)
                                                         tctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                          >  4:830312/1‑149 (MQ=255)
                                                          ctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                         <  4:709660/149‑1 (MQ=255)
                                                          ctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                         <  4:116670/149‑1 (MQ=255)
                                                          ctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                         <  4:358290/149‑1 (MQ=255)
                                                          ctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                          <  2:441234/148‑1 (MQ=255)
                                                          ctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                          >  2:847213/1‑148 (MQ=255)
                                                            gACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgcc                                                                                                                                               <  3:936321/93‑1 (MQ=255)
                                                                ctTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAAt                                                                                    >  2:28207/1‑148 (MQ=255)
                                                                ctTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAAt                                                                                    <  2:128588/148‑1 (MQ=255)
                                                                ctTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                         <  2:659918/95‑1 (MQ=255)
                                                                ctTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                         >  1:659918/1‑95 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                   <  1:89588/100‑1 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                   >  2:89588/1‑100 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATc                                                                                   <  2:364533/148‑1 (MQ=255)
                                                                  tCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                 <  4:893082/149‑1 (MQ=255)
                                                                    ggtaTTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGa                                                                               <  3:634975/146‑1 (MQ=255)
                                                                    gggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATg                                                                                <  1:870430/148‑1 (MQ=255)
                                                                    gggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATg                                                                                <  2:65675/148‑1 (MQ=255)
                                                                    gggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGa                                                                               <  3:455979/149‑1 (MQ=255)
                                                                      gATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGAcc                                                                             >  4:12780/1‑149 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTc                                                                        <  4:53720/149‑1 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGaa                                                                    <  4:650325/149‑1 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGa                                                                     <  1:132033/148‑1 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGa                                                                     <  2:515278/148‑1 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGa                                                                     <  1:363483/148‑1 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGaa                                                                    >  1:609858/1‑148 (MQ=255)
                                                                                   aaGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTg                                                                 <  1:28207/148‑1 (MQ=255)
                                                                                   aaGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGa                                                                <  4:725895/149‑1 (MQ=255)
                                                                                           tgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGa                                                        <  3:90887/149‑1 (MQ=255)
                                                                                               tttACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTg                                                    >  3:48903/1‑149 (MQ=255)
                                                                                               tttACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTg                                                    >  3:290183/1‑149 (MQ=255)
                                                                                                   cTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                         >  4:291721/1‑108 (MQ=255)
                                                                                                   cTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                         <  3:291721/108‑1 (MQ=255)
                                                                                                     tGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGAGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGGCCAGGTTAAAATTGATAGAAACACTTGCCAca                                               >  2:160238/1‑147 (MQ=255)
                                                                                                      gACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCa                                             >  3:296837/1‑149 (MQ=255)
                                                                                                          gATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCa                                          <  2:809700/148‑1 (MQ=255)
                                                                                                               cAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGca                                    <  4:290183/149‑1 (MQ=255)
                                                                                                                    ttCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGcc                                                  <  4:414917/130‑1 (MQ=255)
                                                                                                                    ttCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGcc                                                  >  3:414917/1‑130 (MQ=255)
                                                                                                                      cGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTcc                              <  1:420090/148‑1 (MQ=255)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCg                            >  4:238092/1‑149 (MQ=255)
                                                                                                                        tCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGc                           >  3:663866/1‑149 (MQ=255)
                                                                                                                        tCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGc                           >  4:292292/1‑149 (MQ=255)
                                                                                                                          cTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGc                                     <  1:891971/137‑1 (MQ=255)
                                                                                                                          cTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGc                                     >  2:891971/1‑137 (MQ=255)
                                                                                                                           tGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCa                         <  2:317336/148‑1 (MQ=255)
                                                                                                                           tGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATCGAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCTTGTCCAAGTCAAAAGTGATAGTAAGACCTGCCACCCCGCAATCGCACATTCCCTCCaa                        <  4:143442/149‑1 (MQ=255)
                                                                                                                             tAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAg                       >  1:12916/1‑148 (MQ=255)
                                                                                                                              aaTCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGa                      >  2:920292/1‑148 (MQ=255)
                                                                                                                                 cAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACgg                   >  1:731582/1‑148 (MQ=255)
                                                                                                                                 cAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGt                  <  3:238092/149‑1 (MQ=255)
                                                                                                                                  aaGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCATCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGt                  <  1:808863/148‑1 (MQ=255)
                                                                                                                                  aaGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGtt                 <  3:695923/149‑1 (MQ=255)
                                                                                                                                  aaGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGt                  <  1:688058/148‑1 (MQ=255)
                                                                                                                                    gCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTg                >  1:409810/1‑148 (MQ=255)
                                                                                                                                       ttAAGGGCGGCCGTAGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAAcc                                                                                            >  4:874181/2‑69 (MQ=37)
                                                                                                                                       cTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAAc                                                                                             <  3:279412/68‑1 (MQ=255)
                                                                                                                                        tAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCGCATTCCCGCCAAGACGGTTGAAGa            >  1:264602/1‑148 (MQ=255)
                                                                                                                                        tAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGa            >  1:499978/1‑148 (MQ=255)
                                                                                                                                        tAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGAt           <  3:12780/149‑1 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGc         >  4:576119/1‑149 (MQ=255)
                                                                                                                                           gggCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa        >  4:161231/1‑149 (MQ=255)
                                                                                                                                            ggcggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa        <  2:891213/148‑1 (MQ=255)
                                                                                                                                                 cTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaaaa  <  4:16196/149‑1 (MQ=255)
                                                                                                                                                 cTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaaa   <  2:609858/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAA  >  NC_000913/963615‑963908

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: