Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,418 A→T 10.1% intergenic (‑436/+167) ompF ← / ← asnS outer membrane porin F/asparagine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4180AT10.1% 386.9 / 23.1 137intergenic (‑436/+167)ompF/asnSouter membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (51/72);  new base T (6/8);  total (57/80)
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.98e-01

Only 100 of 138 total aligned reads displayed.

GATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987271‑987567
                                                                                                                                                   |                                                                                                                                                      
 aTTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATaa                                                                                                                                                     <  1:163217/150‑1 (MQ=255)
  tttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAg                                                                                                                                                    <  1:440169/150‑1 (MQ=255)
       ttATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGc                                                                                                                                               <  2:778635/150‑1 (MQ=255)
        tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                              <  2:153772/150‑1 (MQ=255)
        tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                              <  2:644190/150‑1 (MQ=255)
        tataTTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCg                                                                                                                                              <  2:1154400/150‑1 (MQ=255)
         atatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCgg                                                                                                                                             <  2:52844/150‑1 (MQ=255)
          tatTTAAATTAACTTAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGCCACAAAATCTTCTATCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGATATTAGCTAAGAAAAAACTAATCCGCAAAAAGAAAGCGGa                                                                                                                                            >  2:497160/1‑150 (MQ=255)
          tatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGa                                                                                                                                            >  2:1330684/1‑150 (MQ=255)
           atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAt                                                                                                                                           >  1:1044666/1‑150 (MQ=255)
           atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAt                                                                                                                                           <  2:1539369/150‑1 (MQ=255)
           atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAt                                                                                                                                           >  2:1096847/1‑150 (MQ=255)
           atTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAt                                                                                                                                           >  1:649512/1‑150 (MQ=255)
           atTTAAATTAACTAAACGGAAATTCTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAt                                                                                                                                           >  1:1371248/1‑150 (MQ=255)
            tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                          >  1:1425399/1‑150 (MQ=255)
              tAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAg                                                                                                                                        <  1:1293427/150‑1 (MQ=255)
              tAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAg                                                                                                                                        <  2:111975/150‑1 (MQ=255)
                aaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGtt                                                                                                                                     <  2:37194/150‑1 (MQ=255)
                   tAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttt                                                                                                                                  <  2:1207218/150‑1 (MQ=39)
                     aCTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATGCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgc                                                                                                                                >  2:455204/1‑150 (MQ=255)
                     aCTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgc                                                                                                                                <  2:1559266/150‑1 (MQ=255)
                      cTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGttttttct                                                                                                                               <  1:1570746/150‑4 (MQ=39)
                       tAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgcgc                                                                                                                              >  1:437779/1‑150 (MQ=255)
                           cGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTaa                                                                                                                          >  1:473041/1‑150 (MQ=255)
                             gAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAATTCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAAt                                                                                                                        <  2:351804/150‑1 (MQ=255)
                                 ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                    <  1:313389/150‑1 (MQ=255)
                                 ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                    <  2:1044666/150‑1 (MQ=255)
                                  tttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGc                                                                                                                   <  1:1611579/150‑1 (MQ=255)
                                  tttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGc                                                                                                                   <  2:677704/150‑1 (MQ=255)
                                  tttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGc                                                                                                                   >  1:424398/1‑150 (MQ=255)
                                  tttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATACAGAATTCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATTCGGATTACTTTTTGCGCGTAAATTTTga                                                                                                                   >  1:716259/1‑149 (MQ=255)
                                   ttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCa                                                                                                                  >  2:1457316/1‑150 (MQ=255)
                                    tGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCaa                                                                                                                 <  1:780161/150‑1 (MQ=255)
                                    tGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCaa                                                                                                                 >  2:801328/1‑150 (MQ=255)
                                      tttCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATg                                                                                                               <  2:644487/150‑1 (MQ=255)
                                        tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTCGCAATGTa                                                                                                             <  2:1425399/150‑1 (MQ=255)
                                         cagaaggAAACTTTATCGAACTGGCACAAAATCTTCTTTCAGACATCCAGAATGCAAAAAAATAGTAAGAAAATTCTGCTATTAGCTAACAAAAAACTAAACCTCAATAAGAATGAGTATTAGTTTTTGCGCGTAAATGTTTCAATGTaa                                                                                                            >  2:1159831/6‑150 (MQ=255)
                                         cTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTGGTTTTTGCGCGTAAATGTTGCAATGTaa                                                                                                            <  2:1092966/150‑1 (MQ=255)
                                         cTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTaa                                                                                                            <  1:11104/150‑1 (MQ=255)
                                          tGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAat                                                                                                           >  2:316046/1‑150 (MQ=255)
                                           gATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                          >  1:286996/1‑119 (MQ=255)
                                           gATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                          <  2:286996/119‑1 (MQ=255)
                                            aTGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAatat                                                                                                         <  2:473041/150‑1 (MQ=255)
                                              ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                       >  2:594259/1‑150 (MQ=255)
                                              ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                       >  2:207262/1‑150 (MQ=255)
                                              ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                       <  2:303784/150‑1 (MQ=255)
                                              ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                       >  1:353205/1‑150 (MQ=255)
                                               gAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATTCCAAAAAATATTATTAAAATTTTTATATTAGATAATATAAAACTAATCCGGAATAATAATTTTGATTAGTTTTTGCGCGTAAATGTTGAAATGTataattaa                                                                                                      <  2:691977/150‑8 (MQ=255)
                                                 aaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGt                                                                                                    <  1:72273/150‑1 (MQ=255)
                                                    tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAc                                                                                                 <  1:831459/150‑1 (MQ=255)
                                                       aTCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgcg                                                                                                                               <  1:785816/117‑1 (MQ=255)
                                                       aTCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTgcg                                                                                                                               >  2:785816/1‑117 (MQ=255)
                                                               ggCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAatat                                                                                                         <  1:1440986/131‑1 (MQ=255)
                                                               ggCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGAGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGt                                                                                      >  2:52671/1‑150 (MQ=255)
                                                                     aaaTCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGc                                                                                >  1:1351864/1‑150 (MQ=255)
                                                                         cttcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCgggg                                                                            <  1:1137645/150‑1 (MQ=255)
                                                                              ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                       <  1:1187944/150‑1 (MQ=255)
                                                                              ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCa                                                                       <  1:1029214/150‑1 (MQ=255)
                                                                                cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaaa                                                                     >  2:1030631/1‑150 (MQ=255)
                                                                                cAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaaa                                                                     >  2:649042/1‑150 (MQ=255)
                                                                                    cATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGa                                                                 <  2:1142287/150‑1 (MQ=255)
                                                                                     aTCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGAt                                                                >  1:60906/1‑150 (MQ=255)
                                                                                      tCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAAAGGATGATCGTTATGGAGGGGTTGCAAATAGAtt                                                               >  1:823502/1‑150 (MQ=255)
                                                                                           aaTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                          <  2:1471561/150‑1 (MQ=255)
                                                                                            aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCAGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGa                                                         <  1:369722/150‑1 (MQ=255)
                                                                                             tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGAc                                                        <  2:544067/150‑1 (MQ=255)
                                                                                              gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                       >  2:1105301/1‑150 (MQ=255)
                                                                                                aaaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTTCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTa                                                     <  2:79293/149‑1 (MQ=255)
                                                                                                 aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGc                                                                                                                   <  1:1003190/87‑1 (MQ=255)
                                                                                                  aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaac                                                   <  1:934706/150‑1 (MQ=255)
                                                                                                         gTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGa                                            >  2:848163/1‑150 (MQ=255)
                                                                                                             gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAc                                        <  1:1087092/150‑1 (MQ=255)
                                                                                                             gAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAAc                                        <  1:1259515/150‑1 (MQ=255)
                                                                                                                aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                     <  2:1363799/150‑1 (MQ=255)
                                                                                                                aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                     <  1:1090140/150‑1 (MQ=255)
                                                                                                                aaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTg                                     <  1:562610/150‑1 (MQ=255)
                                                                                                                   tgtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACa                                  <  1:989330/150‑1 (MQ=255)
                                                                                                                            aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATa                         <  1:1540921/150‑1 (MQ=255)
                                                                                                                             gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                        >  1:370918/1‑150 (MQ=255)
                                                                                                                              cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                       >  1:363208/1‑150 (MQ=255)
                                                                                                                                 aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaa                    <  1:15276/150‑1 (MQ=255)
                                                                                                                                    aaaaaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                 >  1:610700/1‑150 (MQ=39)
                                                                                                                                    aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                 <  1:526101/150‑1 (MQ=255)
                                                                                                                                     aaaaCTAATCCTCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                <  1:1247583/150‑1 (MQ=255)
                                                                                                                                      aaaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc               >  1:1296156/1‑150 (MQ=255)
                                                                                                                                      aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc               <  2:197315/150‑1 (MQ=255)
                                                                                                                                       aaCTAATCCGCATTCTTATTGCGGATTAGTTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAg                            >  2:1325100/1‑137 (MQ=38)
                                                                                                                                       aaCTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc              >  1:568662/1‑150 (MQ=255)
                                                                                                                                       aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc              >  1:990062/1‑150 (MQ=255)
                                                                                                                                       aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc              <  1:1224068/150‑1 (MQ=255)
                                                                                                                                         cTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg            >  2:554371/1‑150 (MQ=255)
                                                                                                                                         cTAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAg            >  2:1028317/1‑150 (MQ=255)
                                                                                                                                          tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc           <  1:316046/150‑1 (MQ=255)
                                                                                                                                            aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg         >  2:1458164/1‑150 (MQ=255)
                                                                                                                                             tCCGCAATAAGAATTCGGATTAGTTTTTGCGAGTAAATGTTGCAATTTTATATTCAGTAACAGTATGATCGTTATGGCGGGTTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg        <  1:1159831/150‑1 (MQ=255)
                                                                                                                                               cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg      <  1:207262/150‑1 (MQ=255)
                                                                                                                                               cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg      <  1:1457316/150‑1 (MQ=255)
                                                                                                                                                gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGccgc     >  1:482237/1‑150 (MQ=255)
                                                                                                                                                 cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCt    <  1:801328/150‑1 (MQ=255)
                                                                                                                                                   aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg  <  2:1371248/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987271‑987567

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