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→TTGkbaZputative 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

TTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGA  >  NC_000913/3279043‑3279518
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
taacttaTTGAAGCAACCGCCATGAGGGTAGACCAGTTCGGCGGTTATACCGGAAGGACTCCCGCCGATTTTCGCGGCTTTGTGTGTCAGCTCGTCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTg                                                                                                                                                                                                                                       <  2:267220/244‑1 (MQ=255)
          aaacaaactcccaacagggaaaaCCAGTTCGGCGTTGTACCCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCAt                                                                                                                                                                                                                             <  1:130647/229‑1 (MQ=255)
                    cAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCTTTGGTTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCCATGGCCAATGCCGATGGTTTGATTAAAAGCTACGTTGCGGCAGGGTTTAAAAAAATCCCCCTTGATTGCAGCATGTCCTGTCAgg                                                                                                                                                                                                                  >  1:30556/1‑251 (MQ=255)
                        cAGGTAGACCAGTGCGGCGGTTATACCGAAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGa                                                                                                                                                                                                              <  1:419348/251‑1 (MQ=255)
                           gTAGACCAGTTCGGCGGTTATACCGGAATGAAGCCCGCCGCTTTTCACGGGTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGGTTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGGCCATGCCGATGATTTGATTCAAAGCTACGTTGCGGCGGGATTGAAAAAAATCCACCTTGAATTCAGACTGTCCTGTCCGGActgtcc                                                                                                                                                                                                           >  1:535182/1‑246 (MQ=255)
                                           gttaatcccGGAAGGACCCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACtga                                                                                                                                                                                            <  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)
                                                                                          agcgcgccagcTCGTGGAATTTCCCGCAGGATGCGGTGTTGCGGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGgaa                                                                                                                                             <  2:286831/241‑1 (MQ=255)
                                                                                            cTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGCCTGGCGAAAGGGGCggg                                                                                                                                               >  2:362310/1‑247 (MQ=255)
                                                                                             ccgccgcCTCGTTGAATTCCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCTATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGgaagaa                                                                                                                                          <  2:681354/250‑1 (MQ=255)
                                                                                                 cGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGAAATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTggggggaacaaccc                                                                                                                                       >  2:289616/1‑239 (MQ=255)
                                                                                                           aaTTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCCCCGCCCGTCTGGCGAAAGTGGCGGGAAAAAACCGTCTTgtacc                                                                                                                           >  2:167639/1‑247 (MQ=255)
                                                                                                              ttCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACt                                                                                                                         <  1:661502/250‑1 (MQ=255)
                                                                                                                cccGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTg                                                                                                                      <  1:104778/251‑1 (MQ=255)
                                                                                                                     aggctGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCCGCATGGCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGCACGCGCCGCCCGTCTGGTGAAAGTGGCGGGAGGAAcccggctttgacacctttgcgg                                                                                                                   >  2:694274/5‑230 (MQ=255)
                                                                                                                         tGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGgaa                                                                                                                                             <  2:139556/220‑1 (MQ=255)
                                                                                                                         tGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGgaa                                                                                                                                             >  1:139556/1‑220 (MQ=255)
                                                                                                                           cGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATAGTGGCTAAACGCGCCGCCCGTCTGGCGAAAATGGCGGGAGAAACCTGGCtttgaccccattgctggaccca                                                                                                            >  2:439242/1‑230 (MQ=255)
                                                                                                                               gATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGGACACTTTGGCGAAGCCGATCTgg                                                                                                       >  1:20446/1‑251 (MQ=255)
                                                                                                                                 ttCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGa                                                                                                      <  2:225669/250‑1 (MQ=255)
                                                                                                                                  tCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGGAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAAAGCGCCGCCCGTCTGGCGAAAGTGGCGGGAAAAACCTGTCTTGAACACTTTGGCGAAGCCGGTCTggggg                                                                                                    >  1:26032/1‑248 (MQ=255)
                                                                                                                                     gggTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTAt                                                                                                  <  1:765228/250‑1 (MQ=255)
                                                                                                                                         ggtgACCATCTGGGGCCAAACCGCTGGCAAAGCCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCa                                                                                              <  1:202526/250‑1 (MQ=255)
                                                                                                                                         ggtgACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCa                                                                                              <  1:171322/250‑1 (MQ=255)
                                                                                                                                            gACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTgg                                                                                          <  1:616632/251‑1 (MQ=255)
                                                                                                                                                 gcTGGGGCCAAACCGCTGGCAAAACCGGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTTATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTGGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGGATGTCATTGGTACCg                                                                                     <  2:229680/250‑1 (MQ=255)
                                                                                                                                                       gCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGGGGCGGAAGAAACCTGTCTTGAACACTTTGGGGAAGCCGATCTGGGGTATGTCATTGGTACCGAAgggc                                                                               >  1:338356/1‑248 (MQ=255)
                                                                                                                                                        ccAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGgaag                                                                                                                                            <  1:413052/190‑1 (MQ=255)
                                                                                                                                                            aCCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGGCTTGAACACTTTGGCGACGCCGATCTGGAGTATGTCCTTGGGACCGAAGTGCCGGt                                                                           >  1:413754/1‑250 (MQ=255)
                                                                                                                                                              cGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAACATGTCCTGTCAGGACGCTCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTgg                                                                                                                                                   >  1:446753/1‑178 (MQ=255)
                                                                                                                                                                            gccggccgCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCTCCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATg                                                          <  1:498235/251‑1 (MQ=255)
                                                                                                                                                                                gccgCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCCTGAAAc                                                      >  1:474903/1‑251 (MQ=255)
                                                                                                                                                                                    cTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGGAGAAAACTGTCTTGTACACTTTGGCGAAGGCCAACTGGGGGATGTCATTGGGACCGAAGTGCCGGGAccccgcggccggcatgcaaactt                                                   >  2:692109/1‑229 (MQ=255)
                                                                                                                                                                                       aGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTAACTGGCGGCGCGCATGGAACCTTAAgg                                               >  1:410591/1‑250 (MQ=255)
                                                                                                                                                                                         gCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGAGCTGGGGTATGTCATTGGGACCGAAGTGCCGGGACCCGGCGGCGCGCCTGGAACCTTAAGCg                                              >  1:134945/1‑250 (MQ=255)
                                                                                                                                                                                          cAATGGCCAATGCCGATTATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTCCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAg                                            <  2:158468/251‑1 (MQ=255)
                                                                                                                                                                                            aTGGCCACTGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGCTGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAAACTGTATTGAAAACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGGAATGCCGGAACCTGGCGGCGCGCGTGAAACCTTAAGCGAGc                                           >  1:728878/1‑250 (MQ=255)
                                                                                                                                                                                            aTGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGGCATTGGTACCGAAGTGCCGGGAACTGGCGGCGCGCATGAAACCTTAAGCGAGCt                                          >  1:485629/1‑251 (MQ=255)
                                                                                                                                                                                             tGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGCTCTGGAGGATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTg                                         >  1:578686/1‑251 (MQ=255)
                                                                                                                                                                                                 cAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTAGCGAAGCCGATATGGAGAATGTCATTGGTACCGAAGTGCCGGTACCTGGGGGAGCGCATGAAACCTTAAGCGAGCTTGCgg                                     >  1:596603/1‑251 (MQ=255)
                                                                                                                                                                                                    tGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACCCTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTcc                                 >  1:411932/1‑250 (MQ=255)
                                                                                                                                                                                                        gatgatTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTTGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTcaccac                             >  1:286820/1‑251 (MQ=255)
                                                                                                                                                                                                         atgatTTGATTAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACg                            <  1:182816/251‑1 (MQ=255)
                                                                                                                                                                                                            atttgattAAAAGCTACGTTGCGGCAGGATTGAAAAAAATCCAACTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAAACTGTGTTTGAACCATTGTCGACGCCGATCTGGGGTATGGAGTTGTTCCCACAACGCCGGGTCCTCGGGGCGCGCATGAAACCTTAAGCGACCTGGCggggaccccccac                          >  2:585412/1‑240 (MQ=255)
                                                                                                                                                                                                                     aaaGCTACGTTGCGGCAGGATTGAAAAAAATCCACATTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATAGTGGCTGAACGCGCAGCCCGTCTGGCGAAAGTGGCGGAAGAAAACTGTCTTGGACAATTTGGCGAAGACGATCTTGAGTATAGCAATAGGACCGAAGTGCCGGTATCTGGGGGCGGGCATGAAACCTTAAGCGCGCCGGCGGGcaccacacacaaatccgcc                >  1:586644/1‑238 (MQ=255)
                                                                                                                                                                                                                       aGCTGCGTTGCGGCAGGATTGAAAAAATACCGCCTTGATTGCATCAGGGCCGGTCAGGACGATCCGATTCCCTTAACGGATGACATCGTGGCGTAACGCGCCGCCCGTCTGGCGAAAGTGGCTGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCGGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGccgcccg              <  2:26032/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTTGCGAAACCGATCTGGGGGAAGTAATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCCGGCGGGCCACCCACCGGGTGccgcccgc             >  2:448418/1‑251 (MQ=255)
                                                                                                                                                                                                                         cTACGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGGAACCTGTCTTTGACACTTTGGCGGAGCCGATCTGGGGGATGGCAATGGTACCCAAGTGCCGGTAACTGGCGGCGCGCATGAAAACCTAAGCGGGCTGGCgggacccccgccgggtgtcccccgcc             >  2:137527/1‑227 (MQ=255)
                                                                                                                                                                                                                           aCGTTGCGGCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTTGCGAAGCCGATCTGGGGTATGTCATTTGTACCGCAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCCCCCCCGGTGCCGCCCCCGcc          >  2:256384/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGAGCCGGGACCTGGCGGCGCGCATGAAACCTTAAGCGCGCTGGCGGTCACCACCCCCGATGCCCCCCGCGCCACGCTgg   >  1:451754/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cAGGATTGAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGGAAACTTTGGCGAAGCCGATCTGGGGTATGGCATTGGGACCGAAGTGCCGGCAACTGGGGGCGCGACTGAAAACATAAGAGAGCTGGCGGTAACCAACCCGGGTGCCGACCGGGCCCCGCTggg  >  2:274922/1‑249 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
TTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTCAGCTCGCCGACTCGTTGAATTTCCCGCAGGATGCGTTGATTCTGGGTGGTGACCATCTGGGGCCAAACCGCTGGCAAAACCTGCCGGCCGCTCAGGCAATGGCCAATGCCGATGATTTGATTAAAAGCTACGTTGCGGCAGGATTCAAAAAAATCCACCTTGATTGCAGCATGTCCTGTCAGGACGATCCGATTCCCTTAACTGATGACATCGTGGCTGAACGCGCCGCCCGTCTGGCGAAAGTGGCGGAAGAAACCTGTCTTGAACACTTTGGCGAAGCCGATCTGGAGTATGTCATTGGTACCGAAGTGCCGGTACCTGGCGGCGCGCATGAAACCTTAAGCGAGCTGGCGGTCACCACGCCGGATGCCGCCCGCGCCACGCTGGA  >  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: