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

GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  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)
                              aGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCCATGGGACGAATg                                                                                                                                                                                                                   >  1:628201/1‑239 (MQ=255)
                                cTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGAAAGGCAGTAAAGTGCGGGACACCGTGCCGCCGATGGGACGGATGGa                                                                                                                                                                                                                 >  1:144187/1‑239 (MQ=255)
                                        tgcGGGTGGGGATAACGGTGGCAACATTATGCCGTTTGATATGCAACCGGGCCAACTTGCGGTCAACTTTAACGGGTGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCg                                                                                                                                                                                                          <  1:73983/236‑1 (MQ=255)
                                          cGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGATCACCGGGCCGCCGATGGTACGAATGGGTGTATCGCCa                                                                                                                                                                                                       >  1:452702/1‑239 (MQ=255)
                                             tggtgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGGGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGa                                                                                                                                                                                                     <  1:488690/236‑1 (MQ=255)
                                                tgGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGGTCACCGTGCCGCCGGTGGTCCGAATGGTGGTATCGCCCGAtgtt                                                                                                                                                                                                 >  1:74788/1‑239 (MQ=255)
                                                 ggATATCGATGTAAACTTTATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCGCGACGCTGTCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGAtgtt                                                                                                                                                                                                 <  1:506134/238‑1 (MQ=255)
                                                  gttttctttGGCAACTTTATGCCGTGTGATATGCACCCGGGCCAGGTTGGTGTCAACTTTAACGGTAGGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGAtgttg                                                                                                                                                                                                <  2:630951/230‑1 (MQ=255)
                                                           ggaaattttATGCCGTTTGATATGCAACCGAGCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                     <  1:626674/233‑1 (MQ=255)
                                                           gccaccTTTATGCCGTTTGATTTGAAACCGAGCCAACTTGCGGCAAACTTTAACGGTATGCCCCGGACGCTTGCCGGTAAAGTCCGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                     <  2:235225/234‑1 (MQ=255)
                                                                tttttgtccGTTTGTTAGGAACCCGAGCCAGCTTGCGGTCACTTTGAACGGTATGCGCTCGGCGGTTGCCGGTACAGTACGGACCCAGAAGGGTGAAATCTCCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTacac                                                                                                                                                                                <  2:418156/232‑1 (MQ=255)
                                                                cTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGc                                                                                                                                                                                     <  1:592878/234‑1 (MQ=255)
                                                                   tATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCaa                                                                                                                                                                             <  1:233425/239‑1 (MQ=255)
                                                                          tttggattggaaCCGGGCCAGGTTGGGGACAACTTTAAGGGTATGGGCTCGATGCTGGCCGGTACAGTACGGACCCAGGAGGGTGAAAACTACCTGAACGGTGATGCGGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGtt                                                                                                                                                                      <  2:634049/229‑1 (MQ=255)
                                                                           ttGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCCCCGGTAAAGGCGAAAGACAGAAAAAGGCGGCACCCCGTGCCGCCGGAGGCACGAAGTGGGTTATCCCCACATGTTAGTTTCGGGCCTAAACAAAAcacattt                                                                                                                                                                     >  2:248687/1‑233 (MQ=255)
                                                                                     cccGGGCCAGCTTGCTGTCCAGTTTACCGGTAGGCGCCCGACGCTTGCCGGGCAAGTGCGGACCCGGCAGGGTGAAGTGTGCCTGAACGGTGATGCCGACTGGAACCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATg                                                                                                                                                           <  2:63565/238‑1 (MQ=255)
                                                                                      ccGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGGGATGCCGACTGGAGCCAACTTGAAAACTGGCGGGCGCGGGGAACGGCGAAACGCGGGAAAGTGCGGATTCCCGTGCCGCCGGTGGGACGAACGGGTGTATAGCCCGCTGTTGTATTTGAGGCACTACCAAAACTGTATAaccccccggg                                                                                                                                                          >  2:594265/1‑232 (MQ=255)
                                                                                       cggccccttttGCGTTCAACGTATACTTTAGGCGCTAGACGCTTGCCGGTCCAGTACGGACCCAGCGGGGTAAAATCTCCCCGAACGGGGAGGCCGACTGGAGCCAAATTGAAAACGGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGt                                                                                                                                                         <  2:470295/230‑1 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAACCCTGTTTACCCTCCATGGTc                                                                                                                                                        >  1:322604/1‑238 (MQ=255)
                                                                                         aGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                       >  1:428502/1‑239 (MQ=255)
                                                                                          gCCAGGTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTcg                                                                                                                                                       <  1:179860/238‑1 (MQ=255)
                                                                                              gCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGGAAAGGCGCAAGGCAGGAAAAGGAGGGTAACCCTGGCGCCCGTGGTCCGCAAGGATGTCTCGCCCAATAGTGGATTCGAGGCGACACCCAAACGCTTAACCCCCGAGGGGcgggggg                                                                                                                                                  >  2:605181/1‑234 (MQ=255)
                                                                                                ttGCGGTCAACTGTAACGGTATGCTCCCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                <  1:453220/238‑1 (MQ=255)
                                                                                                tggCGGTCAACTTTAACGGTATGCCCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                <  1:196247/236‑1 (MQ=255)
                                                                                                  gCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGAGCACCGGGCCGCCGGTGGTACGAATGGATGTATCGCCAGATGTTGGATTCGAGGCTACACCAAACCCGTTTACCCTCGTTGGTCGCGTGGATg                                                                                                                                              >  1:761240/1‑238 (MQ=255)
                                                                                                           tttAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGTCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGgc                                                                                                                                    <  1:603059/239‑1 (MQ=255)
                                                                                                           tttAACGGTAGGCGCTCGACGCTTGCCGGGACAGGACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACGGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGgc                                                                                                                                    <  2:404145/239‑1 (MQ=255)
                                                                                                                cGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAATGCGGAACACCCTGCCGCCCGTGGGACGAAATGATGATTCGCCCGAGGTGGTATTCGCGGCTACACCAAAACTGATTAACTCCGTGGGGCGCGTGGGATTTCCTTGGGgcggcc                                                                                                                                >  2:215863/1‑236 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCg                                                                                                                            <  1:515984/238‑1 (MQ=255)
                                                                                                                    aTGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCCCCATGGTCGCGTGGGTGTCCCTTGGGGGCGCCTCgt                                                                                                                           >  1:540032/1‑239 (MQ=255)
                                                                                                                     tGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGGTCACCCGGCCGCCGAGGGTACGAAAGGATGTATCGCCAAATTTTGTATTTCAGGGCTAACCCAAACTGTTTTACCTCCGAGGTCGCTGGGGTGTCCCCTggggcccgccaggg                                                                                                                          >  2:515984/1‑228 (MQ=255)
                                                                                                                          tCGACGCTTGCCGGGACAGTACGGACCCAGCAGGGTGAACTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCa                                                                                                                     <  1:82886/239‑1 (MQ=255)
                                                                                                                            gcccgttGCCGGTACAGGACGGACCCGGCAGGGTGAAATCTACCTGAACGGTGATGGCGACTGGAGGCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAAGGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACg                                                                                                                   <  1:681904/234‑1 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCCGATGTTGTATTCGAGGCTACACCAAACCCGTTTACCCTCGATGGTCTCGTGGATGTCCCTTTGGGGCGCATCGTGGTGCACGa                                                                                                                  >  1:49758/1‑239 (MQ=255)
                                                                                                                             aCGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGa                                                                                                                  >  1:871576/1‑239 (MQ=255)
                                                                                                                                         aCAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGGGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAg                                                                                                      <  1:104004/239‑1 (MQ=255)
                                                                                                                                                ggACCCCGTGGGGTTAAATCTACCTGAACGGGGATGCCGATGGGGGCCAAATTGAAAACTGGGGGGCGGGGGTAACGGGGAAAGGGAGTAAAGGGCGGATCGCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTa                                                                                               <  2:174464/239‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                           <  2:210564/116‑1 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGt                                                                                                                                                                                                                           >  1:210564/1‑116 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAgg                                                                                             >  1:708490/1‑239 (MQ=255)
                                                                                                                                                  aCCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGGGACCGATCTGCCGGAAAGCCCCGTAg                                                                                              >  1:710205/1‑238 (MQ=255)
                                                                                                                                                      agcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGGGCACCGTGGGGCACGATCTCCCGGGAAGCGCAGATGGCgc                                                                                          >  1:592720/1‑237 (MQ=255)
                                                                                                                                                       gcagGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGGCTGGAGGTCCCTTGGGCGCGCATCGTGGTGCCCGAATTGCCGGAAAGCGCAGTAGGCGtc                                                                                         >  1:455617/1‑238 (MQ=255)
                                                                                                                                                          gggTAAAATCTACCCGAACGGTGATGCCGGCTGGAACCAAATTGAAAACTGGGGGGCGCGGGTAACTGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTcc                                                                                      <  1:276817/238‑1 (MQ=255)
                                                                                                                                                           ggTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGa                                                                                                                                                <  1:852457/179‑1 (MQ=255)
                                                                                                                                                              gAAGTCTACCGGAACGGGGATGCCGACGGGAGACAAATTGAAAACTGGCGGGGGCGGGTAACGGCGAAAGGCAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGa                                                                                 <  1:441036/239‑1 (MQ=255)
                                                                                                                                                                aaTCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGAtg                                                                               <  1:170233/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAACGGTGATGCCGACTGGAGCCAAATTGAAACCTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  1:374012/239‑1 (MQ=255)
                                                                                                                                                                     aCCTGAAAGGTGAGGCCGAAGGGGGCCAAATTGAAAACGGGCGGGCCGGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATgtggtg                                                                          <  1:324653/239‑1 (MQ=255)
                                                                                                                                                                                      ggatggtggcAAAGTGAAAACGGGGAGGGGGGGGTAACGGGGAAAAGGAGTAAAGGGCGGCTCCCGGTGCCGACGGTGGTACGAATGGAATGAACGGCAGATGTTGGATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCTATGTGGTGATGCTTAACGATAACCt                                                         <  2:733073/230‑1 (MQ=255)
                                                                                                                                                                                        agggAGCAAAATTGAAAACTGGCGGGCGCGGGTAACGGCAAAAGGCAGTAAAGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTg                                                        <  1:605181/236‑1 (MQ=255)
                                                                                                                                                                                         ttgggCCAAATTGAAAATGGGCGGGGGCGGGTAACGGCGAAAGGGAGTAAAGGGCGGATCACCGGGCCGCCGATGGTACGAATGGATGGTTCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGc                                                       <  2:761240/234‑1 (MQ=255)
                                                                                                                                                                                         tGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGAAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGCGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCCGCGCTGTGGTGATGCTTAACGGTAAACTGCa                                                      >  1:303416/1‑239 (MQ=255)
                                                                                                                                                                                         gggtgCCAATTTAAAAACGGGCGGGGGCGGGTAAGGGCGAAAGGCGGTAAGGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCa                                                      <  2:708490/235‑1 (MQ=255)
                                                                                                                                                                                         ggaaGCCAAATTGAAAACTGGCGGGAGGGGGAAACGGCGGAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGGACGAATGGATGTAACGCCAGATGTTGTATTCGAGGCTTCACCAAACCTGTTTAACCTCCATGGTCGCGGGGATGTGCCTTTGGCGCGCACCGTGGGGCCCGAATTGCCGGGAAACCCAGTGAGGGACTCCCCCGAAGCGGGAATGCTTAAAGAAAACCCGcc                                                      >  2:39761/4‑238 (MQ=255)
                                                                                                                                                                                                 aaTTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGGAAAGGGGGGATCAACGTGCCGCCGAGAGTACGAATGGATGTATCGGCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGGTGGTAGCGTGGATGTTCCCTGGGAGCGCATAGTGGTGCCCGATCTGCCCGATAACGCAGTAGGCGTCCCCACGGACGTGGTGATACTTAACGAAAACCAGCAACCGaaaa                                              >  2:352586/1‑238 (MQ=255)
                                                                                                                                                                                                             gcgggcgCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGCGGATGTCCCTTGGGCGCGCATCGTGGGGGACGGTCTGCCGGAAAGGGCAGTAGGCGTCTCCCGCGATGTGGGGGTGCTTAACGAGAACCCCCAACCGGGAGAAGCGAAAAAgg                                  >  2:485532/1‑239 (MQ=255)
                                                                                                                                                                                                                          aCGGCGAAAGGCAGTAAAGTGCGGATCCCCGGGCCGCCGAGGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGAt                     <  1:709207/239‑1 (MQ=255)
                                                                                                                                                                                                                               aaaggggtttAAAGGGGGGTTCACCGGGCCGCCGATGGTACGAATGGATTTATCGCCAGGTGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAAc                 <  2:452702/229‑1 (MQ=255)
                                                                                                                                                                                                                                 aaGGCATTAAAGGGGGGTTCCCCGTGCCGCTGATGGTACAAATGGATGTATCCGCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAg               <  2:144187/238‑1 (MQ=255)
                                                                                                                                                                                                                                        tAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGGGATGCTTAAAGAGAACCTGCAACCGGAAGAGCCGAAAAAGGCGGGGATTCCGTTTAAAAGaaaacgga       >  2:743183/1‑234 (MQ=255)
                                                                                                                                                                                                                                         aaaGTGCGGATCCCCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCGCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGa       <  1:39761/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGt   >  1:758274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGGGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGCCGTTTCCGATTAAAAGAAACCTGATTgac  >  2:33437/1‑237 (MQ=255)
                                                                                                                                                                                                                                             tGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCCCCCGCGATGTGGTGTTGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCTTCAATTCCGATTAAAAGAAAACTAATTGt   >  2:416274/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
GTGATGTGCAAAGCCCGCAGTTGTTTGGTCAGCTTCAGGTTACGGGTGTGGATATCGATGGCAACTTTATGCCGTTTGATATGCAACCGAGCCAGCTTGCGGTCAACTTTAACGGTATGCGCTCGACGCTTGCCGGTACAGTACGGACCCAGCAGGGTGAAATCTACCTGAACGGTGATGCCGACTGGAGCCAAATTGAAAACTGGCGGGCGCGGGTAACGGCGAAAGGCAGTAAAGTGCGGATCACCGTGCCGCCGATGGTACGAATGGATGTATCGCCAGATGTTGTATTCGAGGCTACACCAAACCTGTTTACCCTCGATGGTCGCGTGGATGTCCCTTGGGCGCGCATCGTGGTGCACGATCTGCCGGAAAGCGCAGTAGGCGTCTCCAGCGATGTGGTGATGCTTAACGATAACCTGCAACCGGAAGAGCCGAAAACGGCGTCGATTCCGATTAACAGTAACCTGATTGTC  >  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: