Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,244 C→T L148L (CTC→CTT)  thrB → homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,2440CT98.4% 198.2 / ‑4.5 65L148L (CTC→CTT) thrBhomoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (30/33);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGAT  >  NC_000913/2998‑3484
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
tCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTggtg                                                                                                                                                                                                                                               >  1:195949/1‑251 (MQ=255)
               tttttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGCCACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGAACCGTGTTTTCTTGGTGGTATTCAGTtgatg                                                                                                                                                                                                                                 >  2:358672/1‑250 (MQ=255)
                  ttgtccggagcttgggaaGCAAACTCCCGTGGCGATGACTCTTGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGACTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTTCGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                                                                                                                              <  2:119142/234‑1 (MQ=255)
                  ttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCAGTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                                                                                                                              <  1:246690/250‑1 (MQ=255)
                  ttGCCAGGAGCTGTGCAAGCAAATTCCAGTGGCGATGACTCGGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                                                                                                                              <  2:429023/250‑1 (MQ=255)
                    gCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaa                                                                                                                                                                                                                           <  1:217992/251‑1 (MQ=255)
                       aggagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACCCTGTGGCAAGCCGCTTACTGACAATCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACGACCTGGCCACGTGTTTTATTGGTGGTATGCAGTTGATGATCgaaga                                                                                                                                                                                                                         >  2:192633/12‑250 (MQ=255)
                       aggagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCCTTCATTACGACAACGTGGCCCCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaaga                                                                                                                                                                                                                         >  2:243706/12‑250 (MQ=255)
                               ggcaaGCAAATGCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCTGTGGTCGAGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAATGACGAATCTCCGGCAGCATTCATTACGACAACGCGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAACACGCcat                                                                                                                                                                                                                <  2:154563/248‑1 (MQ=255)
                                           ccAGTGGCGATGACTCTGGAAAAGAATATGCCAAGCGGTTCGGGCTTAGGCTCCAGCGCCTGTGCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATcagccagc                                                                                                                                                                                                     <  2:121006/250‑1 (MQ=255)
                                              gTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAgg                                                                                                                                                                                                 >  1:81131/1‑250 (MQ=255)
                                                  cgttGACTCTGGAAAAGAATATGCCCATCGGTTGGGGCATAGGCACCAGCACCTGTTCGACGGTCGCGGAGCTGACGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGGGCTGGATGGACGATTCTCCGGCAGCATTCATTACGACAACGAGGCACCGTGATTTCTTGGTGGTATGCAATCGCTGATCGAAGAAAACCACATCATCAGCCAGCAAGTggcg                                                                                                                                                                                             >  2:174705/4‑248 (MQ=255)
                                                      gactcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCATTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAggg                                                                                                                                                                                          <  1:143692/246‑1 (MQ=255)
                                                           tGGAAAAGAATATGGAAAGCGGTTAGGGATTAGGCTAAAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGACTGAACACTGTGGCAAGCCGCTTAATGAAACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgat                                                                                                                                                                                    >  1:358998/1‑251 (MQ=255)
                                                            ggAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGGGCCAGGGTTtgatg                                                                                                                                                                                   >  2:338302/1‑251 (MQ=255)
                                                               aaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGt                                                                                                                                                                                <  1:414371/251‑1 (MQ=255)
                                                               aaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGt                                                                                                                                                                                >  1:310392/1‑251 (MQ=255)
                                                                       tgccactCGGCTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGTTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGCGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAGGAAAACGACATCATCAGCCAGCCGGTGCCGGGGTTTGATGAGTGGCTGTgg                                                                                                                                                                        >  2:213645/7‑251 (MQ=255)
                                                                                  tCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGt                                                                                                                                                             >  2:216472/1‑250 (MQ=255)
                                                                                     ggCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATcc                                                                                                                                                         <  1:230096/251‑1 (MQ=255)
                                                                                      gCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCg                                                                                                                                                        >  1:228234/1‑251 (MQ=255)
                                                                                      gCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCg                                                                                                                                                        >  1:347416/1‑251 (MQ=255)
                                                                                       cTTAGGCTCCGGCGCATGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACGGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCGGGAAGGACGAATCTCCGGCAGCGTTCATTACGACAACGTGGCACCGTGTTTTCTTGGAGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCAAGCCAGTGCCAGGGTTTGATGAGTGGCTGTGGGTG‑TGAGCGTATCCgg                                                                                                                                                       >  2:129041/1‑251 (MQ=255)
                                                                                                 agcgCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGt                                                                                                                                             >  1:1796/4‑251 (MQ=255)
                                                                                                 agcgCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGATGTGGGTGCTG‑GCGTATCCGGGGGTTCAAGt                                                                                                                                             >  2:91435/4‑251 (MQ=255)
                                                                                                    gCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGtctc                                                                                                                                          >  1:118650/1‑251 (MQ=255)
                                                                                                     ccTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAATGTCTCg                                                                                                                                         >  1:266239/1‑251 (MQ=255)
                                                                                                        gtgcGGTGGTCGCGGAGCTGAGGGCGTTGAATGAACACAGTGGCAAGCCGCGTAATGACTCTTGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAAGCTCCGGCAGGATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGTCTGTGGGTGCTG‑GCGTTTCCGGGGATTAAAGTCTCGAc                                                                                                                                       <  2:358998/247‑1 (MQ=255)
                                                                                                                          tGATGGCGATGAATGAACACTGTGGCACGCCTCTTACTGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTAATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCCAGCACGTGCCAGGGATTGATGAGTGGCTGTAGGTGCTG‑GCGAATCCGGGGATTAAAGTCTCGACGGCAGACACCAGGGCTAtt                                                                                                                    >  2:78141/1‑251 (MQ=255)
                                                                                                                           gATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAtt                                                                                                                    <  1:255399/250‑1 (MQ=255)
                                                                                                                            aTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGGTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAtttt                                                                                                                  <  1:23863/251‑1 (MQ=255)
                                                                                                                             tGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTCtttt                                                                                                                  >  1:483215/1‑250 (MQ=255)
                                                                                                                                cGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCg                                                                                                              <  1:156832/251‑1 (MQ=255)
                                                                                                                                    gaatgaaCACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGgcg                                                                                                           <  1:456506/250‑1 (MQ=255)
                                                                                                                                          accagttggCAAGCCGGTGAATGACACTCGTTTGCTGGCTTTGAGGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATc                                                                                                    <  2:24703/244‑1 (MQ=255)
                                                                                                                                           cacTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATcg                                                                                                   <  1:406777/251‑1 (MQ=255)
                                                                                                                                                   cAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAgg                                                                                           <  1:326903/251‑1 (MQ=255)
                                                                                                                                                          cgtaaggACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGattgcat                                                                                    <  2:223455/245‑1 (MQ=255)
                                                                                                                                                            tAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGattgcat                                                                                    <  1:453675/249‑1 (MQ=255)
                                                                                                                                                            tAATGACACTCGTGTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGattgcatt                                                                                   <  2:348065/250‑1 (MQ=255)
                                                                                                                                                              aTGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTgcg                                                                                >  1:392669/1‑251 (MQ=255)
                                                                                                                                                                    cTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCCCggg                                                                          >  2:479116/1‑251 (MQ=255)
                                                                                                                                                                                      tGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttc                                                        >  1:304675/1‑251 (MQ=255)
                                                                                                                                                                                       gggCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCa                                                                                                                                                                                                    <  2:61604/111‑1 (MQ=255)
                                                                                                                                                                                       gggCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCa                                                                                                                                                                                                    >  1:61604/1‑111 (MQ=255)
                                                                                                                                                                                        ggCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTGACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttcat                                                      <  1:259270/251‑1 (MQ=255)
                                                                                                                                                                                          cGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCttcattc                                                    <  1:46101/251‑1 (MQ=255)
                                                                                                                                                                                                              ccGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCg                                    <  1:243706/247‑1 (MQ=255)
                                                                                                                                                                                                                 gcagcaTTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCt                             <  2:411007/251‑1 (MQ=255)
                                                                                                                                                                                                                    gcaTTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACgg                                                                           <  2:224927/202‑1 (MQ=255)
                                                                                                                                                                                                                    gcaTTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACgg                                                                           >  1:224927/1‑202 (MQ=255)
                                                                                                                                                                                                                     cattcattTCGACAACGTGGCACCGGGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGGGCCAGGGTTTGATGAGGTGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGGCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGc                         <  2:366256/251‑1 (MQ=255)
                                                                                                                                                                                                                          attACGACAACGTGGAACCGGGTTTTCTTGGTGGTATGAAGTTGATGATAGAAGAAAACGCCATCATAAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGATG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCAAGGGCTATTTTACCGGCGACGCAACGTCACCGGGTTTGCAATGCGAACGGGCGCAATCTGGCTGGCTTCCTTCACGCCTACTATTCCCCTCAGCCTGAGCTTGcc                    >  2:251027/1‑251 (MQ=255)
                                                                                                                                                                                                                            tACGACAACGTGGCACCGTGTTTTCTGGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCcg                   <  2:81131/250‑1 (MQ=255)
                                                                                                                                                                                                                                aCAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAg              <  1:136292/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCCCGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCt            >  2:387818/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaCGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGAGGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCGAGTGCGATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCt            <  1:207515/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCGAGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCCAGCAAGTGACAGGGTTTGATGAGTGACTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTa                                                                                                                 >  1:29014/1‑150 (MQ=255)
                                                                                                                                                                                                                                  aaCGAGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTa                                                                                                                 <  2:29014/150‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACg                                                                                                                                      >  2:346694/1‑124 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACg                                                                                                                                      <  1:346694/124‑1 (MQ=255)
                                                                                                                                                                                                                                          cccGTGTTTTCTTGGTGGTAGGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGTGTGCGG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAg    <  2:195949/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGAt  <  1:216472/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttttCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATcgccgc                                                                                               <  2:439196/154‑1 (MQ=255)
                                                                                                                                                                                                                                                ttttCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATcgccgc                                                                                               >  1:439196/1‑154 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
TCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTGCCGCGAAGCTGATGAAAGAT  >  NC_000913/2998‑3484

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