Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,279,276 C→G F121L (TTC→TTG)  kbaZ → putative tagatose‑1,6‑bisphosphate aldolase 1 chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,279,2760CG96.5% 190.1 / ‑6.0 59F121L (TTC→TTG) kbaZputative tagatose‑1,6‑bisphosphate aldolase 1 chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base G (33/22);  minor base T (1/1);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTACTGATTGAAGCAACCTCCAATCAGG‑TAGACCAGTTCGGCGGTTATACC‑GGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑TCACCACGCCGGATGCCGCCCGCGCCACGCTGGA  >  NC_000913/3279043‑3279518
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
taacttaTTGAAGCAACCGCCA‑TGAGGGTAGACCAGTTCGGCGGTTATACC‑GGAAGGACTCCCGCCGATTTTCGCGGCTTTGTGTGTCAGCTCGTCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTg                                                                                                                                                                                                                                       <  2:267220/244‑1 (MQ=255)
          aaacaaactcccaacagggaaaaCCAGTTCGGCG‑TTGTACCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCAt                                                                                                                                                                                                                             <  1:130647/229‑1 (MQ=255)
                    cAATCAGG‑TAGACCAGTTCGGCGGTTATACC‑GGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCTTTGGTTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCCATGGCCAATGCCGATGGTTTGATTAAAAGCTACGTTGCGGCAGGGTTTAAAAAAATCCCCCTTGATTGCAGCATGTCCTGTCAgg                                                                                                                                                                                                                  >  1:30556/1‑251 (MQ=255)
                        cAGG‑TAGACCAGTGCGGCGGTTATACC‑GAAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGa                                                                                                                                                                                                              <  1:419348/251‑1 (MQ=255)
                           g‑TAGACCAGTTCGGCGGTTATACC‑GGAATGAAGCCCGCCGCTTTTCACGGGTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGGTTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGGCCATGCCGATGATTTGATTCAAAGCTACGTTGCGGCGGGATTGAAAAAAATCCACCTTGAATTCAGACTGTCCTGTCCGGActgtcc                                                                                                                                                                                                           >  1:535182/1‑246 (MQ=255)
                                           gttaatccc‑GGAAGGACCCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtga                                                                                                                                                                                            <  1:168328/244‑1 (MQ=255)
                                              atttaccggattctcgcccgccGTTTTTCGCGCCTTTGTTTGTCAGCTCTCCGACTCGGTGAATTTCCCGCAGGATGCGTTTATTCTGGGTGGTGTCCTTCTGGTGCCAACCCGCTGGCAACACCGGCCGGCCGCCCACGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTTAAAAAAATTCTCCTTGATTTCAGCATGTCCTGTCAGGACGATCCGATTCCCTGAACttat                                                                                                                                                                                           <  1:456067/236‑4 (MQ=255)
                                                      gAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTg                                                                                                                                                                                  >  1:225669/1‑251 (MQ=255)
                                                          ggcgcccgccGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTg                                                                                                                                                                              <  1:69991/249‑1 (MQ=255)
                                                                cgccgATTTTCGCGGCTGTGTTTGCCAGCTCGCCGACTCGTTGAATTGCCCCCAGGATGCGTTGATGCTGGGTGGTGACCATCGGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTTATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCGTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAAcgcg                                                                                                                                                                        <  2:352703/251‑1 (MQ=255)
                                                                           gcgGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGGTTGAAAAAAATCCACCTTGTTTGCAACATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGCGATCGTGGCTGCACGCGCCGCCCggccg                                                                                                                                                             >  2:554933/1‑247 (MQ=255)
                                                                                    tttgtCAGCTCGCCGACTCGTTGAATTTACCGCAGGATGCGTTGATTCTAGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAAGCAACGGCCAATGCCGATGATTTGATTAAAAGCTACGGTGCGGCAGGGTTGAACAACATCCAACTTGGTTGCATCATGTCCTGGGAGGACGGTCCGATTCCCCTAACTGATGTCATCGTTGCTGCACCTGCCGCCCGGCTGGCGACAGTg                                                                                                                                                    >  2:222399/1‑251 (MQ=255)
                                                                                     ttgtaggCTCGACGACGCGTGGAATCACCCGGAGGATGCGTTGAACCTGGGTGGTGACCAGCTGGGGCAAAACCGCTGGCAAAACCTGCCGGACGCTCAGGCAATGGCCAATGCCGAGGATTTGATTAAAAGCGACGTTGCGGCAGGATTGAAAAAAATCCACCATGATTGGAGCATGTCCTGTCAGTCCGAACCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCTGTCTGGCGAAAGTg                                                                                                                                                    <  2:371665/244‑1 (MQ=255)
                                                                                     ttgtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCCCCTTTATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTTATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAgggg                                                                                                                                                   >  2:128946/1‑248 (MQ=255)
                                                                                       gtCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGCAAGTGGCg                                                                                                                                                 >  1:241223/1‑251 (MQ=255)
                                                                                          agcgcgccagcTCGTGGAATTTCCCGCAGGATGCGGTGTTGCGGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGg‑aa                                                                                                                                             <  2:286831/241‑1 (MQ=255)
                                                                                            cTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGCCTGGCGAAAGGGGCggg                                                                                                                                               >  2:362310/1‑247 (MQ=255)
                                                                                             ccgccgcCTCGTTGAATTCCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCTATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGg‑aagaa                                                                                                                                          <  2:681354/250‑1 (MQ=255)
                                                                                                 cGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGAAATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTggggggaacaaccc                                                                                                                                       >  2:289616/1‑239 (MQ=255)
                                                                                                           aaTTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCCCCGCCCGTCTGGCGAAAGTGGCGGGAAAAAACC‑GT‑CTTgtacc                                                                                                                           >  2:167639/1‑247 (MQ=255)
                                                                                                              ttCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACt                                                                                                                         <  1:661502/250‑1 (MQ=255)
                                                                                                                cccGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTg                                                                                                                      <  1:104778/251‑1 (MQ=255)
                                                                                                                     aggctGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCCGCATGGCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGCACGCGCCGCCCGTCTGGTGAAAGTGGCGG‑GAGGAAcccggctttgacacctttgcgg                                                                                                                   >  2:694274/5‑230 (MQ=255)
                                                                                                                         tGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGg‑aa                                                                                                                                             <  2:139556/220‑1 (MQ=255)
                                                                                                                         tGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGg‑aa                                                                                                                                             >  1:139556/1‑220 (MQ=255)
                                                                                                                           cGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATAGTGGCTAAACGCGCCGCCCGTCTGGCGAAAATGGCGG‑GAGAAACCTGG‑Ctttgaccccattgctggaccca                                                                                                            >  2:439242/1‑230 (MQ=255)
                                                                                                                               gATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGGACACTTTGGCGAAGCCGATCTgg                                                                                                       >  1:20446/1‑251 (MQ=255)
                                                                                                                                 ttCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGa                                                                                                      <  2:225669/250‑1 (MQ=255)
                                                                                                                                  tCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGGAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAAAGCGCCGCCCGTCTGGCGAAAGTGGCGG‑GAAAAACCTGT‑CTTGAACACTTTGGCGAAGCCGGTCTggggg                                                                                                    >  1:26032/1‑248 (MQ=255)
                                                                                                                                     gggTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTAt                                                                                                  <  1:765228/250‑1 (MQ=255)
                                                                                                                                         ggtgACCATCTGGGGCCAAACCGCTGGCAAAGCCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCa                                                                                              <  1:202526/250‑1 (MQ=255)
                                                                                                                                         ggtgACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCa                                                                                              <  1:171322/250‑1 (MQ=255)
                                                                                                                                            gACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTgg                                                                                          <  1:616632/251‑1 (MQ=255)
                                                                                                                                                 gcTGGGGCCAAACCGCTGGCAAAACCGGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTTATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTGGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGGATGTCATTGGTACCg                                                                                     <  2:229680/250‑1 (MQ=255)
                                                                                                                                                       gCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGGGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGGGAAGCCGATCTGGGGTATGTCATTGGTACCGAAgggc                                                                               >  1:338356/1‑248 (MQ=255)
                                                                                                                                                        ccAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGg‑aag                                                                                                                                            <  1:413052/190‑1 (MQ=255)
                                                                                                                                                            aCCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGG‑CTTGAACACTTTGGCGACGCCGATCTGGAGTATGTCCTTGGGACCGAAGTGCCGGt                                                                           >  1:413754/1‑250 (MQ=255)
                                                                                                                                                              cGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAACATGTCCTGTCAGGACGCTCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTgg                                                                                                                                                   >  1:446753/1‑178 (MQ=255)
                                                                                                                                                                            gccggccgCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCTCCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATg                                                          <  1:498235/251‑1 (MQ=255)
                                                                                                                                                                                gccgCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCCTGAAAc                                                      >  1:474903/1‑251 (MQ=255)
                                                                                                                                                                                    cTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑GAGAAAACTGT‑CTTGTACACTTTGGCGAAGGCCAACTGGGGGATGTCATTGGGACCGAAGTGCCGGGAccccgcggccggcatgcaaactt                                                   >  2:692109/1‑229 (MQ=255)
                                                                                                                                                                                       aGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTAACTGGCGGCGCGCATGGAACCTTAAgg                                               >  1:410591/1‑250 (MQ=255)
                                                                                                                                                                                         gCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGAGCTGGGGTATGTCATTGGGACCGAAGTGCCGGGACCCGGCGGCGCGCCTGGAACCTTAAGCg                                              >  1:134945/1‑250 (MQ=255)
                                                                                                                                                                                          cAATGGCCAATGCCGATTATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTCCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAg                                            <  2:158468/251‑1 (MQ=255)
                                                                                                                                                                                            aTGGCCACTGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGCTGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAAACTGT‑ATTGAAAACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGGAATGCCGGAACCTGGCGGCGCGCGTGAAACCTTAAGCGAGc                                           >  1:728878/1‑250 (MQ=255)
                                                                                                                                                                                            aTGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGGCATTGGTACCGAAGTGCCGGGAACTGGCGGCGCGCATGAAACCTTAAGCGAGCt                                          >  1:485629/1‑251 (MQ=255)
                                                                                                                                                                                             tGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGCTCTGGAGGATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTg                                         >  1:578686/1‑251 (MQ=255)
                                                                                                                                                                                                 cAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTAGCGAAGCCGATATGGAGAATGTCATTGGTACCGAAGTGCCGGTACCTGGGGGAGCGCATGAAACCTTAAGCGAGCTTGCgg                                     >  1:596603/1‑251 (MQ=255)
                                                                                                                                                                                                    tGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACCCTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑Tcc                                 >  1:411932/1‑250 (MQ=255)
                                                                                                                                                                                                        gatgatTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTTGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑Tcaccac                             >  1:286820/1‑251 (MQ=255)
                                                                                                                                                                                                         atgatTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑TCACCACg                            <  1:182816/251‑1 (MQ=255)
                                                                                                                                                                                                            atttgattAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCAACTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAAACTGTGTTTGAAC‑CATTGTCGACGCCGATCTGGGGTATGGAGTTGTTCCCACAACGCCGGGTCCTCGGGGCGCGCATGAAACCTTAAGCGACCTGGCggggaccccccac                          >  2:585412/1‑240 (MQ=255)
                                                                                                                                                                                                                     aaaGCTACGTTGCGGCAGGATTGAAAAAAATCCACATTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATAGTGGCTGAACGCGCAGCCCGTCTGGCGAAAGTGGCGG‑AAGAAAACTGT‑CTTGGACAATTTGGCGAAGACGATCTTGAGTATAGCAATAGGACCGAAGTGCCGGTATCTGGGGGCGGGCATGAAACCTTAAGCGCGCCGGCGG‑Gcaccacacacaaatccgcc                >  1:586644/1‑238 (MQ=255)
                                                                                                                                                                                                                       aGCTGCGTTGCGGCAGGATTGAAAAAATACCGCCTTGATTGCATCAGGGCCGGTCAGGACGATCCGATTCCCTTAACGGATGACATCGTGGCGTAACGCGCCGCCCGTCTGGCGAAAGTGGCTG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCGGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑TCACCACGCCGGATGccgcccg              <  2:26032/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTTGCGAAACCGATCTGGGGGAAGTAATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCCGGCGGGCCACC‑CACCGGGTGccgcccgc             >  2:448418/1‑251 (MQ=255)
                                                                                                                                                                                                                         cTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGGAACCTGT‑CTTTGACACTTTGGCGGAGCCGATCTGGGGGATGGCAATGGTACCCAAGTGCCGGTAACTGGCGGCGCGCATGAAAACCTAAGCGGGCTGGCgggacccccgccgggtgtcccccgcc             >  2:137527/1‑227 (MQ=255)
                                                                                                                                                                                                                           aCGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTTGCGAAGCCGATCTGGGGTATGTCATTTGTACCGCAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑TCACCCCCCCCGGTGCCGCCCCCGcc          >  2:256384/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGAGCCGGGACCTGGCGGCGCGCATGAAACCTTAAGCGCGCTGGCGG‑TCACCACCCCCGATGCCCCCCGCGCCACGCTgg   >  1:451754/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGGAAACTTTGGCGAAGCCGATCTGGGGTATGGCATTGGGACCGAAGTGCCGGCAACTGGGGGCGCGACTGAAAACATAAGAGAGCTGGCGG‑TAACCAACCCGGGTGCCGACCGGGCCCCGCTggg  >  2:274922/1‑249 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
TTACTGATTGAAGCAACCTCCAATCAGG‑TAGACCAGTTCGGCGGTTATACC‑GGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGG‑AAGAAACCTGT‑CTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGG‑TCACCACGCCGGATGCCGCCCGCGCCACGCTGGA  >  NC_000913/3279043‑3279518

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: ‑