Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,238,232 C→T P200L (CCT→CTT)  yjbH → putative lipoprotein YjbH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,238,2320CT95.9% 136.0 / ‑2.0 50P200L (CCT→CTT) yjbHputative lipoprotein YjbH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (26/21);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCG  >  NC_000913/4237997‑4238473
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
cAAGTGGCGGTTGGGGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTgg                                                                                                                                                                                                                                     <  2:59778/251‑1 (MQ=255)
     ggCGGTTGGTGCGCGGGATAGCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCGTGCCTCACTGTTTGgcggcg                                                                                                                                                                                                                                <  2:137676/251‑1 (MQ=255)
         gTTGGTGCGCGGGATATCGGCGGTACGGGGCTTTTTGATGCGGAATATCTATTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTGGGCCTGGGCTGGGGGTATTTGGGCACCAGCGTTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGa                                                                                                                                                                                                                            <  1:21339/251‑1 (MQ=255)
                              ggTACGGGGCTGTTTGATGGGGAATATGTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAAcc                                                                                                                                                                                                       <  2:342222/251‑1 (MQ=255)
                                 aCGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCt                                                                                                                                                                                                    <  1:22808/251‑1 (MQ=255)
                                  cGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTg                                                                                                                                                                                                   <  2:248116/251‑1 (MQ=255)
                                   ggggCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATATATCGACGGGAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTgc                                                                                                                                                                                                  >  2:86939/1‑251 (MQ=255)
                                             gATGCGGCATAACTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTAGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGCAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTTTCGACGTGAGCCAGATGTTCCACGGTCTTGCCTCCCTGTTTGGGGGCGTGGAATAACAGAAGCCCTGGCAACCGCTGCGCCTGACACt                                                                                                                                                                                        >  2:69593/1‑251 (MQ=255)
                                                   gAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGATCGAGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGACAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATGGACGGAAGCCGGATGTTCCACGGTATTGCCTCACTGTTTGGCGGCGTGGAATCCCAGACGCCCTGGCAACCGCTGTGCCTTAAACTGGAGt                                                                                                                                                                                   >  1:316635/1‑250 (MQ=255)
                                                    aaTATCTAGTTGCCTGCAAAGACAGGGGGCCGTTCGCTTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCCGCCAGTGGTCAATATAGCTATCGCGTTACCAGCTACAAACAGGCGGGATCTCACGCCGGGAGACAGCGAGTCCACGGTCTTGACTCACTGGTTGGCAGGGTTGAATACAGGACGCCCCGGCAAACGCTTGGCCTGAAACtcggggat                                                                                                                                                                                 >  2:173154/1‑244 (MQ=255)
                                                            gttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCa                                                                                                                                                                          >  1:42819/1‑250 (MQ=255)
                                                            gttgCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTAGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAATCAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGCCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCa                                                                                                                                                                          <  2:157432/250‑1 (MQ=255)
                                                                ccAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCaataat                                                                                                                                                                     <  1:253322/251‑1 (MQ=255)
                                                                   gCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTAt                                                                                                                                                                  <  1:97245/251‑1 (MQ=255)
                                                                    cAAAGCCTGGGGGCCGTTAGATTTTACGATCGGCTTGGGCGTGGGGTACTTGGGCGCCATCGGTAATGTGGAAAATCCACTCTGTTGCGACAGTGATAAATATTGCTATCGCGATAACAGCTACAAACCGGCGGGATCTATCGACGGGAGCCCGAAATTCCAAGGTCTTGCCCCACTGGTTGGCGGCGTGGACAACCCGACGCCCAGACAACCGCTACGCCTGGACCTGGAGtgtgcaggaaaaaatt                                                                                                                                                                    >  2:341030/1‑233 (MQ=255)
                                                                     aaaGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAACCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCCACCGCTGCGCCTGACACTGGAGTATGAAGGCAATAATTATc                                                                                                                                                                 >  2:219046/1‑250 (MQ=255)
                                                                         caggggggCCGTGCGATTGTACGCGCGGCCTGGGCTGGGGGTAGTTGGGCACCCGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTAGcagcag                                                                                                                                                            <  2:77772/248‑1 (MQ=255)
                                                                            gggggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGATCTGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGCTCTATCGACGGTAACCAGATGTTACACGGTCTTGCCTCCCTGTTTGTCGGCGTGGAAAACCAGACGCCATGGCAACCGCGGCGCCCGAAACTGGAGTATGAAGGCAATAATTATcaaccgggt                                                                                                                                                         >  2:206659/1‑244 (MQ=255)
                                                                             ggggCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGAtt                                                                                                                                                        >  1:360992/1‑251 (MQ=255)
                                                                                   gTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCg                                                                                                                                                  >  1:95679/1‑251 (MQ=255)
                                                                                       gATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGa                                                                                                                                               >  1:437177/1‑250 (MQ=255)
                                                                                                cccGGCCTGGGCTGGGGGAAGTTGGGCACCAGCGGTAATGGGAAAAACCCGCTCTGTTCAGCCTGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCTTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCa                                                                                                                                     <  2:275602/249‑1 (MQ=255)
                                                                                                          gctgggGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGATGGAGCAAAAaagtaag                                                                                                                           <  1:219046/251‑1 (MQ=255)
                                                                                                            tgggGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGtt                                                                                                                         >  1:409469/1‑251 (MQ=255)
                                                                                                             gggggTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGGGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGGGTATGAAGGCAATAAATATCAGCAGGATTTTGCCGGGAAGCTGGCGCAAAAAAGTAAGttt                                                                                                                        >  2:96135/1‑251 (MQ=255)
                                                                                                              ggggTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGTAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTa                                                                                                                       <  1:334619/251‑1 (MQ=255)
                                                                                                                ggTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAAc                                                                                                                     >  1:353434/1‑251 (MQ=255)
                                                                                                                  tATTTGGGCCCCAGCGGGAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATCAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGCTCTATAGTCGGTAGCCAGATGTTCCACGGTCTTGGCTCACTGTTTGGCGGCGTGGAATAACAGACGCCCAGGCAACCGCTGCACCTGAGACAGGGGTATGAATGCAAAAAATATCAGCAGAATTATGCCGGGAAGATGGAGCAAAAAAGTAAGTTTAACgc                                                                                                                   >  2:360430/1‑250 (MQ=255)
                                                                                                                            ccAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGAtt                                                                                                        >  1:290517/1‑251 (MQ=255)
                                                                                                                                 ggTAATGAGAAAAATCAGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTAACGCAGGTAGCCAGAAGTTCCACGGTCTTGCCTCACTGGTTGGCGGCGGTGACAACCAGACGCCCTGGCAACCGCTGCGCCTGAACATGGGGTATGAAGGCAATAATTATCCGCAGGCATTTGCCGGGAAGCTGGGGCAACAACGGAAGTTTACCGTCggagggtattatcg                                                                                                    >  2:398918/1‑239 (MQ=255)
                                                                                                                                     aggtgAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGtt                                                                                               <  2:457717/249‑1 (MQ=255)
                                                                                                                                       gtgAAAAATCCGCTATGTTCAGCCAGTGATAAATATTGCTAGCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGGAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGGCGGTGCGATTTATCGCGTTAc                                                                                             <  2:419777/251‑1 (MQ=255)
                                                                                                                                         gAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGGATGAAGGCAATAATTATCAGCAGGGTTTTGCCGGGAAGCAGGGGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCACGTTAccc                                                                                           >  2:288528/1‑250 (MQ=255)
                                                                                                                                                       ggtctgCCAGCGATAAATATTGCTAGCGCGATATCAGCTTCTAACAGGCGGGATCTTTCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACTCCCAGGCAACCGCTGCGCCCTAAACGGGAGTATGAAGGCAATAATTATCATCAGGATTTTGCCTGGAAGCTGTAGCAAAAAAGTAAGATTAGCGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGAc                                                                                <  2:223124/243‑1 (MQ=255)
                                                                                                                                                       gTTCAGCCAGTGATAAAAATTGCTATCGCAATAACAGCTACAAACAGTCGGGATCTAGCGACGGTAGCCAGATGTTCCATGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTGTCGCGGTACCTATTGTGCCGTCGtt                                                                             <  2:60305/251‑1 (MQ=255)
                                                                                                                                                               aGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCtt                                                                       >  1:493970/1‑249 (MQ=255)
                                                                                                                                                                 tGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGTCGGTAGCCAGATGTTCCACGGACTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTTCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAATAAGTTTAACGTCTGTGCGATTTATCGCGTTACCGATTGGGCCGGCGGTAACCTTAGcc                                                                   >  1:259902/1‑250 (MQ=255)
                                                                                                                                                                     aaaTATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATg                                                                <  1:96135/250‑1 (MQ=255)
                                                                                                                                                                               aTCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCa                                                       <  2:74365/249‑1 (MQ=255)
                                                                                                                                                                                   cgATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTTTGAACGTGGCAACACCt                                                 >  1:172045/1‑251 (MQ=255)
                                                                                                                                                                                           gCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGttt                                         >  1:165472/1‑251 (MQ=255)
                                                                                                                                                                                                    aGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGCCGTTCACCTTAGCTTTGAACGTGGCAACACCTTTATGTTTGGCGTTACg                                >  2:181264/1‑251 (MQ=255)
                                                                                                                                                                                                               aTCGACGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGGGTATGAAGGCAAGAATTATCAGCAGGTTTTGGCAGGGAAGCTGGAGCAAAAAAGTCAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGCCGTTAACATTAGCAATGAACGTGACAACACCTTTATGTTTTGCGTAACGTTgcgaccca                      >  2:262667/1‑245 (MQ=255)
                                                                                                                                                                                                                   ccGGTAGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACttt                 <  2:42819/250‑1 (MQ=255)
                                                                                                                                                                                                                    cgggaGCCAGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGAGGCAACACCTTTATGTTTGGCGTTACGTTGCGAACCCACttt                 >  2:301345/5‑250 (MQ=255)
                                                                                                                                                                                                                     ggTAGCCAGATGTTCCACGGTCTTGCCTAACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTaa               <  1:334151/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTgg                                                                                                                                         <  2:10037/123‑1 (MQ=255)
                                                                                                                                                                                                                            aGATGTTCCACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTgg                                                                                                                                         >  1:10037/1‑123 (MQ=255)
                                                                                                                                                                                                                                gtcccACGGTCTTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACTTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGcc    <  1:299570/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCACGGTCTTGCCGCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGGAAGTTTAACGTCGGCGCGTATTTTCGCGTTACCGATTGGGCCGACGTTAAACTTAGCTATGACCGTGGCAACACCTTTATGTGTGGCGGTACGGTGCGCACCAACTTTAAAGATcatctcccg  >  2:47412/1‑243 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
CAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGGATCTATCGACGGTAGCCAGATGTTCCACGGTCCTGCCTCACTGTTTGGCGGCGTGGAATACCAGACGCCCTGGCAACCGCTGCGCCTGAAACTGGAGTATGAAGGCAATAATTATCAGCAGGATTTTGCCGGGAAGCTGGAGCAAAAAAGTAAGTTTAACGTCGGTGCGATTTATCGCGTTACCGATTGGGCCGACGTTAACCTTAGCTATGAACGTGGCAACACCTTTATGTTTGGCGTTACGTTGCGCACCAACTTTAACGATCTGCGCCCG  >  NC_000913/4237997‑4238473

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: