Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,421 A→T 6.7% intergenic (‑439/+164) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4210AT 6.7% 348.2 / 11.7 120intergenic (‑439/+164)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (56/56);  new base T (2/6);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

Only 100 of 125 total aligned reads displayed.

TTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTC  >  NC_000913/987273‑987570
                                                                                                                                                    |                                                                                                                                                     
    cTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATg                                                                                                                                                  <  2:638373/150‑1 (MQ=255)
      tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                                <  1:415507/150‑1 (MQ=255)
       atatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGATAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCgg                                                                                                                                               <  2:1339756/150‑1 (MQ=255)
       atatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCgg                                                                                                                                               >  1:1684562/1‑150 (MQ=255)
        tatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGa                                                                                                                                              <  2:607/150‑1 (MQ=255)
         atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTAATCGACCTGGCACAAAATCTTCTTTCAGAAATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCAGCAATAAGAATGCGGAt                                                                                                                                             >  1:746213/1‑150 (MQ=255)
          tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                            >  1:1855195/1‑150 (MQ=255)
          tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTATCAGACATCCAGAGTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                            <  1:1218817/150‑1 (MQ=255)
             aaaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGt                                                                                                                                         >  1:1733719/1‑150 (MQ=255)
             aaaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGt                                                                                                                                         <  1:1350218/150‑1 (MQ=255)
              aaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTTGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGtt                                                                                                                                        <  2:1406330/150‑1 (MQ=255)
                ttAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGtttt                                                                                                                                      <  2:1125823/150‑1 (MQ=39)
                 tAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttt                                                                                                                                     <  1:1823133/150‑1 (MQ=39)
                 tAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttt                                                                                                                                     <  2:1529959/150‑1 (MQ=39)
                  aaCTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGtttttt                                                                                                                                    <  2:1972695/150‑2 (MQ=39)
                    cTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttttct                                                                                                                                  <  2:893069/150‑4 (MQ=39)
                      aaaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGTATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGcgcg                                                                                                                                >  2:755653/1‑150 (MQ=255)
                           gAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAAt                                                                                                                           <  2:1536183/150‑1 (MQ=255)
                               ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                       <  2:1325909/150‑1 (MQ=255)
                                 ttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCa                                                                                                                     <  2:1752403/150‑1 (MQ=255)
                                  tGTTTCTGATGGAAACTTTATCGACCTGGCACATAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATCAGAATGCGGATTAGTTTTTGCGCGTAAATTTTGCaa                                                                                                                    <  1:1221503/150‑1 (MQ=255)
                                   gTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAAt                                                                                                                   >  1:575742/1‑150 (MQ=255)
                                     ttCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCAGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGt                                                                                                                 >  1:97923/1‑150 (MQ=255)
                                     ttCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATACAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGt                                                                                                                 >  1:1037873/1‑150 (MQ=255)
                                      tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGAGATCCAGAATGCCAAATAATAGTTTGAAATTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                                <  2:1290950/150‑1 (MQ=255)
                                      tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGAATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                                >  2:761556/1‑150 (MQ=255)
                                       cTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCATAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCTATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTaa                                                                                                               <  2:1731200/150‑1 (MQ=255)
                                        tGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                              <  1:1404796/150‑1 (MQ=255)
                                        tGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                              >  2:1840123/1‑150 (MQ=255)
                                        tGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                              >  2:1230675/1‑150 (MQ=255)
                                          aTGGAAACTTTATCTACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAatat                                                                                                            >  2:516990/1‑150 (MQ=255)
                                             gAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTcc                                                                                                         >  2:1839187/1‑150 (MQ=255)
                                               aaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGt                                                                                                       <  2:38102/150‑1 (MQ=255)
                                               aaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGt                                                                                                       <  2:31550/150‑1 (MQ=255)
                                               aaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGt                                                                                                       <  1:440548/150‑1 (MQ=255)
                                                  tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAc                                                                                                    <  2:1084983/150‑1 (MQ=255)
                                                    tATCGCCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGAAATATTCCGTAACAg                                                                                                  >  2:1345039/1‑150 (MQ=255)
                                                    tATCGACATGGCACAAAATCTTCTTTCAGACATACAGAATGCCAACAAATAGTATGAATATTGGGCTTTTAGCTAAGACAAAACTAATCCGAAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAg                                                                                                  <  1:304358/150‑1 (MQ=255)
                                                            tGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATATTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCg                                                                                          <  2:1752157/150‑1 (MQ=255)
                                                            tGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCg                                                                                          <  2:1754737/150‑1 (MQ=255)
                                                            tGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCg                                                                                          >  1:1267329/1‑150 (MQ=255)
                                                             ggCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGAAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATAGt                                                                                         >  1:735068/1‑150 (MQ=255)
                                                              gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGtt                                                                                        >  2:1981467/1‑150 (MQ=255)
                                                              gCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGtt                                                                                        >  2:883140/1‑150 (MQ=255)
                                                               cacaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTa                                                                                       >  2:1039854/1‑150 (MQ=255)
                                                                      tcttctTTCAGACATCCATAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCggg                                                                                <  1:210564/150‑1 (MQ=255)
                                                                      tcttctTTCAGACATCCAGAACGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCggg                                                                                <  1:340990/150‑1 (MQ=255)
                                                                        ttcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGGTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGAGTAAATGTTGCAATGTAATATTCCGTAACAGGATGAGCGTTATGGCGGGGt                                                                              <  2:1219139/150‑1 (MQ=255)
                                                                         tctaTCAGACATCCAGAATGACAAAAAATAGTATGAAAATTATGCTATTAGCTAAGAAAAAACTAATCCTCAATTAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAAAGTAATATTCCGTAACAGGATGATAGTTATGGCTGGGtt                                                                             >  2:1040088/1‑150 (MQ=255)
                                                                          cttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGAGCGTAAATGTTGCAATGTAATATTCCGTAACAGTATGATCGTTATTGCTGGGTTg                                                                            >  2:529507/1‑150 (MQ=255)
                                                                            ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGTGGTTGCa                                                                          >  1:1429139/1‑150 (MQ=255)
                                                                            ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                          >  2:1534007/1‑150 (MQ=255)
                                                                            ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                          <  1:1477909/150‑1 (MQ=255)
                                                                             tCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGTATTAGTTTTTGCGCGTAAAAGTTTCAATGTAATATTCCGTAACAGGATGATCGTTATGGCCGGGTTGCaa                                                                         >  1:1518811/1‑150 (MQ=255)
                                                                                gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                      <  1:734111/150‑1 (MQ=255)
                                                                                gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGATTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATa                                                                      <  1:1030002/150‑1 (MQ=255)
                                                                                 aCATCCAGAATGCCAAAAAATAGTATTAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAg                                                                     <  2:1800424/150‑1 (MQ=255)
                                                                                 aCATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAg                                                                     <  1:1839187/150‑1 (MQ=255)
                                                                                 aCATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAg                                                                     <  2:28295/150‑1 (MQ=255)
                                                                                  cATCCAGTATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGa                                                                    <  1:755653/150‑1 (MQ=255)
                                                                                  cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATTTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGa                                                                    <  2:1253569/150‑1 (MQ=255)
                                                                                  cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGa                                                                    >  2:670668/1‑150 (MQ=255)
                                                                                    tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAtt                                                                  >  2:69907/1‑150 (MQ=255)
                                                                                    tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAtt                                                                  >  1:1190144/1‑150 (MQ=255)
                                                                                      cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGc                                                                >  1:1589254/1‑150 (MQ=255)
                                                                                        gAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCtt                                                              <  1:441576/150‑1 (MQ=255)
                                                                                         aaTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                             <  1:911709/150‑1 (MQ=255)
                                                                                            gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                          >  1:1601040/1‑150 (MQ=255)
                                                                                            gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                          >  1:1856659/1‑150 (MQ=255)
                                                                                             ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                         >  1:350767/1‑150 (MQ=255)
                                                                                             ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                         >  1:1629919/1‑150 (MQ=255)
                                                                                             ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTAGACtt                                                         >  2:1288003/1‑150 (MQ=255)
                                                                                                aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaac                                                      <  1:1847319/150‑1 (MQ=255)
                                                                                                 aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaaca                                                     <  1:1072240/150‑1 (MQ=255)
                                                                                                     tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaaca                                                 <  1:1230675/150‑1 (MQ=255)
                                                                                                      aGTATGAAAATTTTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATCCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAg                                                <  2:288889/150‑1 (MQ=255)
                                                                                                         aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGAGTAAATGTTGCAATGTAATATTCCTTAACATGATGATCGTTATGGCGGGGTTGCAAATATATTGCTTGACTTAACAAACAGATa                                             >  1:1299007/1‑150 (MQ=255)
                                                                                                              aaTTGTGATATTATCTAAGAAAAAACTAATCAGCAATAATAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATAGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                        <  1:230842/150‑1 (MQ=255)
                                                                                                                ttGTGCTATTAGCTAAGAAAAAACTAATCCGAAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGAc                                      <  1:953720/150‑1 (MQ=255)
                                                                                                                  gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAg                                    <  1:282268/150‑1 (MQ=255)
                                                                                                                       aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAg                               <  1:1522478/150‑1 (MQ=255)
                                                                                                                        ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGa                              >  1:1709156/1‑150 (MQ=255)
                                                                                                                           gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                           >  1:596998/1‑150 (MQ=255)
                                                                                                                            cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                          >  2:1899717/1‑150 (MQ=255)
                                                                                                                            cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                          <  1:1756777/150‑1 (MQ=255)
                                                                                                                              aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                        <  2:1267329/150‑1 (MQ=255)
                                                                                                                                gcaaaaaCTAATCCGCATTCTTAGTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATTGATTGCTTGACTTAACAAACAGATAACTTGACAGAATAGATaaaaaaa                      >  1:903482/3‑150 (MQ=38)
                                                                                                                                 aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAATTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAag                     >  1:386938/1‑150 (MQ=255)
                                                                                                                                 aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAag                     >  1:380969/1‑150 (MQ=255)
                                                                                                                                   aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                   <  2:1990628/150‑1 (MQ=255)
                                                                                                                                    aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                  >  2:113555/1‑150 (MQ=255)
                                                                                                                                         aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATTATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAATAGACAGAAAAGATAAAAAAAGAGCCAGCg             >  2:153385/1‑150 (MQ=255)
                                                                                                                                          aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg            >  1:447235/1‑150 (MQ=255)
                                                                                                                                              gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgc        >  1:529282/1‑150 (MQ=255)
                                                                                                                                               cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt       <  1:1734628/150‑1 (MQ=255)
                                                                                                                                                 aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg     <  2:1796967/150‑1 (MQ=255)
                                                                                                                                                    aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATCTTCCGTAACAGGCTGCTCGTTATGGCGGGGTTGCAAATAGATAGCTTGACTAAACAAACAGATAACTAGACACAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc  <  2:324963/150‑1 (MQ=255)
                                                                                                                                                    aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc  >  2:511899/1‑150 (MQ=255)
                                                                                                                                                    aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc  >  2:1249323/1‑150 (MQ=255)
                                                                                                                                                    aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc  >  1:586125/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTC  >  NC_000913/987273‑987570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: