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,7140CT91.5% 121.9 / ‑0.6 49G184G (GGC→GGT) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (21/22);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCA  >  NC_000913/12474‑12960
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCaccc                                                                                                                                                                                                                                              >  2:7442/1‑248 (MQ=255)
                                   ctgctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGggtggtg                                                                                                                                                                                                                <  1:138587/246‑1 (MQ=255)
                                      ctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTc                                                                                                                                                                                                        >  1:439649/1‑251 (MQ=255)
                                       tgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                       >  1:462373/1‑251 (MQ=255)
                                        gAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                       <  1:372025/250‑1 (MQ=255)
                                             gttGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACGTAATCTGTTAGCTCCGTCCCGGCATACTTTAACGAGGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTGCGATAttt                                                                                                                                                                                                 <  2:475883/249‑1 (MQ=255)
                                             gCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTTCTATCTCGGTTTATGTCCTGGTTGGTGGTACTTTCGATAttt                                                                                                                                                                                                 <  1:604183/251‑1 (MQ=255)
                                                 aaaaaaaTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAACCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGGGGTGGTACTTTCGATATTTCtat                                                                                                                                                                                             >  1:449273/1‑251 (MQ=255)
                                                     aaaTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATc                                                                                                                                                                                         <  1:7442/251‑1 (MQ=255)
                                                              aaaCCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATcgac                                                                                                                                                                                >  1:111663/1‑251 (MQ=255)
                                                               aaCCGCTGAAGATTACCCGGGAGAACCGTTAACTGAAGCTGTTACAACCGTACCGGCATACTTTAACGATGATCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAATGAACCGAACGCAGCTGAGCAGGCTTACGGACTGAACAAAGGTACTGGAAACCGTACTATCGCGGTTTATGACCTGGGTGGAGGTACTTTAGATCTTTCTATTAACGaaactgacg                                                                                                                                                                               >  2:27867/1‑245 (MQ=255)
                                                                 ccGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGAATCGCTGGTCTGGAAGTAAAACGTATCATAAACGAACCGACCGCAGCTGAGCTGGCTTACGGACTGGACAAAGAGACTGGCAACCGTACTATCGAGGTTGATGACATGGGAGGAGGTACTTTCGATATTTATATTATCGTAATcggcga                                                                                                                                                                              >  2:659627/1‑250 (MQ=255)
                                                                      gAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTg                                                                                                                                                                         >  2:35365/1‑250 (MQ=255)
                                                                      gAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGa                                                                                                                                                                        <  1:150810/251‑1 (MQ=255)
                                                                                   gTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTc                                                                                                                                                                                                        >  1:57514/1‑206 (MQ=255)
                                                                                   gTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAGACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTc                                                                                                                                                                                                        <  2:57514/206‑1 (MQ=255)
                                                                                       aCCGGTAACTGAAGCAGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCACTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGGGGTGGTACTTTCGCTATTTCTATTATCGAAATCGACGACGTTGAGGGCGAACAAAACTTCg                                                                                                                                                       >  2:4341/1‑251 (MQ=255)
                                                                                         cGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGa                                                                                                                                                      >  1:646815/1‑250 (MQ=255)
                                                                                                           tCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACAGTCTGGACAAAGGAACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGAACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCACCAAACggg                                                                                                                                   >  1:542179/1‑250 (MQ=255)
                                                                                                             aCCGTCCCGGCATACTTTAACCATGCTCAGCGTAAGGCAACCAAAGACGCAGGTCGTATCGCGGGTCTGGAAGTAAAACGTATCAACAACGAACCGACCGCAGCTGCGCTGGCGTACGGTATGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATAGTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAACACCTTCGAATTTCTGGCAACCAACGGTGa                                                                                                                                 <  2:164416/251‑1 (MQ=255)
                                                                                                               cggaCCGGCATTCTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTAGCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATa                                                                                                                               <  1:137966/248‑1 (MQ=255)
                                                                                                                 tACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAACTCGACGAAGTTGACGGCGAAAAGACCTTCGAAGTTCTGGCAACCAACGGTGATAcc                                                                                                                             >  1:330456/1‑251 (MQ=255)
                                                                                                                   ccGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCTACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCa                                                                                                                           <  1:440836/251‑1 (MQ=255)
                                                                                                                             ttAACGAAGCTCAGCGACAGGCCACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCACCGAAACGCCCGCAGCTGCACTGGCTTACGGACTGGACACAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGCCGGCGCAACAACCTTCGAAGTTCTGGCAACCAACGGAGAT‑CCCAccctgggtggg                                                                                                                >  2:49258/1‑242 (MQ=255)
                                                                                                                                           cGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTAACATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAg                                                                                                   <  2:646815/251‑1 (MQ=255)
                                                                                                                                              cAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCTACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCCCCTGGGTGGTGAAGCCTTCGACAGcc                                                                                                 <  2:88220/250‑1 (MQ=255)
                                                                                                                                              cAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGcc                                                                                                 <  1:359242/250‑1 (MQ=255)
                                                                                                                                               aGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGt                                                                                               >  1:217711/1‑251 (MQ=255)
                                                                                                                                                     ccAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATc                                                                                         >  2:56840/1‑251 (MQ=255)
                                                                                                                                                      cAAAGACGCCGGCCGTCTCGCTGGTCTGGCAGTTAAATGTATCCTCAACTAACCGATCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGATTCTGGCAACCGTACGATCGCTGTTTATGTCCTGGGTGGTGGTACTTTCCATATTTCTATTATCGAAATCGACGAAGTCTACGGCGAAAATACCGTCGAACTTCTGGCAACCAACGGTGATCCCCACCTTGGTGGTGAATACTTCGACAGCCTTCTGATCa                                                                                        <  2:542179/251‑1 (MQ=255)
                                                                                                                                                                       ccGCTGGTCGGGAATTAAAACGTAACATCAACGAACCGACCGCAGCTGCGCTGGTTTACGGTCTGGAAAAAGGTACTGGCAACCGGACTATCGCGGTTTATGCCCTGGGTGGTGGTACTTTCTATATTTCTATTATCTAAATCGACGAAGTTGACGGCGAAAAAACCTTCGACGTTCTGTCAACCAACGGTGATACCCACCTGGTTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTtaagaa                                                                       <  2:449273/250‑1 (MQ=255)
                                                                                                                                                                        cacgtgtctagaacgaaaacggttcatcaTCGAACCGACCGCAGCTGCGCTGTCTTACTGTCTGTAAACATGTACGTGCAACCTTACCATCGCGGTTTAGGACCTGGGTTGGGGTACATACGATATTTCTATTACCGAAATCGACGAAGTAGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGTTGACATCCCCCTGGGTGGTGAATACGTCGACAGTCGGCTTAGCAACTTCCTGGTTgaagaa                                                                       <  2:330456/227‑1 (MQ=255)
                                                                                                                                                                            ggTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCa                                                                   >  1:358995/1‑250 (MQ=255)
                                                                                                                                                                             gTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTCTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAACCCTTAGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGAGGTGAAGACTTCGACAACCGTCTGATCACCTACATGGTTGACGAAttaa                                                                   >  2:442594/1‑247 (MQ=255)
                                                                                                                                                                                tGGAAGTAAAACGTATCATCCACGAACCGACCGCAGCTGCGCTGGCATAAGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGCTATTTCTATTATCGAAATCGACGACGATGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGATGATACCAACCTGGGTGGTGAAAACTCCGACAGCCGGCTGATAAACAACCTGATTGAAGAATGCaagaa                                                               >  2:481796/1‑250 (MQ=255)
                                                                                                                                                                                             tATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATAGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGCAGTTCTGGAAACAAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATAAACTACCAGGTGGAAGGATTAAAGAAAGATCAAGGCATTg                                                 >  2:456975/1‑251 (MQ=255)
                                                                                                                                                                                                catcaACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGAcc                                              <  2:111663/251‑1 (MQ=255)
                                                                                                                                                                                                 agcaACGAACCGACCGCAGCTGCGGTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTAGCGCGGTTTATGACCTGGGTGGTGGTACTTTAGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCGGGTTGAATAATTCAAGAAAGATCAGGGCATTGACCt                                             <  2:322370/249‑1 (MQ=255)
                                                                                                                                                                                                    aaCGAACCGTCCTCAGCTGCGCTGGTTTACGGTCTGGCCAAAGGTACTTGCAACCGTACTATCGCGGTTTATGGCCTGGGTGGTGGTATTTACGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACGTCGACAGCCTTCTGAGCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTTCATTGGCCTgcg                                          <  2:652953/251‑1 (MQ=255)
                                                                                                                                                                                                      cGAACCGACCGCAGCTGCGCTGGCTTACGGGCTGGACAAAAGTACTGGCACCCGGACTCTCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGGCGGCGAAAAACCCTTCGAAGTTCTGGCACCCCACGGTGATACCCACCTGGGTGGGGACGACTTCGACCGCCATCTGATCCACTACCTGGTTGAAGTATTCaggaagtataagggacatgacctgcgc                                         >  2:174271/1‑228 (MQ=255)
                                                                                                                                                                                                       gAACCGACCGCAGCCGCGCGGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTCCGACAGCCGTCTGATCAACTACCTGGTTGAAGAAt                                                                      <  1:101373/220‑1 (MQ=255)
                                                                                                                                                                                                         accgaccgCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACg                                     >  1:170583/1‑251 (MQ=255)
                                                                                                                                                                                                                       cgcTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATg                       >  1:457894/1‑251 (MQ=255)
                                                                                                                                                                                                                        gcTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGc                      >  1:426652/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggcgtgcGGTCTGGACAATGGTACTGGCAAACGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTTACGGCTAAAAAACCTTCGAAGGTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCACCTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAgc                   <  2:678034/245‑1 (MQ=255)
                                                                                                                                                                                                                           ggCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAg                    <  1:56840/250‑1 (MQ=255)
                                                                                                                                                                                                                             cTTACGGTCTGGACAAAGGTACTGGCAAACGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTAACGACATCGACGAAGTTGGCGGCGAAAAAACCTTCGAAGTTCTGGCCACCAACGGTGATACACACCTGGGTGGTGAAGACTTAGACCGCCGTCTGATCAGATACCTGGTTGTAGAATGCAAGGAAGGTCAAGACATTGACCTGAGCAACGacccactgaccaagcaactc                 >  2:312103/1‑232 (MQ=255)
                                                                                                                                                                                                                                         aCAAAGGTACTGGCAACCGTAATATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGATTTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCg     >  1:101112/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaGGTACTGGCAACCGTACTATCGCGGTTTATGCCCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCa  <  1:28539/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
AGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCA  >  NC_000913/12474‑12960

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