Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,714 C→T G184G (GGC→GGT dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,7140CT95.6% 192.3 / NA 70G184G (GGC→GGTdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/31);  minor base A (3/0);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCA  >  NC_000913/12472‑12960
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCa                                                                                                                                                                                                                                                 >  1:102943/1‑250 (MQ=255)
 aagtaaaagcgaaaaaaaGGGCCCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAATCCGCTGAAGATTACCTGGGGGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGGCAGGCCACCAAAGACGCAGGCCGTATCTCTGGTCTGGAAGTAAAACGTATCATCAACGATCCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGTCAAc                                                                                                                                                                                                                                               <  2:445829/237‑1 (MQ=255)
 aaGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAAc                                                                                                                                                                                                                                               >  1:557460/1‑251 (MQ=255)
    ttAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGt                                                                                                                                                                                                                                            >  1:55209/1‑251 (MQ=255)
     aaaaGGCAAAAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTa                                                                                                                                                                                                                                           <  2:263312/250‑1 (MQ=255)
        aGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACt                                                                                                                                                                                                                                         >  1:5110/1‑250 (MQ=255)
           ccAGAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATcg                                                                                                                                                                                                                                     <  1:204849/251‑1 (MQ=255)
                   aggCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAAACGTACTATCGCGGtttct                                                                                                                                                                                                                             >  2:166672/2‑249 (MQ=255)
                    ggCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATg                                                                                                                                                                                                                            >  1:349468/1‑251 (MQ=255)
                       aCCGCCGCAGATCTCTGCTGAAGTGCTGCAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAATCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGACTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACCAAGGTACTGGCAACCGTACTATCGCGGTTTATGAcc                                                                                                                                                                                                                         >  2:482719/1‑251 (MQ=255)
                                     atgctgAAGTGCTGAAAAAAATGAAGAAAACCGCGGAAGATTACCTGGGTGAACCGGTAACTGAAGTTGTTATCACCGTACCGGCATACTTTAACGATGCTCATCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTAc                                                                                                                                                                                                            <  1:363655/249‑1 (MQ=255)
                                      tgctgAAGGGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACtt                                                                                                                                                                                                          >  1:249344/1‑251 (MQ=255)
                                         tgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGCACCGGTAACTGAAGCTGTTATCACCGTACCGGAATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                       >  1:62050/1‑251 (MQ=255)
                                          gAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                       >  1:459329/1‑250 (MQ=255)
                                          gAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTAACTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTCTATGACCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                       <  2:459329/250‑1 (MQ=255)
                                           aaGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGat                                                                                                                                                                                                     >  2:198382/1‑251 (MQ=255)
                                           aaGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGat                                                                                                                                                                                                     <  1:572300/251‑1 (MQ=255)
                                             gTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCTACCGCAGCTGCGCTGGCTTACGGGCTGGACAAAGGGACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGatat                                                                                                                                                                                                   <  1:36059/251‑1 (MQ=255)
                                               gCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTTCCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCTCAGCTGCTCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGTGTGGTGGGACTTTCTATAtt                                                                                                                                                                                                  <  2:438960/250‑1 (MQ=255)
                                                  gAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTGAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCta                                                                                                                                                                                              <  2:5110/251‑1 (MQ=255)
                                                      aaaaTGAAGAAAACAGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGCTATTTCtattat                                                                                                                                                                                          >  2:116228/1‑251 (MQ=255)
                                                         aTGAAGAAAACCGCTCAAGATTACCGGGGTGAACCGGTAGCTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAATAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACTGTCTGGACAAAGGTTCTGGCAACCGTACTTGCGCGGTTTATGACCTGGGTGGTTTTACTTTCGATATCTCTATTATCGa                                                                                                                                                                                       <  2:312717/251‑1 (MQ=255)
                                                                 aaCCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGTTCTTTCTATTATCGAAATcga                                                                                                                                                                                 >  2:352166/1‑249 (MQ=255)
                                                                       tGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTg                                                                                                                                                                         >  1:53247/1‑251 (MQ=255)
                                                                       tGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTg                                                                                                                                                                         >  2:170743/1‑251 (MQ=255)
                                                                                    ggTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAAc                                                                                                                                                            <  1:548879/251‑1 (MQ=255)
                                                                                     gTGAACCGGTAACTGAAGCTGTTAGCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTAGCGAAATCGACGAAGTTGACGGCGAAAAAAcc                                                                                                                                                           <  2:417011/251‑1 (MQ=255)
                                                                                      tGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAACGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGACAAAACCt                                                                                                                                                          >  1:563709/1‑251 (MQ=255)
                                                                                          ccGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCg                                                                                                                                                       >  1:300814/1‑250 (MQ=255)
                                                                                             gTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGt                                                                                                                                                   <  1:331672/251‑1 (MQ=255)
                                                                                                     aGCTGTTATAACCGAACCGGCATACTTTAACGATAAGCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGACTGGAAGAAAAACGTATCATCAACGTACAGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACAGTAATATCGCGGTTTATGACCTGGGGGGTGGTACTTTAGATATTTATATTATCGAAATCGACGAAGTTGGCGGCAAACAAACATTCGAAGTTCTGGcac                                                                                                                                           >  2:412604/1‑250 (MQ=255)
                                                                                                         gTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTACCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGTACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATAGTTCTATTATCGAAATCGACGAAGTTGACGGCTAATAAGCCTTCGAAGGTCTTGcaacc                                                                                                                                         <  2:243176/249‑1 (MQ=255)
                                                                                                                ccGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGAt                                                                                                                                >  1:80329/1‑251 (MQ=255)
                                                                                                                    aCCGGCATACTTTACCGATGCTCAGCGTCCGGCGACCACCGACGCCGGCCGTATCGCTGGTCTGGAAGTAAAACGAATCATCAACGAACTGACCGCGGCTGGGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACATGGGTCGGGGAACTTTCGATATTTCTATTATCGCAATCGCCGAAGTGGACGGCGTAAAATAATTCGAAGGTCTGGaaacgagaggtgttaccc                                                                                                                            >  2:86622/1‑237 (MQ=255)
                                                                                                                    aCCGGCATACTTTAACGAGGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAccc                                                                                                                            <  2:487127/251‑1 (MQ=255)
                                                                                                                         cATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTg                                                                                                                       >  1:57527/1‑251 (MQ=255)
                                                                                                                         cATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTg                                                                                                                       <  1:231090/251‑1 (MQ=255)
                                                                                                                            aCTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTGGACGGCGAACACACCTTCGAAGTTCTGGCAACCAACGGTGCTACCCACCTgggt                                                                                                                    >  2:535891/1‑250 (MQ=255)
                                                                                                                                 aaCGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGgtggtg                                                                                                                >  2:482590/1‑250 (MQ=255)
                                                                                                                                     gtgcacAGCGTCAGGCCACCAAAGATGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACTGTTTGGACAAAGGTACTGTCAACCGTACTAGCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGAc                                                                                                           <  2:23396/246‑1 (MQ=255)
                                                                                                                                       gCTCAGAGTCAGGCAACCAAAGACGCAGACCGTATCGCTGGTCTGGAAGTAAAACGTATCATAAACGAACCGACCGCAGATGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGGACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATCCCCACCTGGGTGGGGAAGACtt                                                                                                         >  2:550981/1‑251 (MQ=255)
                                                                                                                                          cAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCg                                                                                                       >  2:543804/1‑250 (MQ=255)
                                                                                                                                              gTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGc                                                                                                  >  2:41494/1‑251 (MQ=255)
                                                                                                                                                     aaCCAAAGACGCAGGCCGTATCGCTGGTCTTGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACATAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATCAGCGAAATCGACGAAGTTGACGGCGAAAAAACCTTc                                                                                                                                                        <  2:445118/190‑1 (MQ=255)
                                                                                                                                                     aaCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGGACTGGCACCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGAGGGCGAAAAAACCTTc                                                                                                                                                        >  1:445118/1‑190 (MQ=255)
                                                                                                                                                      aCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGAt                                                                                          >  1:390946/1‑251 (MQ=255)
                                                                                                                                                      aCCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGAAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGAACGCCGTCTGAt                                                                                          >  1:250862/1‑251 (MQ=255)
                                                                                                                                                         aaaGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGggt                                                                                                                                                                                                                    >  2:74263/1‑126 (MQ=255)
                                                                                                                                                         aaaGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGggt                                                                                                                                                                                                                    <  1:74263/126‑1 (MQ=255)
                                                                                                                                                           aGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACt                                                                                     >  1:265496/1‑251 (MQ=255)
                                                                                                                                                                  ggCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGt                                                                              <  1:543804/251‑1 (MQ=255)
                                                                                                                                                                      gTATCGCTGGTCTGGAAGTTAAACGTATCATCAACGAACCGACAGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTga                                                                           >  1:485689/1‑250 (MQ=255)
                                                                                                                                                                               gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCTCCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTCCCTGGTTGAAGAATTCaag                                                                 <  2:18440/251‑1 (MQ=255)
                                                                                                                                                                                   ggAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTTCTGTCTGGACAAAGGTTCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGACATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGTCTTCGACAGCCCTCTTATCAACTACCTGGTTGAACAATTCaagaa                                                               <  1:378525/249‑1 (MQ=255)
                                                                                                                                                                                    gAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTTAAGACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCaagaaaga                                                            <  1:47119/251‑1 (MQ=255)
                                                                                                                                                                                      aGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATc                                                          <  1:198382/251‑1 (MQ=255)
                                                                                                                                                                                       gTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATc                                                          >  1:625801/1‑250 (MQ=255)
                                                                                                                                                                                         aaaaCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTGCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAgg                                                       <  2:55209/251‑1 (MQ=255)
                                                                                                                                                                                                  catcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGAc                                               <  2:390946/250‑1 (MQ=255)
                                                                                                                                                                                                     caACGAACCGACCGCAGCTGCGCTGGGTTACGGTCTTGACAAAGGTACTGGCAACCGTACTATCGCGGTTGATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTGCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTg                                            <  1:84621/250‑1 (MQ=255)
                                                                                                                                                                                                      aaCGACCCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCCACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATAGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAACCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTgcg                                          >  2:206573/1‑251 (MQ=255)
                                                                                                                                                                                                         gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGAACTGGCACCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGCGAAGACTTCGACAGCCGTCCGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCa                                        >  2:443930/1‑250 (MQ=255)
                                                                                                                                                                                                              gaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCCACGGTGATa                                                                                                                               >  1:4162/1‑158 (MQ=255)
                                                                                                                                                                                                              gaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACGATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGAGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATa                                                                                                                               <  2:4162/158‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCt                              <  1:612287/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGGATTGACCTGCGCAACGATCCGCt                              <  2:102943/251‑1 (MQ=255)
                                                                                                                                                                                                                               cTTACGGTCTCGACAAAGGTACTGGCAACCGTACTATCGCGGTTTGTGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGGTACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGCTGAAGAATTCAAGAAAGCTCAGGGGATTGACCTGCGCAACGCTCCGCTGGCAATGCAgcgc                 >  2:274912/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGCCGAAGTTGACGGCGAACAAACCTTCGAAGTTCTGGCCACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGCATTCAAGAAAGATCAGGGTATTGACCTGAGCAACGATCCGCTGGCAATGCAGCGCCTgca            >  2:543454/1‑247 (MQ=255)
                                                                                                                                                                                                                                          gACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGc      <  2:567615/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGTATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCa  <  2:557460/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCA  >  NC_000913/12472‑12960

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