Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,446,969 T→G 17.6% V953G (GTG→GGG)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,446,9690TG17.6% 124.9 / 2.3 70V953G (GTG→GGG) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (22/33);  new base G (9/6);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  NC_000913/4446732‑4447207
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
gggttgtgAAAAGCCCGGTTTTGTTTGGGCAGTTTCCGGGTAGGGGTGGGGATATGGATGGCAACGTTATGCCGTTTGATATGACACCGAGCCAGCTTGCGGTAAACGTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTg                                                                                                                                                                                                                                                  <  2:78651/235‑1 (MQ=255)
      tgCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGGCCGGGTAACGGCGAAAGGCGGGAAAGGCCGGAt                                                                                                                                                                                                                                             >  1:418156/1‑238 (MQ=255)
      tgCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAACGGCGAAAGGCCGTAAAGGGCGGagc                                                                                                                                                                                                                                            >  1:302092/1‑237 (MQ=255)
           aGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGAGCAAAGCGGTAAAGGGCGGGTAAccat                                                                                                                                                                                                                                       >  1:442683/1‑237 (MQ=255)
                gCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGGGAAAGGCCGTAAAAGGCGGATCACCGTgccg                                                                                                                                                                                                                                   >  1:174464/1‑238 (MQ=255)
                       tttGGTGTGTTTGAGGTTGCGGGTTGGGATTTCGAGGGCACATTTATGCCGTTTGAATTGAAACCGTGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGCCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTGCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                           <  2:562946/239‑1 (MQ=255)
                       tttGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGGGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGAGGATCACCGGGCCGCCGAGGGt                                                                                                                                                                                                                           >  1:489619/1‑239 (MQ=255)
                            tCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAACGGCGAAAGGCAGTAAAGTGGGGGTCACCGTGCCGCCGATGGTACGaa                                                                                                                                                                                                                      >  1:849588/1‑239 (MQ=255)
                              aGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCCATGGGACGAA‑Tg                                                                                                                                                                                                                   >  1:628201/1‑239 (MQ=255)
                                cTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGAAAGGCAGTAAAGTGCGGGACACCGTGCCGCCGATGGGACGGA‑TGGa                                                                                                                                                                                                                 >  1:144187/1‑239 (MQ=255)
                                        tgcGGGTGGGGATAACGGTGGCAACATTATGCCGTTTGATATGCAACCGGGCCAACTTGCGGTCAACTTTAACGGGTGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCg                                                                                                                                                                                                          <  1:73983/236‑1 (MQ=255)
                                          cGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGATCACCGGGCCGCCGATGGTACGAA‑TGGGTGTATCGCCa                                                                                                                                                                                                       >  1:452702/1‑239 (MQ=255)
                                             tggtgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGGGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGa                                                                                                                                                                                                     <  1:488690/236‑1 (MQ=255)
                                                tgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGGTCACCGTGCCGCCGGTGGTCCGAA‑TGGTGGTATCGCCCGAtgtt                                                                                                                                                                                                 >  1:74788/1‑239 (MQ=255)
                                                 ggATATCGATGTAAACTTTATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCGCGACGCTGTCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGAtgtt                                                                                                                                                                                                 <  1:506134/238‑1 (MQ=255)
                                                  gttttctttGGCAACTTTATGCCGTGTGATATGCACCCGGGCCAGGTTGGTGTCAACTTTAACGGTAGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGAtgttg                                                                                                                                                                                                <  2:630951/230‑1 (MQ=255)
                                                           ggaaattttATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGc                                                                                                                                                                                     <  1:626674/233‑1 (MQ=255)
                                                           gccaccTTTATGCCGTTTGATTTGAAACCGAGCCAACTTGCGGCAAACTTTAACGGTATGCCCCGGACGCTTGCCGGTAAAGTCCGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGc                                                                                                                                                                                     <  2:235225/234‑1 (MQ=255)
                                                                tttttgtccGTTTGTTAGGAACCCGAGCCAGCTTGCGGTCACTTTGAACGGTATGCGCTCGGCGGTTGCCGGTACAGTACGGACCCAGAAGGGTGAAATCTCCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTacac                                                                                                                                                                                <  2:418156/232‑1 (MQ=255)
                                                                cTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGc                                                                                                                                                                                     <  1:592878/234‑1 (MQ=255)
                                                                   tATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCaa                                                                                                                                                                             <  1:233425/239‑1 (MQ=255)
                                                                          tttggattggaaCCGGGCCAGGTTGGGGACAACTTTAAGGGTATGGGCTCGATGCTGGCCGGTACAGTACGGACCCAGGAGGGTGAAAACTACCTGAACGGTGATGCGGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGtt                                                                                                                                                                      <  2:634049/229‑1 (MQ=255)
                                                                           ttGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCCCCGGTAAAGGCGAAAGACAGAAAAAGGCGGCACCCCGTGCCGCCGGAGGCACGAAGTGGGT‑TATCCCCACATGTTAGT‑TTCGGGCCTAAACAAAAcacattt                                                                                                                                                                     >  2:248687/1‑233 (MQ=255)
                                                                                     cccGGGCCAGCTTGCTGTCCAGTTTACCGGTAGGCGCCCGACGCTTGCCGGGCAAGTGCGGACCCGGCAGGGTGAAGTGTGCCTGAACGGTGATGCCGACTGGAACCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATg                                                                                                                                                           <  2:63565/238‑1 (MQ=255)
                                                                                      ccGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGGGATGCCGACTGGAGCCAACTTGAAAACTGGCGGGCGCGGGGAACGGCGAAACGCGGGAAAGTGCGGATTCCCGTGCCGCCGGTGGGACGAA‑CGGGTGTATAGCCCGCTGTTGTA‑TTTGAGGCACTACCAAAACTGTATAaccccccggg                                                                                                                                                          >  2:594265/1‑232 (MQ=255)
                                                                                       cggccccttttGCGTTCAACGTATACTTTAGGCGCTAGACGCTTGCCGGTCCAGTACGGACCCAGCGGGGTAAAATCTCCCCGAACGGGGAGGCCGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGt                                                                                                                                                         <  2:470295/230‑1 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAACCCTGTTTACCCTCCATGGTc                                                                                                                                                        >  1:322604/1‑238 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                       >  1:428502/1‑239 (MQ=255)
                                                                                          gCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                       <  1:179860/238‑1 (MQ=255)
                                                                                              gCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGCAAGGCAGGAAAAGGAGGGTAACCCTGGCGCCCGTGGTCCGCA‑AGGATGTCTCGCCCAATAGTGGA‑TTCGAGGCGACACCCAAACGCTTAACCCCCGAGGGGcgggggg                                                                                                                                                  >  2:605181/1‑234 (MQ=255)
                                                                                                ttGCGGTCAACTGTAACGGTATGCTCCCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGa                                                                                                                                                <  1:453220/238‑1 (MQ=255)
                                                                                                tggCGGTCAACTTTAACGGTATGCCCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGa                                                                                                                                                <  1:196247/236‑1 (MQ=255)
                                                                                                  gCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGAGCACCGGGCCGCCGGTGGTACGAA‑TGGATGTATCGCCAGATGTTGGA‑TTCGAGGCTACACCAAACCCGTTTACCCTCGTTGGTCGCGT‑GGATg                                                                                                                                              >  1:761240/1‑238 (MQ=255)
                                                                                                           tttAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGTCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGgc                                                                                                                                    <  1:603059/239‑1 (MQ=255)
                                                                                                           tttAACGGTAGGCGCTCGACGCTTGCCGGGACAGGACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACGGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGgc                                                                                                                                    <  2:404145/239‑1 (MQ=255)
                                                                                                                cGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGAACACCCTGCCGCCCGTGGGACGAA‑ATGATGATTCGCCCGAGGTGGTA‑TTCGCGGCTACACCAAAACTGATTAACTCCGTGGGGCGCGTGGGATTTCCTTGGGgcggcc                                                                                                                                >  2:215863/1‑236 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCg                                                                                                                            <  1:515984/238‑1 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCCCCATGGTCGCGT‑GGGTGTCCCTTGGGGGCGCCTCgt                                                                                                                           >  1:540032/1‑239 (MQ=255)
                                                                                                                     tGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGGTCACCCGGCCGCCGAGGGTACGAA‑AGGATGTATCGCCAAATTTTGTATTTCAGGGCTA‑ACCCAAACTGTTTTACCTCCGAGGTCGCTG‑GGGTGTCCCCTggggcccgccaggg                                                                                                                          >  2:515984/1‑228 (MQ=255)
                                                                                                                          tCGACGCTTGCCGGGACAGTACGGACCCAGCAGGGTGAACTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCa                                                                                                                     <  1:82886/239‑1 (MQ=255)
                                                                                                                            gcccgttGCCGGTACAGGACGGACCCGGCAGGGTGAAATCTACCTGAACGGTGATGGCGACTGGAGGCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑GGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACg                                                                                                                   <  1:681904/234‑1 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCCGATGTTGTA‑TTCGAGGCTACACCAAACCCGTTTACCCTCGATGGTCTCGT‑GGATGTCCCTTTGGGGCGCATCGTGGTGCACGa                                                                                                                  >  1:49758/1‑239 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGa                                                                                                                  >  1:871576/1‑239 (MQ=255)
                                                                                                                                         aCAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGGGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAg                                                                                                      <  1:104004/239‑1 (MQ=255)
                                                                                                                                                ggACCCCGTGGGGTTAAATCTACCTGAACGGGGATGCCGATGGGGGCCAAATTGAAAACTGGGGGGCGGGGGTAACGGGGAAAGGGAGTAAAGGGCGGATCGCCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTa                                                                                               <  2:174464/239‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                           <  2:210564/116‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                           >  1:210564/1‑116 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAgg                                                                                             >  1:708490/1‑239 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGGGACCGATCTGCCGGAAAGCCCCGTAg                                                                                              >  1:710205/1‑238 (MQ=255)
                                                                                                                                                      agcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGTGCCGCCGATGGGACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGGGCACCGTGGGGCACGATCTCCCGGGAAGCGCAGATGGCgc                                                                                          >  1:592720/1‑237 (MQ=255)
                                                                                                                                                       gcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGGCT‑GGAGGTCCCTTGGGCGCGCATCGTGGTGCCCGAATTGCCGGAAAGCGCAGTAGGCGtc                                                                                         >  1:455617/1‑238 (MQ=255)
                                                                                                                                                          gggTAAAATCTACCCGAACGGTGATGCCGGCTGGAACCAAATTGAAAACTGGGGGGCGCGGGTAACTGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTcc                                                                                      <  1:276817/238‑1 (MQ=255)
                                                                                                                                                           ggTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGa                                                                                                                                                <  1:852457/179‑1 (MQ=255)
                                                                                                                                                              gAAGTCTACCGGAACGGGGATGCCGACGGGAGACAAATTGAAAACTGGCGGGGGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGa                                                                                 <  1:441036/239‑1 (MQ=255)
                                                                                                                                                                aaTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGAtg                                                                               <  1:170233/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAACCTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  1:374012/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAAAGGTGAGGCCGAAGGGGGCCAAATTGAAAACGGGCGGGCCGGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  1:324653/239‑1 (MQ=255)
                                                                                                                                                                                      ggatggtggcAAAGTGAAAACGGGGAGGGGGGGGTAACGGGGAAAAGGAGTAAAGGGCGGCTCCCGGTGCCGACGGTGGTACGAA‑TGGAATGAACGGCAGATGTTGGA‑TTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCTATGTGGTGATGCTTAACGATAACCt                                                         <  2:733073/230‑1 (MQ=255)
                                                                                                                                                                                        agggAGCAAAATTGAAAACTGGCGGGCGCGGGTAACGGCAAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTg                                                        <  1:605181/236‑1 (MQ=255)
                                                                                                                                                                                         ttgggCCAAATTGAAAATGGGCGGGGGCGGGTAACGGCGAAAGGGAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAA‑TGGATGGTTCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGc                                                       <  2:761240/234‑1 (MQ=255)
                                                                                                                                                                                         tGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGAAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGC‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCCGCGCTGTGGTGATGCTTAACGGTAAACTGCa                                                      >  1:303416/1‑239 (MQ=255)
                                                                                                                                                                                         gggtgCCAATTTAAAAACGGGCGGGGGCGGGTAAGGGCGAAAGGCGGTAAGGTGCGGATCCCCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                      <  2:708490/235‑1 (MQ=255)
                                                                                                                                                                                         ggaaGCCAAATTGAAAACTGGCGGGAGGGGGAAACGGCGGAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAA‑TGGATGTAACGCCAGATGTTGTA‑TTCGAGGCTTCACCAAACCTGTTTAACCTCCATGGTCGCGG‑GGATGTGCCTTTGGCGCGCACCGTGGGGCCCGAATTGCCGGGAAACCCAGTGAGGGACTCCCCCGAAGCGGGAATGCTTAAAGAAAACCCGcc                                                      >  2:39761/4‑238 (MQ=255)
                                                                                                                                                                                                 aaTTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGGGGGGATCAACGTGCCGCCGAGAGTACGAA‑TGGATGTATCGGCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGGTGGTAGCGT‑GGATGTTCCCTGGGAGCGCATAGTGGTGCCCGATCTGCCCGATAACGCAGTAGGCGTCCCCACGGACGTGGTGATACTTAACGAAAACCAGCAACCGaaaa                                              >  2:352586/1‑238 (MQ=255)
                                                                                                                                                                                                             gcgggcgCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGC‑GGATGTCCCTTGGGCGCGCATCGTGGGGGACGGTCTGCCGGAAAGGGCAGTAGGCGTCTCCCGCGATGTGGGGGTGCTTAACGAGAACCCCCAACCGGGAGAAGCGAAAAAgg                                  >  2:485532/1‑239 (MQ=255)
                                                                                                                                                                                                                          aCGGCGAAAGGCAGTAAAGTGCGGATCCCCGGGCCGCCGAGGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGAt                     <  1:709207/239‑1 (MQ=255)
                                                                                                                                                                                                                               aaaggggtttAAAGGGGGGTTCACCGGGCCGCCGATGGTACGAA‑TGGATTTATCGCCAGGTGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAAc                 <  2:452702/229‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGCATTAAAGGGGGGTTCCCCGTGCCGCTGATGGTACAAA‑TGGATGTATCCGCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg               <  2:144187/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGGGATGCTTAAAGAGAACCTGCAACCGGAAGAGCCGAAAAAGGCGGGGATTCCGTTTAAAAGaaaacgga       >  2:743183/1‑234 (MQ=255)
                                                                                                                                                                                                                                         aaaGTGCGGATCCCCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGa       <  1:39761/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGt   >  1:758274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGGGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGCCGTTTCCGATTAAAAGAAACCTGATTgac  >  2:33437/1‑237 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCCCCCGCGATGTGGTGTTGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCTTCAATTCCGATTAAAAGAAAACTAATTGt   >  2:416274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAA‑TGGATGTATCGCCAGATGTTGTA‑TTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGT‑GGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  NC_000913/4446732‑4447207

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑