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,8580TG91.7% 144.2 / 1.5 48K129N (AAA→AACfbaAfructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (24/20);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGT  >  NC_000913/3070615‑3071100
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
cAGTGCGGTGTATGCGTAATAAACGGCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGGCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCAc                                                                                                                                                                                                                                             <  2:334170/251‑1 (MQ=255)
          tATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAAAAGCGTCCATGTGACTGTTGTCCACGCGGTCTTCTTCACCACCGGTGCAACCCGGTTCGTTTTCCAGAGCCATGCCGTTTTTGGACATGCGCTCCAGGAATTTAGACCAGATTTCGATTTTCTTTTGCAGAGATTCTTCAGACGGGTGGATCTTGTGAGAAGAGAACCGCGGCTTACCGTTAGCTGCGGAGGGTTTTTCACCCGAGTCCa                                                                                                                                                                                                                                    >  2:145663/1‑250 (MQ=255)
            tGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCCGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGGGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCGCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGTACATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTTGTTTTAACCCGCGTCCAAc                                                                                                                                                                                                                                  >  2:123254/1‑250 (MQ=255)
              cGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAAc                                                                                                                                                                                                                                  >  1:492522/1‑248 (MQ=255)
                tAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGAcc                                                                                                                                                                                                                             >  1:1356/1‑251 (MQ=255)
                    cAACGTCTTCCGGCTGGGTGTACAGGGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCCGAGGCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGCGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCg                                                                                                                                                                                                                         >  2:184715/1‑251 (MQ=255)
                            tCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACg                                                                                                                                                                                                                  <  1:505838/250‑1 (MQ=255)
                            tCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACg                                                                                                                                                                                                                  <  1:184715/250‑1 (MQ=255)
                                gCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCTTCCAACAGACCGTCGATCCACGgcagc                                                                                                                                                                                                             <  1:449096/251‑1 (MQ=255)
                                         aCAGTGCAGAAGCGTCCATGTGTCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTttcttc                                                                                                                                                                                                    <  1:145663/251‑1 (MQ=255)
                                         aCAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTttctt                                                                                                                                                                                                     <  1:121678/250‑1 (MQ=255)
                                                   aGCGTCCATGTGGCTGTTGTCCACGCCGTCTGTTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGGGGTTTTCACCCGCGTCCAACAGATCGGCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGt                                                                                                                                                                                          <  2:304289/251‑1 (MQ=255)
                                                          aTGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTTGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGtgtg                                                                                                                                                                                   <  1:15919/251‑1 (MQ=255)
                                                                  gttgtCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCAAGGCCGCAGTTTCTTAGCGCATTGGTCAGTGAGCAGGATa                                                                                                                                                                            >  1:331643/1‑250 (MQ=255)
                                                                    tgtCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAAc                                                                                                                                                                          <  2:251632/250‑1 (MQ=255)
                                                                                    ttcACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGtt                                                                                                                                                         >  2:232542/1‑251 (MQ=255)
                                                                                         caccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGCTGTTCTCTTGCCGAGATTCTTCAGACACGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCCCGGCAGCAGTTTCTTCGCGCAATTGTCATTGTGCAGGCTAACCGGCACACCCTAATGTTCAGcc                                                                                                                                                    >  2:38232/1‑251 (MQ=255)
                                                                                                  aaCCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTg                                                                                                                                            >  1:28985/1‑250 (MQ=255)
                                                                                                  aaCCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTg                                                                                                                                            >  1:349205/1‑250 (MQ=255)
                                                                                                   aCCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGGCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAAGGTGCAGCCATCTGGTGaa                                                                                                                                          <  2:438019/251‑1 (MQ=255)
                                                                                                       aGTTCGATTTCCATAGTCATGCCGATTTTGGACAGGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGGCCGTCGATCCACGGCAGCAGTTGCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTg                                                                                                                                      <  2:121222/251‑1 (MQ=255)
                                                                                                            gATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATgcg                                                                                                                                 <  2:157945/251‑1 (MQ=255)
                                                                                                             aTTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATgcgc                                                                                                                                <  1:232542/251‑1 (MQ=255)
                                                                                                                tCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAcc                                                                                                                             >  1:282777/1‑251 (MQ=255)
                                                                                                                  cAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAcc                                                                                                                             >  1:254192/1‑249 (MQ=255)
                                                                                                                             cGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGccc                                                                                                                >  1:420911/1‑251 (MQ=255)
                                                                                                                              gATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGccc                                                                                                                >  1:495148/1‑250 (MQ=255)
                                                                                                                                    ggACATGCGCTCCAGGGATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATa                                                                                                          <  1:334286/250‑1 (MQ=255)
                                                                                                                                        aTGCGCTCCAGGTATTTAGAGCAGATTTCGATGTGCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTGCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATagcagc                                                                                                     <  1:24356/251‑1 (MQ=255)
                                                                                                                                             cgccAGGGATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCt                                                                                                <  2:28985/249‑1 (MQ=255)
                                                                                                                                              tCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCCGTGTGCAGGATAACCGGAACACCATAATGTTCAGCAATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCCGGATAGCAGCACCCTg                                                                                               >  2:483917/1‑251 (MQ=255)
                                                                                                                                                  ggTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCgg                                                                                            >  1:194250/1‑250 (MQ=255)
                                                                                                                                                      tttAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGCTCATGTGAGAAGCGAACAGCGTCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGCTCCACGGCAGCAGTTTCTTCGCGCCGTGGTCAGCGAGCAGGCTAACCGGAACACCATACTGTCCGACCATCTGAAGACCGTGATGCGCACCAGAGATCGCGCCCGGAATAGACGCACCCAGCGGAACgg                                                                                       >  2:243854/1‑250 (MQ=255)
                                                                                                                                                                    cGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAGCACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGACCGTCAGATTTCACGCCt                                                                         >  2:403379/1‑251 (MQ=255)
                                                                                                                                                                        gTTCTTTTGCAGAGATTGAACAGACAGGTCGATCATGTGAGAAGAGACCAGCGGCTTACCGAAAGCCGCGAAGTGGCTTTCACCCGCGTACAACAGACCGTCGATCAGCGGCAGCAGTTTCTTCGCGCAGTGGTCAGCGTGCAGGATAACAGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCCACAGAGATCGCGCCCAGGCTAGCAGCACCCTGCGGCACGTCAGTTTTCACGCCtttcc                                                                     >  2:453718/1‑249 (MQ=255)
                                                                                                                                                                             ctTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGCACGTCAGATTTCACGCCTTTACCAGCg                                                                >  1:182881/1‑251 (MQ=255)
                                                                                                                                                                               tGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGAt                                                              >  1:486456/1‑251 (MQ=255)
                                                                                                                                                                               tGCAGAGAGTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCGTCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCGCGCCTTTACCAGCg                                                                <  1:56204/249‑1 (MQ=255)
                                                                                                                                                                                             aGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCATCCATCTGGTGAACGTGATGCGCCCCAGAGATCGCGCCCAGGATCGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGACGcacccc                                                >  2:390148/1‑249 (MQ=255)
                                                                                                                                                                                                cAGGTCTATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACTGCAGCAGTTTCGTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGt                                             <  1:499809/251‑1 (MQ=255)
                                                                                                                                                                                                cAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCTGTAGCTGCGAAGTGGTTTTCTCCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGt                                             <  1:144726/251‑1 (MQ=255)
                                                                                                                                                                                                 aGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATa                                                             >  1:127125/1‑234 (MQ=255)
                                                                                                                                                                                                 aGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATa                                                             <  2:127125/234‑1 (MQ=255)
                                                                                                                                                                                                     ggATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCAGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGgag                                        <  1:433679/250‑1 (MQ=255)
                                                                                                                                                                                                          aTGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTg                                   >  2:317844/1‑251 (MQ=255)
                                                                                                                                                                                                          aTGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTg                                   >  1:354816/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTac               >  1:438110/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCTGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCACAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGt  <  2:182881/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
CAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGT  >  NC_000913/3070615‑3071100

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