Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009641,246,9910TG32.0% 95.2 / 13.0 51V52G (GTG→GGG) thiDphosphomethylpyrimidine kinase, 4‑amino‑5‑hydroxymethyl‑2‑methylpyrimidine and 4‑amino‑5‑hydroxymethyl‑2‑methylpyrimidine pyrophosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base G (14/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 3.52e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAGAGCTTTGATCCATGGAACAAATGATATTGCCCGGGCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGGAGTTGCTCATTCCAAGAAGCTA  >  NC_000964/1246747‑1247237
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ggCAGAGCTTTGATCCATGGAACAAATGATATTGCCCGGGCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCCC‑AAAAAAAGCTT‑GGCGGTTACGGGGGGCATc                                                                                                                                                                                                                                                    >  1:168880/1‑250 (MQ=255)
    gaacTTTGGTCCATGGAACAAATGATATTGCACGGGCAAAATCCATTTAGCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGTCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGAGGTATTCGGAATGTCTGCGATTACGGCGTTAACCGCTC‑AAAATACGCTT‑GGCGTTCTCGGGGTGCATCCGCTg                                                                                                                                                                                                                                               <  2:37079/248‑1 (MQ=255)
          tGATCTATGGAACAACTTAGATTGCCCGGGTACAATCCATATATCATAAATGGATTGGGTAATAAAATGAGTGTAGTGTAATGAGTAGATATGAAGCATTGACAAGCGCCGTATACGATTCAGGCGGTGGCGCAGGGATACAGGCTGATATTAAAACCGTTCAAGAGCTTGATGGAGTCGGAATGTCTGCGTTTACGGCGGTTACCGCTC‑AAAAGGCTCTT‑GGTGTTCGCGGTAGGCATCCGCTTACTGt                                                                                                                                                                                                                                          <  2:99290/250‑1 (MQ=255)
            aTCCATGGAACAAATGATATTGCCCGGGCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTaga                                                                                                                                                                                                                                       <  1:147647/251‑1 (MQ=255)
                             tatTGCCCGGGCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACACCTC‑CAAATAAGCTT‑GGCGTTCACCGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGGCAgggc                                                                                                                                                                                                                      >  2:32733/1‑247 (MQ=255)
                                      ggCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAAACGCTC‑CAAATACGCTT‑GGCGTTTACGGGGGGGATCCGCTGACTGTAGGGACG‑GTGAGGCAGCAAAATTACGCTg                                                                                                                                                                                                            >  1:54494/1‑251 (MQ=255)
                                           aaTCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCt                                                                                                                                                                                                       <  1:134594/251‑1 (MQ=255)
                                                  ttATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGGAACAGCCCCAAAAT‑CGCTT‑GTGGCTCTCGGGGGGCATCCGTTGCCTGTTGAGACG‑GTTAGCAAGCAAATTGatgccgtttctgttgatt                                                                                                                                                                                                >  2:31290/1‑234 (MQ=255)
                                                    aTCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTg                                                                                                                                                                                              <  2:161320/251‑1 (MQ=255)
                                                       aTAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTgaga                                                                                                                                                                                           <  1:240990/251‑1 (MQ=255)
                                                                  gggTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑CAAATACGCTT‑GGCGTTCCCGGGGGGGATCCGCTGGCTGGAGAGACG‑TTGAGACAGCAAAATGGCGCTTGTGCTGGGGGTTTGGGACCGGATGCAGt                                                                                                                                                                                >  1:2463/1‑251 (MQ=255)
                                                                   ggTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑TGCGGTCACGGGGGGGATCCGCTGGCTGTAGAGACG‑GTGGGACAGCAAAATGACGGTGTTGCTGGGGGTTTGGGACCGGATGCAGtt                                                                                                                                                                               >  1:211791/1‑251 (MQ=255)
                                                                         aaaaGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCCCTC‑CAAATACGCTT‑GGCGGTCACCGGGGGCATCCGCTGACTTTAGAGACG‑GTGGGACAGGAAATTGACGCTGTTGCTGAGGGTTTGAGACCGGGTGCAggttagacc                                                                                                                                                                         >  1:174612/1‑243 (MQ=255)
                                                                             ggAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGaa                                                                                                                                                                     <  1:286767/251‑1 (MQ=255)
                                                                               agtgtagtgtAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATg                                                                                                                                                                   <  1:281037/251‑1 (MQ=255)
                                                                                   tagtgtaATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAAAAGCCC‑AAAATAAACTTTGGCGTCCA‑GGGGGGCCAACGCTGGATGGTGAGACGGTTG‑GACAGACAATTGTCGCTGTTGCTGTGGGTTTGCGAACGGGTGCGGTTAAGACGGGaagtctct                                                                                                                                                               >  2:296001/1‑245 (MQ=255)
                                                                                                    ataAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCCCTC‑AAAAAACGATT‑GGCGTTTACCGGGGGCATCCCCTGACAGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGTGGGTTTGTGAACGAGTGCAGTTAAGACGGGAAAGCTCTGGAAAGCCGATATGaat                                                                                                                                              >  1:72329/1‑249 (MQ=255)
                                                                                                        aGCAGTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTgaag                                                                                                                                          <  2:183283/251‑1 (MQ=255)
                                                                                                                    cGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGa                                                                                                                               <  2:285041/250‑1 (MQ=255)
                                                                                                                    cGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGa                                                                                                                               <  2:200202/250‑1 (MQ=255)
                                                                                                                               aTTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAg                                                                                                                   <  1:342950/251‑1 (MQ=255)
                                                                                                                                  cAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTa                                                                                                                 <  1:64569/250‑1 (MQ=255)
                                                                                                                                      cggcggcgCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGAt                                                                                                            <  2:315862/251‑1 (MQ=255)
                                                                                                                                         cggcgCATGGATACAGGCTGATATTAAAGCCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGTCATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTTCGGATTCa                                                                                                         <  1:37429/251‑1 (MQ=255)
                                                                                                                                            cgcAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATc                                                                                                      <  2:174612/251‑1 (MQ=255)
                                                                                                                                              cAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCgt                                                                                                    <  1:15328/251‑1 (MQ=255)
                                                                                                                                              cAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCCC‑AAAATACGCCT‑GGCGGTCACCGGGTGGATCCGCTGGCTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGGTTTGAGACCGGGTGCCGTTAAGACGGGAATGCTCTGGGAAGCCGGCATGATTGAAGAAGTGGCAAGAAAAATTTACGAGTACGGATTCAATCgt                                                                                                    >  1:1994/1‑251 (MQ=255)
                                                                                                                                                gggATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGGGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTTAGACCGGGTGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGGTTCAATCgtgt                                                                                                  >  1:169352/1‑251 (MQ=255)
                                                                                                                                                     aCAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAAtt                                                                                              <  1:62715/250‑1 (MQ=255)
                                                                                                                                                        ggCTGATATTAAAACCGTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGGTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAAttgttg                                                                                          <  2:279993/251‑1 (MQ=255)
                                                                                                                                                          cTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAAACGCTC‑AAAATACGCTT‑GGCGTTCCCGGGGGGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGGTGCAGTTAAGACGGGGATGCTCTGGAATGCCGACATGATTGAAGAAGTGGGAAGAAAAATTGACGAGTACGGATTCAATCGTGTAAttgttggt                                                                                        >  1:2133/1‑249 (MQ=255)
                                                                                                                                                             atatTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCg                                                                                     <  1:99056/251‑1 (MQ=255)
                                                                                                                                                               atTAAAACCTTTCAAGAGCTTGAAGTATTCGGAATGTCTGCGATTACGGCGGTAAACGCTC‑AAAATACACTT‑GGCGGGCCCGGGGGGCATCCGCTGACGTTTGAGAAG‑GTGAGACAGCAAAGTGACGCTGTTGGTGTGGATATGAGACCGGATGAAGATAAAAAGGGGGGGGTCTGGAAGGCCGACAGGGTTGAAGAAGAGGAAAGCAAAATTGAAGAGAAAGGGTTCAATTGTGTGAATGTTGATACGGt                                                                                   >  2:37745/1‑251 (MQ=255)
                                                                                                                                                                  aaaaCCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCCC‑AAAATACGCCT‑TGCGTTCACGGGGGGCATCCCCTGGCTGTAGAGACG‑TTTGGACAGcc                                                                                                                                                                                                                       >  1:168356/1‑115 (MQ=255)
                                                                                                                                                                             agagCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAGCCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAgg                                                                     <  1:179778/251‑1 (MQ=255)
                                                                                                                                                                                gCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAgg                                                                     <  2:279925/248‑1 (MQ=255)
                                                                                                                                                                                gCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAAACGCTC‑AAAATACGCCT‑GGCGGTCTCGGGGGGCATCCCCTGACTGTAGAGACG‑TTGAGACAGCAAATTGAAGCTGTTGTTGAGGATTTTAGACCGGGTGCAGTTAAGACGGGGATGCTCTGGAATGACGACATGATTGTAGAAGGGGGAAGAAAAATTGACGAGTACGGATTTAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAgggg                                                                   >  1:279925/1‑249 (MQ=255)
                                                                                                                                                                                    gATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCgggg                                                               <  1:145152/250‑1 (MQ=255)
                                                                                                                                                                                       gTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGTGGGGCAt                                                            <  1:299089/250‑1 (MQ=255)
                                                                                                                                                                                              gAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCGGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTtgctg                                                    <  1:256745/251‑1 (MQ=255)
                                                                                                                                                                                                              ggcggTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGAGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTg                                    >  2:184118/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGcc                                  >  1:177143/1‑250 (MQ=255)
                                                                                                                                                                                                                     aCCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGGGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCTCCtt                             <  2:102509/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccGCTC‑AAAATGCGCTT‑GGCGTGCACGGGGTGCTTCCGCTGGCTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGGCCGGATGCAGTTAAGACGTGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGAGCCGGTTATGTTCGCAAAGGGCGGTGCGTCCTGGCTGCGTGATGAGTCAGTTTCCACCTTa                            <  2:275611/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGGCG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGGCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAAGCGTGTAATTGTTGTGCCGGTTATGATCGCAAAGGGCGGGGCAGCCTTGCTGCGTGATGAGTCatgtgccggcttaaa                          <  2:305387/251‑15 (MQ=255)
                                                                                                                                                                                                                         cTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTaa                           >  1:151778/1‑249 (MQ=255)
                                                                                                                                                                                                                           c‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAgg                        <  1:19273/250‑1 (MQ=255)
                                                                                                                                                                                                                               caTACGCTT‑GGCGTTCACGGGGGGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATTCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGGAGtt                    >  2:72961/2‑251 (MQ=255)
                                                                                                                                                                                                                                      tt‑GGCGTTCCCGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGGAGTTGCTCAtt             >  2:208453/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGTAGTTGCTCATTCCaa         <  2:192473/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggggTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGGAGTTGCTCATTCCAAGAAGCtc  >  2:169467/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GGCAGAGCTTTGATCCATGGAACAAATGATATTGCCCGGGCAAAATCCATTTATCATAAATGGATTGGGTAATAAAAGGAGTGTAGTGTAATGAGTATATATAAAGCATTGACAATCGCCGGATCAGATTCAGGCGGCGGCGCAGGGATACAGGCTGATATTAAAACCTTTCAAGAGCTTGATGTATTCGGAATGTCTGCGATTACGGCGGTAACCGCTC‑AAAATACGCTT‑GGCGTTCACGGGGTGCATCCGCTGACTGTAGAGACG‑TTGAGACAGCAAATTGACGCTGTTGCTGAGGATTTGAGACCGGATGCAGTTAAGACGGGAATGCTCTGGAATGCCGACATGATTGAAGAAGTGGCAAGAAAAATTGACGAGTACGGATTCAATCGTGTAATTGTTGATCCGGTTATGATCGCAAAGGGCGGGGCATCCTTGCTGCGTGATGAGTCAGTTGCCACCTTAAAGGAGTTGCTCATTCCAAGAAGCTA  >  NC_000964/1246747‑1247237

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