Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 908,889 C→T E136K (GAA→AAA)  ltaE ← low‑specificity L‑threonine aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,8890CT96.6% 198.6 / NA 59E136K (GAA→AAA) ltaElow‑specificity L‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/29);  minor base G (0/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAG  >  NC_000913/908648‑909128
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
cGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAg                                                                                                                                                                                                                                        <  1:579179/251‑1 (MQ=255)
    tAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCTACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTaa                                                                                                                                                                                                                                     <  2:303701/250‑1 (MQ=255)
      atggaaTCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCGTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGTGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAAttt                                                                                                                                                                                                                                  <  2:283584/247‑1 (MQ=255)
       tGTAATCACGATTACCGACGAGTAATGAAACGACTGACGTAGCAATACCTTTCGACAGTCAAATGGTGAACAAATCACAATATTGCGTGATTTATTTCAGTTAGCAGCCGGAAGCCACCACGGCATTAAAGATGAGCGCACCGTCAACATTCAGCTCCAGATTGCGATCGCGTGTAACTCCACATCCTTCTTTACGGTATTCCCGCGGACCCACTTTGCCGTTGGGGGTGTTTTTCATACTGAGTAATTTg                                                                                                                                                                                                                                 >  2:367333/1‑251 (MQ=255)
          aaTCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGt                                                                                                                                                                                                                               <  1:459907/250‑1 (MQ=255)
                           aGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATAtc                                                                                                                                                                                                             <  2:525525/251‑1 (MQ=255)
                           aGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCAGTTGTGGGTGTTTTTCAGACTGAGTAATTTGTTGCGGGCGAAATGGATAtc                                                                                                                                                                                                             >  1:555662/1‑251 (MQ=255)
                                 gAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCgg                                                                                                                                                                                                        <  1:418085/250‑1 (MQ=255)
                                          ggcgccccacgtccTTTCGAAAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCCCCGTCAACATGCAGCGCAAGATTGCGCGCGCTGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGAttt                                                                                                                                                                                               <  1:520040/238‑1 (MQ=255)
                                          ggCGTCCCTAGCCCTTTCGTCAGGCAAAAGGTGAACGAATCACAATATTGCGTGATCTCGTTCATTTCGCAGCCGTAAGCCACCGCGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATTCTTCTTTCAGGTATTCCCGCGTCAACACTGTGCCGTTCTGGGTGTTTTTCAGGTTGAGTAATTTGGTGCGGGCGAAATTGATATCGTCTGGTTTGAtgtg                                                                                                                                                                                              <  2:329402/251‑4 (MQ=255)
                                             gTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCAt                                                                                                                                                                                           <  1:268639/251‑1 (MQ=255)
                                                ccAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGc                                                                                                                                                                                        >  1:574802/1‑251 (MQ=255)
                                                 cAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGcc                                                                                                                                                                                       >  2:632616/1‑251 (MQ=255)
                                                  aaGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCa                                                                                                                                                                                      <  2:231085/251‑1 (MQ=255)
                                                  aaGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCa                                                                                                                                                                                      <  1:551540/251‑1 (MQ=255)
                                                                     aTGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGt                                                                                                                                                                   >  1:349570/1‑251 (MQ=255)
                                                                        gTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGcc                                                                                                                                                                >  2:502620/1‑251 (MQ=255)
                                                                        gTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTAATCGCCACTTTATCCAGCGGTAGCGTGcc                                                                                                                                                                >  2:147262/1‑251 (MQ=255)
                                                                         tGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGCTATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCg                                                                                                                                                               >  2:391658/1‑251 (MQ=255)
                                                                             cGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCgg                                                                                                                                                           >  2:486450/1‑251 (MQ=255)
                                                                                aTCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAg                                                                                                                                                        <  1:592110/251‑1 (MQ=255)
                                                                                     aatattgagttatagctTTCAGTTTGCTGCCGTAAGCCACCGCGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGGGGGTAAATTCCCATGCTTCTTTCAGGTATTCCGGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCTGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCATCcgcg                                                                                                                                                   <  2:512420/236‑1 (MQ=255)
                                                                                        atTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCt                                                                                                                                                <  1:491550/251‑1 (MQ=255)
                                                                                           gCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATg                                                                                                                                             <  1:570197/251‑1 (MQ=255)
                                                                                                    ctTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGtt                                                                                                                                         <  1:312031/246‑1 (MQ=255)
                                                                                                      ttCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGtt                                                                                                                                   <  1:290759/250‑1 (MQ=255)
                                                                                                       tCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGAATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGCCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGCGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTg                                                                                                                                  >  1:281377/1‑250 (MQ=255)
                                                                                                          gTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAAt                                                                                                                               >  2:204549/1‑250 (MQ=255)
                                                                                                                         ccaccaCGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCAccgcc                                                                                                               >  2:338912/1‑251 (MQ=255)
                                                                                                                                    ttAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGc                                                                                                    >  1:539505/1‑251 (MQ=255)
                                                                                                                                      aaaGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCtt                                                                                                  <  2:136196/251‑1 (MQ=255)
                                                                                                                                       aaGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTc                                                                                                 >  1:569014/1‑251 (MQ=255)
                                                                                                                                           tGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAAc                                                                                             >  1:124309/1‑251 (MQ=255)
                                                                                                                                                    cGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCTCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGTGCAAAATGGATGTCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGTTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTa                                                                                    <  1:63271/251‑1 (MQ=255)
                                                                                                                                                       caccatgaaacgcCAGAATACGCTAGCGGGTAAATTCCCCTGCTTATTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTCTCAGACTGAGTAATTTGTCGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCAGTCGGCAGCCGCGTCTATGGGTTGCGGCTGAATACTGCCCAGCACCGGCGCGCCAACGGCTTCAAACAGATAGTTCTgc                                                                                 >  2:358786/11‑251 (MQ=255)
                                                                                                                                                           aTGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGATATGCGCgg                                                                             >  1:482449/1‑251 (MQ=255)
                                                                                                                                                                     aGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTATCGTGCCGTCGGCAGCCGCGCCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAACCAGATAGTTATGCGCGGCATTACCGAc                                                                   >  2:461189/1‑251 (MQ=255)
                                                                                                                                                                     aGATTGCGCTCGCGGATAAATTCCCATGCTGCTATCAGGTATCCCCGCGGCGACCGTTTGCTGTTGTGGGTGTTTTGCAGCCTGAGTAATTTGGTGCGGGCGAAATGGATATCGGCGGGTGTGTTTTCCATCGTCACTTTATCCGGCGGTAGTGTGCCGTCGGCAGCCGCGTCTATGGGGTGCGGTTGAATATTGCCGAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGAc                                                                   <  1:432339/251‑1 (MQ=255)
                                                                                                                                                                         tGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAat                                                                >  1:84798/1‑250 (MQ=255)
                                                                                                                                                                            gcTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATAc                                                            >  1:162025/1‑251 (MQ=255)
                                                                                                                                                                               cgcgGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATActct                                                         <  2:200328/251‑1 (MQ=255)
                                                                                                                                                                                       aaTTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATggg                                                                                                                                           <  2:635576/161‑1 (MQ=255)
                                                                                                                                                                                       aaTTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATAGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATggg                                                                                                                                           >  1:635576/1‑161 (MQ=255)
                                                                                                                                                                                        atgcccATGCTTCTTTCGGGTATGCCCGCGGCAAGACTGTGCCGTTGTGGGTGTTTTGCAGACTGAGTAATTTGGTGCGGGGGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTGATCGAGCGGTAGCTTGCCGGCGGCAGCGGCGTCTATGGGTTGCGGTTGAATACTTCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCGTCGCCGCGt                                                <  2:443568/248‑1 (MQ=255)
                                                                                                                                                                                        aTTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGt                                                <  1:229362/251‑1 (MQ=255)
                                                                                                                                                                                         ttCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTTAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCGTACCGACAATAACCTCTTCGCCGCGtt                                               >  2:506745/1‑251 (MQ=255)
                                                                                                                                                                                              aTGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAg                                          <  1:37076/251‑1 (MQ=255)
                                                                                                                                                                                                gCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTg                                        >  1:195737/1‑251 (MQ=255)
                                                                                                                                                                                                      ttCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTCATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAACGGTAGCGTGCCGTCGGCAGCCGCGTCTATTGGTTGCGGTTGACTACTGCCCACCACCGCCGCGCCCCCGGCTTCAAACAGATAGATATTCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCCGTGACTg                                    >  2:74171/1‑249 (MQ=255)
                                                                                                                                                                                                      ttCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACt                                     <  1:632616/248‑1 (MQ=255)
                                                                                                                                                                                                             aTTCCCGCGGCAACACTTTGACGTTGTGGGTGTTTTTCAGTCTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAACGGTAGCGTGCCGTCGGGAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGGGCCACCGGCTTCAAACAGATAGTTAAGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTGCGCatcgtctaagcaggtcc                           >  2:7501/1‑235 (MQ=255)
                                                                                                                                                                                                                     ggCAACACTTTGCCTTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGtt                   <  1:51661/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggCAACACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGGGACTGAGCAGAGCGACAANGtt                   >  2:592722/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaCACTTTGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGc                >  2:366738/1‑251 (MQ=255)
                                                                                                                                                                                                                             tttGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCCGTGACTGAGCag                               >  2:463741/1‑231 (MQ=255)
                                                                                                                                                                                                                             tttGCCGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCag                               <  1:463741/231‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGcc       <  2:195737/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGTTGTGGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGcc       <  1:239610/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gtgGGTGTTTTTCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAg  <  2:569014/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
CGTTTAATGTAATCACGATTACCGACGAGTAATGAACCGACTGGCGTCCCAAGACCTTTCGACAGGCAAATGGTGAACGAATCACAATATTGCGTGATCTCTTTCAGTTCGCAGCCGTAAGCCACCACGGCATTAAAGATGCGCGCACCGTCAACATGCAGCGCCAGATTGCGCTCGCGGGTAAATTCCCATGCTTCTTTCAGGTATTCCCGCGGCAACACTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAG  >  NC_000913/908648‑909128

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