Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,571,895 A→G 20.2% L536P (CTC→CCC)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,571,8950AG20.2% 122.1 / 3.2 61L536P (CTC→CCC) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base A (14/33);  new base G (7/7);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCA  >  NC_000913/3571680‑3572110
                                                                                                                                                                                                                         |                                                                                                                                                                                                                          
ggttcaaatccgcccccAGCCGCGGGGCACCGGGGGGCTAGTTTTCGCCGCCTTCCAAGAGAGGCCAGTCGAGGCGGGCGTCAGGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttttt                                                                                                                                                                                                                <  1:132870/213‑1 (MQ=255)
 gttccggtccccccccggccGCTGGACCGCGTTGTGCCGATTAACGCCGCCGTCCCACGGAAGCCAGTTGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttttt                                                                                                                                                                                                                <  1:106389/208‑1 (MQ=255)
    cAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGGGACTCGCTTTTTTAc                                                                                                                                                                                                             >  2:516290/1‑227 (MQ=255)
     cgctCGCGCACCAGACGCTGGACGCCGTGGTGCCAGTTATCGCCGCCTGCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTGCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCg                                                                                                                                                                                                           <  1:113988/225‑1 (MQ=255)
     agtcccccccaaaggcGCTGGACCCCGGGTTTCCAGTTACCTCCGCCTTCCAACAGATGCCAGTCGAAGGTGGCGTCATGGTTCCACTCGCGGCCCGGGGCAAATTCGTTAACCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATTCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCg                                                                                                                                                                                                           <  1:573732/213‑1 (MQ=255)
      gATCGCGCACCAGACGCTGGACTCCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCg                                                                                                                                                                                                           <  2:231371/227‑1 (MQ=255)
        tCGGGCACCAGGCGCTGGTCACCGTGGTGCCAGTTAGCGCCGCCTTCCAACAGATGCCACTCGAGGCGGGCGTCATGTTTCCACTCGCGGCACTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCTAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTg                                                                                                                                                                                                         <  1:19363/227‑1 (MQ=255)
           cgcACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGCATTTCTGCCATGCGTCGCCCGGCATGCGGTCGGGAATCGAGTTTTTACCGTgggc                                                                                                                                                                                                      >  2:408149/1‑225 (MQ=255)
           cgcACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGGCGCCCGGCATGCGGGCGGGAATCGATTTTTTACCGTGGAc                                                                                                                                                                                                      >  2:372439/1‑227 (MQ=255)
                          ggacggtggggCCAGTTTTCGCCGCCTTTCAACAGATCCACGTCGAGGCTGGCGGCAAGGTTCCACCCGCGGCCCTGGGAAAATTCGTTGCCCATGAACAGGTGTTTCGTGCCCGGGAATGCCCACAACCAACCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGCCGCCCGGCATGCGGTCGAGAAGCGATTTTGTACCGTGGACCACTGCATCATGCGAc                                                                                                                                                                                      <  1:220418/223‑1 (MQ=255)
                              cGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTTTGCCATGCGTCTCCCGGCATGCGGTCGGGAATCGTTTTTTTACCGTGGACACCTTCACCCTCCGACTACg                                                                                                                                                                                  >  2:288860/1‑228 (MQ=255)
                                   tgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTCCCCCTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCCGGAATGCCCACATCCAGCCATCGTAGGCGCGCAGGTTCGCGCATGTCTGCCCTGCGTCGCCCGTCAAGCGGTGGGGAACCGATTTTTTACCCTGGGCCCCTtcaacctacgccaccggaag                                                                                                                                                                              >  1:253184/1‑210 (MQ=255)
                                        gTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttt                                                                                                                                                                                                                  >  2:591456/1‑186 (MQ=255)
                                             cgccccCTTCCAACAGAGGCCAGTCGAGGCCGGCGTCATGGTTCCCCTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGGTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGtttt                                                                                                                                                                   <  1:372439/228‑1 (MQ=255)
                                              ccccccctCCCACAGGTGCCAGTCGAGGCTGGGTGGCTGGTTCCCCTCGCGGGCGTGGGCAAATTCGTTACCCATGAAAAGGAGTTTCTTGCCCGGGGATGCCCACATCCAGCCTTAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGGGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGtttt                                                                                                                                                                   <  1:316700/220‑1 (MQ=255)
                                                 gccTTCCAACAGGTGCCATTCGGGGCGGGCGCCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAg                                                                                                                                                                <  1:241770/227‑1 (MQ=255)
                                                     tCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAgtg                                                                                                                                                              <  2:67166/225‑1 (MQ=255)
                                                            tggtccatttgAGGCTGGCGTCATTGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACAGCCAGGCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTaga                                                                                                                                                   <  1:286725/218‑1 (MQ=255)
                                                             aTGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTACCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGACGGCCCAAGTTTTCAgtgaaatattaaag                                                                                                                                                   >  1:246893/1‑217 (MQ=255)
                                                              tGCCAGTCGGGGCTGGGGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAgaga                                                                                                                                                 <  2:3006/228‑1 (MQ=255)
                                                               gCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGGGAATGGATATTTTACCGTGGACCCCTTCATCATGGGACACAGCCAGGCAGAAGTTTTCATGTAGGTTGTAAAGaa                                                                                                                                                >  1:629430/1‑228 (MQ=255)
                                                                  aGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCCTCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATcc                                                                                                                                             >  2:371435/1‑228 (MQ=255)
                                                                         ggttgCGTCATGGTTCCCCCCGGGGCCCGGGGAAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAgg                                                                                                                                       <  1:113470/223‑1 (MQ=255)
                                                                              cGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGCCGGGCATCGATTTTTTACCGTGGACCCCTTCCTCCTGCGGAACCGCCCGAGCGAAGTTTTCGGTGCAGTTGTGAGGAATCCCGTAGGTCaagt                                                                                                                                 >  1:388860/1‑225 (MQ=255)
                                                                                 tttGGTTCCCCCCGGGGCCCGGGGCAAATTCGTTACCCATGAACGGTAGGTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTAt                                                                                                                              <  1:589327/226‑1 (MQ=255)
                                                                                 tcgttttccctccccgGCCCTGGGCAAATTCGTTCCCCTTGAACAGTGCTTTCTTGCCCGGGAATTCCCCAATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTAt                                                                                                                              <  1:294009/214‑1 (MQ=255)
                                                                                            tCGGGGCCCTGGGCAAATTCGTTACCGAGGAACAGTAGTTGCTTGCCCGGGAATGCCCACATCCAGCCATAGTGAGCGAGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATAc                                                                                                                   <  2:584656/228‑1 (MQ=255)
                                                                                            ccggggCCCTGGGCAAATTCGTAACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCAAACCCCGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATAc                                                                                                                   <  1:653815/224‑1 (MQ=255)
                                                                                               cgGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGATTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGGGTAGTTGTGGCGAATCCCGAAGGTCAGTTTATCGTGATGATACTGa                                                                                                                >  2:63259/1‑228 (MQ=255)
                                                                                                 gCCCTGGGCAACTTCGTTACCCATGAACAGGGGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTGCGCGAATGTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAAGCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACg                                                                                                              <  1:371435/228‑1 (MQ=255)
                                                                                                            aTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCATTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGt                                                                                                   >  2:460283/1‑228 (MQ=255)
                                                                                                              ttgttacccttggcgcgtagtTTATTGCCCGGGAATGCCCTCCTCCGGCCATTGTAGGCGCGAGGGTTCGCGAATTTCTGCCATGCGTCGCCAGGGAGGCGGTCGGGAAGCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTc                                                                                                  <  1:115010/211‑1 (MQ=255)
                                                                                                              ttgtgttccttGAAGAGTAGTTTCTTGCCCGGGTACCGCCCATGCCAGCCCGAGGGGGCGCGCAGGGTTGCTAATTTCGGGCATGCGGCGCCCGGCATGCGGGCGGAAAACGATTTTTTACCGTGGACCACTTCAACATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTc                                                                                                  <  1:186438/217‑1 (MQ=255)
                                                                                                              tCGTTACCCATGAACAGTAGGTTCTGCACCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTc                                                                                                  <  1:408149/227‑1 (MQ=255)
                                                                                                                 ttACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGCCAACGGCCGGACGAAGGTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGGGATGATACTGCCGAAAAACCGGGTCGaa                                                                                               >  1:544053/1‑226 (MQ=255)
                                                                                                                   gacagtGAACGGTTGTTCTTTGCCCGGGAACGCCCACCGCCAGCATGGGTAGGGCGCGATGTTCGGGAATTTCTGCCCTGCGTTGCCCGGGATGGGGTCGGTAAACAATTTTTGACCGGGGACCCCGTCATCATGCGTCAACGTCAGGACGAAGTTGTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAACCCGGGTCGAGCt                                                                                             <  1:608139/222‑1 (MQ=255)
                                                                                                                   caccAGGAACAGTGGTTTCTTGCCCGGTAATGCCCACACCCAGCCAATGTAGGCGCGCAGGTGCGCGAATTTTCGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCt                                                                                             <  1:29711/225‑1 (MQ=255)
                                                                                                                         gAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGTGCGGGCTTCATGt                                                                                       >  2:547732/1‑227 (MQ=255)
                                                                                                                            tgtttttttctttcCCGGGCATGCCCCCCTCCAGCCATAGTCGGCGCCCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTc                                                                                  <  1:221483/222‑1 (MQ=255)
                                                                                                                             cggtgtTTTTTTCCCGGGAATGCCCACCTCCCGCCATAGTAGGGGGTCAGGTTCGGAAATTTCCGCCATGCGTCGCCCGGCATTCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCAACATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTcc                                                                                <  1:288860/224‑1 (MQ=255)
                                                                                                                              tttgtTTCTTGCCCGGGAGTGCCCACATCCAGCCATAGTGGGCGCGCAGGTTCGCGAATTTCTGCCGTGGGGCGCCCGGCATGCGGTCGGGAATCGATTTTTTACCGTGGACCACGTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCa                                                                               <  1:333685/225‑1 (MQ=255)
                                                                                                                              gtagttctgtGCCCGGGAATTCCCACATCCGGGCATAGTAGGCCCGCAGGTTCCCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCa                                                                               <  2:126937/219‑1 (MQ=255)
                                                                                                                                  tgttgtGCCCAGGAATGCCCGCAGCCAGCCATAGTAGGCGGGTAGGTTCGCGAATTTCGGCCAGTCGTCGCCCGGAAGGCGGGCGAGAAGCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTGCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTAGAGCTTCATGTAGTCCAGGgt                                                                           <  1:453859/223‑1 (MQ=255)
                                                                                                                                         cccGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTAcc                                                                                                                                                                                                            >  2:287869/1‑95 (MQ=255)
                                                                                                                                          ccGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCg                                                                                                                                                                                                           <  1:287869/95‑1 (MQ=255)
                                                                                                                                              gAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTTATGAGAATCCCCAAGGTCAGTTTTTCGTGATGATACTGGCGATAAACCCGGTCGAGCTTCATGCAGTCCAGGGGGTCCGGCATcc                                                                >  1:571536/1‑227 (MQ=255)
                                                                                                                                                        gctCCAACCATAGTAGGCGCGCAGGTTCGCGAATTTCCGCCATGCCGCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGtt                                                      <  2:388860/226‑1 (MQ=255)
                                                                                                                                                            ccgccattgtggggcctcAGGTTCGGGAATTTCCGCCATGCGTCCCCCGGCCTGCGGTCGGGAATCGATTTTTTACCGGGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCa                                                   <  1:516290/210‑1 (MQ=255)
                                                                                                                                                              gCCATTGTGGGGGCCGAGGTTCGCGAATTTCTGCCAGGCGTCGCCCGGCATGGGGTCGAGAATCGATTTTTTCCCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGGGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACtt                                                <  1:460283/228‑1 (MQ=255)
                                                                                                                                                                   agtagGCGCGCAGGTTCGCGAATTTCTTCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCATTGTAGTTGTAGAGAATCCCGAAGCTCAGCTTATCGTGATGATACTGCAGATAAACCGGGCCGAGCTTCCTTAGATCACAGGGTTCATGCATCCAACCAAGACTCCACTGGTAc                                            >  1:502562/1‑227 (MQ=255)
                                                                                                                                                                   agtagGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGCCAACGGCAGGACGAAGTTTTCAGTGTAGTTGTCGAGAATCCCGAAGGTCAGTTTTTCGTGATGATACTGACGATAAACCGGGCCGAGCTTCCTGGAGTCCAGGGTGGCATGCAACCAGACAAGGTTCCACtttaac                                            >  1:387101/1‑223 (MQ=255)
                                                                                                                                                                    tgagGCGCGCAGGTTGGCGAAGTTCTGCGCTGCGTTGGCCGGCATGCGGCGGAGAATCGATTTGTTCCCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTGGTAGAGAATCCCGAAGGTCAGGTTATCGTGATGATACTGACGAGAAACCGGGTCGAGCTTCATGTAGTCCATGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCa                                          <  1:372501/226‑1 (MQ=255)
                                                                                                                                                                      gaggggcgacggTTCGGGAATTTCCGCCCGGCGTCCCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAg                                         <  1:482686/218‑1 (MQ=255)
                                                                                                                                                                       ggCGCGCAGGCTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGCAGTTTTCAGTGCAGTTGTAGAGAATCCCGAAGGCAAGGTTATCGTGATGATACTGCCGCAAAAACGGGGCTGACTTCATGTCGTCCAGGGTGTCATGCTTCCAGCCGAGGTTCCACtttcaaccaa                                        >  1:508227/1‑219 (MQ=255)
                                                                                                                                                                         cgcgcAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGAAGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGTAGGTCCGTTTATCGGGAGGATACTGACGATAAACAGGGTCCAGCTTCATGTCGTCCAGGGGGACAGGAATCAAGACGAGGGTCCCCTTGAACCAGAAGc                                     >  1:207829/1‑228 (MQ=255)
                                                                                                                                                                                 gTTGGCGAATGTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGAc                               <  2:207829/227‑1 (MQ=255)
                                                                                                                                                                                                  aagcctCGCCCGGCATGCGGTCGGACATCGTTTTTTTTCCGTGGACCACTTCATCAGGCGACTACCGCAGGACGAAGTTTTCAGTTTAGTTGTAGAGAAGCCCGATGGTCAGTTTATCATGAGGATACTGACGATAAGCCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGa              <  2:207601/222‑1 (MQ=255)
                                                                                                                                                                                                  aTGCGTCGCCCGGCAGGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCAGGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGa              <  2:172145/227‑1 (MQ=255)
                                                                                                                                                                                                        cGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCATTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGaaa       <  2:94305/228‑1 (MQ=255)
                                                                                                                                                                                                        cGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGCTACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGGTGTTCCACTTGTACCAGAAGCCCAGACCGCCCCTATCCTGCGGCCgagaca       >  2:170576/1‑226 (MQ=255)
                                                                                                                                                                                                             ggCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCa  <  2:387101/228‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                          
GGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCA  >  NC_000913/3571680‑3572110

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