Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,070,858 T→G K129N (AAA→AAC fbaA ← fructose‑bisphosphate aldolase class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,070,8580TG97.8% 137.8 / ‑3.8 45K129N (AAA→AACfbaAfructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (25/19);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

CAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGG  >  NC_000913/3070615‑3071099
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTc                                                                                                                                                                                                                                                <  1:53174/249‑1 (MQ=255)
    tCGGTTTATGCGTAATCAACTGTTTCAGGCTGGTGGTACAGGGCAGAAGCGCCCATGTGGCTGTTGTCCCCGCCGTATTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTAATGCCGATTTGGGACATGCGCTCCAGGTATTTAGATCAGATTGCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCTGTAGCTGCTAAGTGGTTTTCACCcg                                                                                                                                                                                                                                           <  2:259929/250‑1 (MQ=255)
           aTGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAAc                                                                                                                                                                                                                                   <  1:212800/251‑1 (MQ=255)
             gCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAg                                                                                                                                                                                                                                 <  1:187947/251‑1 (MQ=255)
                 aaTCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCg                                                                                                                                                                                                                             >  1:295585/1‑251 (MQ=255)
                       cGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATc                                                                                                                                                                                                                       >  1:259004/1‑251 (MQ=255)
                                 gCTGGGTGTACAGGGCAGAAGCGTCCATGTGGCGGTTGTCCACGCCGTCTGCGTCACCACAGGTGCAACCCAGGTCGATTTCCAGAGGCATGGCGAGGTGGGACATGCGCACCAGGTATTTAGAGCAGATTGCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGgcagc                                                                                                                                                                                                              <  2:316988/251‑1 (MQ=255)
                                            aGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTcgc                                                                                                                                                                                                   <  1:179279/251‑1 (MQ=255)
                                            aGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCGTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTcgc                                                                                                                                                                                                   <  1:127877/251‑1 (MQ=255)
                                                           aTGTGGCAGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCAACCGCGTCCAAAAGACCGTCGAGCCACGGCAGCAGTTTCTTAGCACAGTGGTccgcgtg                                                                                                                                                                                   >  2:137708/1‑245 (MQ=255)
                                                             gtgGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGc                                                                                                                                                                                  <  1:114713/250‑1 (MQ=255)
                                                                  tgttgtCCACGCCGTCTTCTTCAACACCGGTGCAACCAAGTACGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCACGTACATCAGACCGTCGATCCCCGGCAGCAGTTTCTTCGCACGGTGGTACATGTGCAGGAt                                                                                                                                                                             >  2:267550/1‑250 (MQ=255)
                                                                                         ccaccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGc                                                                                                                                                     >  1:308953/1‑251 (MQ=255)
                                                                                          caccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTGCAGcc                                                                                                                                                    <  1:353285/251‑1 (MQ=255)
                                                                                                     cccAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTAAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACAGGTAGCTGCGAAGTGGTTTTAACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGACGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGACATCTGGTGAAc                                                                                                                                         >  2:318796/1‑251 (MQ=255)
                                                                                                            agATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGAAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGATTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCACCAGACCGTCAATCCACGGCAGCAGTTTCTTCGCGCAGAGGTCAGTGTGCAGGATACCCGGAACCCAATACTGTTCAGCCATCTGGTGAACGTGATgc                                                                                                                                  >  2:242564/2‑251 (MQ=255)
                                                                                                                       gTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCag                                                                                                                           >  1:178605/1‑247 (MQ=255)
                                                                                                                           tGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATCTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATcgcg                                                                                                                   >  1:272004/1‑251 (MQ=255)
                                                                                                                                    tGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGAt                                                                                                           >  1:337260/1‑250 (MQ=255)
                                                                                                                                     ggACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATag                                                                                                         <  1:318796/251‑1 (MQ=255)
                                                                                                                                                  aGGTATTTAGAGCAGCTTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGACCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGCTAACCGGAACACCATAATGTTCAGCCATCTGGGGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACACTGCgg                                                                                            >  2:47357/1‑251 (MQ=255)
                                                                                                                                                     gatttagatacgatgttgaggTTCTCTTGCAGAGTGTCTTCAGACAGGTCGAGTATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCTAAGTGGTTTTCACCCGCGTCCAACAGCCCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCAGAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACTCTGCGGAAc                                                                                         <  2:162161/231‑1 (MQ=255)
                                                                                                                                                          agagCAGATTGCGACGTTCTCTTGCAGAGATTCTTCATAAAGGTCTATCATGTGAGAAGAGAACAGCGGCTTACCGGGAGCTGCGAAGTGGGTTTCCCCCGCGTCCAACATACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCATCACCCTGCGGAACGTCa                                                                                     <  2:306541/250‑1 (MQ=255)
                                                                                                                                                            agCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTc                                                                                      >  1:101397/1‑247 (MQ=255)
                                                                                                                                                             gCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGAtt                                                                                 >  1:88889/1‑251 (MQ=255)
                                                                                                                                                                gATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACATACCGTCGATCCACGGCAGCAGTTTCGTCGCGCTGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGGTCAGTCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCa                                                                              <  1:181768/251‑1 (MQ=255)
                                                                                                                                                                  tttCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGACGTGGTTTTCACCCGCGTCCAACAGACAGTCGATCCACGGCAGCAGTTTCTTCGCACAGTAGTCAGTGTGCAGGATAACCGGTACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCCGAGATCGCGCCCAGGCTAGCAGCACCCTACGGAACGTCAGATTTCa                                                                              >  1:57118/1‑249 (MQ=255)
                                                                                                                                                                  tttCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACg                                                                            >  1:140189/1‑251 (MQ=255)
                                                                                                                                                                  tttCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACAGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGACATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCCACCTGCGGAACGTCAGCTTTCACg                                                                            >  2:226086/1‑251 (MQ=255)
                                                                                                                                                                             cctttcatagaTGCTTCAGATCGGTCGATCATTTGGGAAGATAACAGCGTCTTACCGGTAGCTCTGAATTGGTTTTCACCCGCGTCCAACAGTCCGTCTATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCAGCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTCCTCGCCGTTACCAGc                                                                 <  2:195466/243‑1 (MQ=255)
                                                                                                                                                                                  cAGAGATTCTTAAGAAAGGTCGATCATATGAGAAGAGAACAGCGGCTTACCGATAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGACGATCCACTGCAGCAGTTTCTTCGCGCAGAGGTCAGTGTGCAGGATAACCGGAACACAATAATGTTCAGCCATCTGTTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCATCACCCACCGGAACGTCAGCTTTAACGCCTTTACCAgaccta                                                             >  2:18504/1‑245 (MQ=255)
                                                                                                                                                                                   agagATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATaa                                                            >  1:169662/1‑250 (MQ=255)
                                                                                                                                                                                       aTTCTTCAGACAGGGCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGaa                                                       <  1:383754/251‑1 (MQ=255)
                                                                                                                                                                                             cAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGcacc                                                  <  1:142824/250‑1 (MQ=255)
                                                                                                                                                                                                    gTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGCAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGAGCACCAGAAATCGCGCCCCGGCTAGCAGCACCCTGCGGAAAGTCAGATTTCACGCCTTTACCAGCGAGAAAGGAAGACCCACCGTTgg                                         >  2:202900/1‑251 (MQ=255)
                                                                                                                                                                                                      cGATCATTTGAGAAGAGAACTGCGGCTGTCCGGTAGCGTCGAATTTGTTTTCTCCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCATTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACACTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGgag                                       <  2:267297/251‑1 (MQ=255)
                                                                                                                                                                                                       gATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAAAGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTgg                                         >  2:314120/1‑248 (MQ=255)
                                                                                                                                                                                                         tCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGAAGTGGTCAGTGTGCCGGATAACCGGAACACAATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATCAAGGACGCACACACGTTGGAGAAc                                    >  2:46722/1‑251 (MQ=255)
                                                                                                                                                                                                          cATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACt                                   >  1:43445/1‑251 (MQ=255)
                                                                                                                                                                                                                  aagagaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGTCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATa                           <  1:153830/251‑1 (MQ=255)
                                                                                                                                                                                                                        aCAGCGGCTTACCGGTAGCTGCGAAGTAGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTTCCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGgc                     <  2:259004/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCtt                 >  1:220263/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTa               <  1:361929/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACtt           >  1:263716/1‑248 (MQ=255)
                                                                                                                                                                                                                                           tGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTCACTTTAGCAGCgg  >  1:2808/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGG  >  NC_000913/3070615‑3071099

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: