Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,044 A→G 6.0% intergenic (‑41/+13) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0440AG 6.0% 712.3 / 24.1 215intergenic (‑41/+13)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (122/80);  new base G (9/4);  total (131/84)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 215 total aligned reads displayed.

CTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCG  >  NC_000913/216898‑217190
                                                                                                                                                  |                                                                                                                                                  
 tACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCt                                                                                                                                                >  2:854999/1‑150 (MQ=255)
 tACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCt                                                                                                                                                >  1:2281576/1‑150 (MQ=255)
   ctgctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCttt                                                                                                                                              >  2:550865/1‑150 (MQ=255)
    tgctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                                                                                             >  1:1027114/1‑150 (MQ=255)
    tgctgTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTc                                                                                                                                             >  2:919956/1‑150 (MQ=255)
       tgtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                                                                          >  2:1182598/1‑150 (MQ=255)
       tgtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                                                                          >  1:1497426/1‑150 (MQ=255)
        gtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                                         <  1:1195764/150‑1 (MQ=255)
        gtgGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                                         >  2:801009/1‑150 (MQ=255)
          ggCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCAt                                                                                                                                       >  1:345232/1‑150 (MQ=255)
              gATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGc                                                                                                                                   >  1:1141062/1‑150 (MQ=255)
               aTGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                                                                  >  1:2654353/1‑150 (MQ=255)
                tGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                                                                 >  2:101043/1‑150 (MQ=255)
                 gggTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCATGAAATTCAGGACCCTTTCTGGCATCAGCCtt                                                                                                                                >  1:2362963/1‑150 (MQ=255)
                       ggctggctGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCt                                                                                                                          <  2:2273347/150‑1 (MQ=255)
                          tggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                                                                       >  2:1770175/1‑150 (MQ=255)
                           ggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                                      <  2:2147387/150‑1 (MQ=255)
                           ggctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                                      <  2:2620331/150‑1 (MQ=255)
                             ctgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCa                                                                                                                    <  1:838597/150‑1 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCttt                                                                                                                               >  2:317550/1‑138 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTacac                                                                                                                   >  1:2088601/1‑150 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                                   >  2:2324405/1‑150 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                                   >  2:1221442/1‑150 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                                   >  2:2453565/1‑150 (MQ=255)
                              tgATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                                                                                   >  1:1452500/1‑150 (MQ=255)
                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                                                                  >  2:2038307/1‑150 (MQ=255)
                               gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                                                                  <  1:534663/150‑1 (MQ=255)
                                     aaCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGatat                                                                                                            <  1:2217372/150‑1 (MQ=255)
                                       cGGTACACGCCGAACGTAAAAGATAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCTGGACCCTTTCTGGCATCAGCCTTTAAGCTGTTTCACCTGATATTc                                                                                                          <  2:1028421/150‑1 (MQ=255)
                                       cGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTc                                                                                                          <  2:1498800/150‑1 (MQ=255)
                                         gTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcga                                                                                                        <  2:185338/150‑1 (MQ=255)
                                          tACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                                       >  2:223146/1‑150 (MQ=255)
                                            cAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                     <  1:773771/150‑1 (MQ=255)
                                                ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTc                                                                                                 >  2:2009996/1‑150 (MQ=255)
                                                ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTc                                                                                                 >  1:2484831/1‑150 (MQ=255)
                                                ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTc                                                                                                 >  2:429135/1‑150 (MQ=255)
                                                 cGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                                                                                >  1:2157115/1‑150 (MQ=255)
                                                    aaGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCa                                                                                             <  1:943819/150‑1 (MQ=255)
                                                     aGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTATCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAg                                                                                            >  1:1280781/1‑150 (MQ=255)
                                                            gaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                                                                     <  1:1200228/150‑1 (MQ=255)
                                                            gaaagCAGCGCACAAACACGCGGCACAACGCCATTGGAATAACCAAAAAAAACAATCCAACTAAAAAAAGCGCCCTCAAATACACCACCCTTTCTGGCAACAGCCTATAATCTGTTTCACCAGATATTCGACGATCTCACCAGTCtaaat                                                                                     <  2:2667586/150‑1 (MQ=255)
                                                              aagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatta                                                                                   >  2:2707726/1‑150 (MQ=255)
                                                               agCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCGGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaattaa                                                                                  <  2:502147/150‑1 (MQ=255)
                                                                  tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                               >  1:1609240/1‑150 (MQ=255)
                                                                  tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                               >  2:308391/1‑150 (MQ=255)
                                                                  tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                               >  2:1691567/1‑150 (MQ=255)
                                                                  tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                                                                               >  1:695075/1‑150 (MQ=255)
                                                                     ggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGttt                                                                            >  1:331399/1‑150 (MQ=255)
                                                                      gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                                                                           >  2:1965713/1‑150 (MQ=255)
                                                                         aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                                                                        >  2:886901/1‑150 (MQ=255)
                                                                         aaaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                                                                        >  1:201000/1‑150 (MQ=255)
                                                                           aCAGGCGGCAGTACGCCTTTTGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACAAGTCTTAATTAACTGTTTCTCGcc                                                                      <  2:1555166/150‑1 (MQ=255)
                                                                           aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGcc                                                                      >  2:2579468/1‑150 (MQ=255)
                                                                           aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGcc                                                                      >  2:134323/1‑150 (MQ=255)
                                                                            cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCg                                                                     >  2:907139/1‑150 (MQ=255)
                                                                             aGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGt                                                                    >  1:310372/1‑150 (MQ=255)
                                                                             aGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTGAATTAACTGTTTCTCGCCGt                                                                    >  2:2180082/1‑150 (MQ=255)
                                                                               gcggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                                                                  >  2:2392431/1‑150 (MQ=255)
                                                                                  gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTCTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga                                                               <  1:1593549/150‑1 (MQ=255)
                                                                                      tACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGat                                                           <  1:2324405/150‑1 (MQ=255)
                                                                                      tACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGat                                                           >  2:2089039/1‑150 (MQ=255)
                                                                                           cTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta                                                      >  1:636914/1‑150 (MQ=255)
                                                                                                gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTc                                                 <  2:725745/150‑1 (MQ=255)
                                                                                                gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTc                                                 >  2:571096/1‑150 (MQ=255)
                                                                                                      ccATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtc                                           <  1:1004449/150‑1 (MQ=255)
                                                                                                       cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcg                                          >  1:1627187/1‑150 (MQ=255)
                                                                                                       cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcg                                          >  2:463417/1‑150 (MQ=255)
                                                                                                         ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcgtc                                        >  1:1908060/1‑150 (MQ=255)
                                                                                                             ataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTg                                    <  1:2240305/150‑1 (MQ=255)
                                                                                                                cactCCATGTAAAAAAAGGGCCCGGAAATTCAGGACCCTTTCTGGCATCAGCCTTTCATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCg                                 <  2:1276302/147‑1 (MQ=255)
                                                                                                                cAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCg                                 <  2:1048479/150‑1 (MQ=255)
                                                                                                                cAATCCATGTAAAAAAAGGGCCCGGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACAATATTTATATTCGATATCGTCGTTGTCg                                 <  1:1342259/150‑1 (MQ=255)
                                                                                                                  aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAg                               <  2:2484831/150‑1 (MQ=255)
                                                                                                                       tGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTAc                          <  2:483998/150‑1 (MQ=255)
                                                                                                                        gTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTaag                         >  2:2536710/1‑148 (MQ=255)
                                                                                                                        gTAAAAAAAGGGCCCTGAAATTCAGGACCCTGTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCTACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAGTTATATGCGATATCGTCGTTGTCGCGGttcag                         >  2:1742818/1‑147 (MQ=255)
                                                                                                                           aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTc                      >  2:2498950/1‑150 (MQ=255)
                                                                                                                           aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTc                      >  1:285489/1‑150 (MQ=255)
                                                                                                                           aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACGGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTc                      <  2:289803/150‑1 (MQ=255)
                                                                                                                            aaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCg                     >  1:996569/1‑150 (MQ=255)
                                                                                                                              aaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGcc                   <  1:1905125/150‑1 (MQ=255)
                                                                                                                                 gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTAACCAGTCTTAATTAACTTTTTCTCGCCGTTGCGACGATATTTATATTCTATATCGTCGTTGTCGAGGTTAAGGTCGCCCAg                >  1:944039/1‑150 (MQ=255)
                                                                                                                                 gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAg                <  1:1815950/150‑1 (MQ=255)
                                                                                                                                  ggtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                     >  1:539264/4‑64 (MQ=37)
                                                                                                                                  ggtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                     <  2:539238/61‑1 (MQ=37)
                                                                                                                                  ggtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                                     <  2:539297/61‑1 (MQ=37)
                                                                                                                                  ggCCCTGAAAGTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGc               <  2:2683283/150‑1 (MQ=255)
                                                                                                                                   gtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGca              >  2:1437360/3‑150 (MQ=255)
                                                                                                                                   gtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTATCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGca              >  1:787154/3‑150 (MQ=255)
                                                                                                                                   gCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGca              >  2:2121532/1‑150 (MQ=255)
                                                                                                                                     ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcaca            >  2:1149967/1‑150 (MQ=255)
                                                                                                                                     ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcaca            >  2:1941311/1‑150 (MQ=255)
                                                                                                                                      cggAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACaa           <  1:216063/148‑1 (MQ=255)
                                                                                                                                      cTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACaa           <  1:1747888/150‑1 (MQ=255)
                                                                                                                                       tGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATACGACGATGTCACCAGTCTTAATTAACAGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcacact          >  2:1664694/1‑148 (MQ=255)
                                                                                                                                           aTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAgtg      >  2:2345095/1‑150 (MQ=255)
                                                                                                                                            ttCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAgtgt     >  2:1417352/1‑150 (MQ=255)
                                                                                                                                              cAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGc   <  1:2618178/150‑1 (MQ=255)
                                                                                                                                               aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCg  >  1:2389548/1‑150 (MQ=255)
                                                                                                                                               aGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCg  >  1:213129/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTACTGCTGTGGCTGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCACAATAGTGTGCG  >  NC_000913/216898‑217190

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