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,8580TG98.7% 224.9 / ‑4.3 74K129N (AAA→AACfbaAfructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (34/39);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGT  >  NC_000913/3070622‑3071100
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
gtgtATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTc                                                                                                                                                                                                                                      <  2:143959/251‑1 (MQ=255)
 tgtATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTcc                                                                                                                                                                                                                                     >  1:462298/1‑251 (MQ=255)
          aaTCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGAcc                                                                                                                                                                                                                             <  1:28453/250‑1 (MQ=255)
          aaTCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCg                                                                                                                                                                                                                            >  1:656276/1‑251 (MQ=255)
             cAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTc                                                                                                                                                                                                                          <  1:268232/250‑1 (MQ=255)
                  tCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGCCCGTCGATcc                                                                                                                                                                                                                     >  2:7829/1‑250 (MQ=255)
                      gcGGCTGGGTGTTCAGTGCAGAAGCGGCCATGTGGCTGTTGTCCATGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACAGGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTGAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTGACCGGTAGCTGCGAAGTGGTTTTCTCCCGCGTCCAACAGACCGTCGATCCACgg                                                                                                                                                                                                                 <  1:100211/249‑1 (MQ=255)
                       cGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGgc                                                                                                                                                                                                                >  1:329227/1‑250 (MQ=255)
                        ggCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGgc                                                                                                                                                                                                                >  1:669045/1‑249 (MQ=255)
                           tGGGTGTACAGTACAGAAAAGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGGAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCCGGTATTTAGCGCAGATTTCGCTGTTCTCTTGCAGTGATTCTTCAGACAGGTCGATCGTGTGAGAAGAGAACAGCTGCTTACCGGTAGCTGCGAATTGGTTTTCACCCGCTTCCAACAGCCCGTCGATCCCCGgcagc                                                                                                                                                                                                             >  2:408873/1‑249 (MQ=255)
                            gggTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTGCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGgcagca                                                                                                                                                                                                            <  1:363612/249‑1 (MQ=255)
                            gggTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTACACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                                                                                                                                                                                          >  1:497305/1‑251 (MQ=255)
                             ggTGTACAGTGCAGAAGTGTCCATGTGGCTGTTGTCCTCGCCGGCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCTCCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGtt                                                                                                                                                                                                         <  1:394815/251‑1 (MQ=255)
                              gtgtACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGGCGATCCACGGCAGCAGttt                                                                                                                                                                                                        >  2:263043/1‑251 (MQ=255)
                                  aCAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTttctt                                                                                                                                                                                                     <  1:7829/250‑1 (MQ=255)
                                    aGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCCTCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGGttcttagc                                                                                                                                                                                                  >  2:654037/1‑248 (MQ=255)
                                          gAAGCGCCCATGTGGCTGTGGTCCACGCCGTCTTCTTCCCCCCCGGTGCAACCCGGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCTAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTGGATCCACGGCAGCAGGTTCTTCGCGCAGTg                                                                                                                                                                                            <  2:263575/251‑1 (MQ=255)
                                           aaGCGTCCATGTGTCTGGTGACCACGCCGACTTCTTCACCACCGGTGCAACCCAGTGCGATGTCCATAGTCATGCCGAGTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGTTGTTCTCTTGCAGAGAGTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGATCGTCGATCCACGGCAGCATTGTCTGCGTGCAGTgg                                                                                                                                                                                           <  2:279310/251‑1 (MQ=255)
                                               gTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                                                                                                                                                                                         >  1:106954/1‑249 (MQ=255)
                                                 ccATGTGGCTGTTGTCCACGCCGTCTTCTTCACCCCCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTTGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAgtg                                                                                                                                                                                     <  2:177361/251‑1 (MQ=255)
                                                  cATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAgtgt                                                                                                                                                                                    <  1:181359/251‑1 (MQ=255)
                                                  cATGTGGCTGTTGTCCACGCCGTATTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTACGATGTTCTCTTGCCGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCAGCTTACCGGTAGCTGCGAAGGGGTTTTCACCCGCGTCCAACAGACAGTCGAGCCAAGGCAGCAGTTTCTTCGCGCAGTGGTCAgggt                                                                                                                                                                                    >  2:447658/1‑248 (MQ=255)
                                                   aTGTGGCTGTTGTCCACGCCTGCTTCTTCTCCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCATATTGCGATGGTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTTCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTGCGCGCAGTGGTCAGtgtg                                                                                                                                                                                   <  2:612999/251‑1 (MQ=255)
                                                     gtgGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGCAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGc                                                                                                                                                                                  >  1:279120/1‑250 (MQ=255)
                                                       ggCTGTTGTCCACGCCATCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGATAGCTGCGAGGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGAAGTTTCTTCGCGCAGAGGTCAGTGTGCAgg                                                                                                                                                                               >  2:374968/1‑251 (MQ=255)
                                                            ttgtCCACGCCGTCTTCTTCACCACCGGTGCAACCCCGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCAACGGCAGCAGTTTCTTCGCGCAGTGGTCAGGGTGCAGGATAAc                                                                                                                                                                          >  2:545089/1‑251 (MQ=255)
                                                                              tcACCACCGGTGCAACCCAGTTCGCTTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACTATAATg                                                                                                                                                           <  2:457851/248‑1 (MQ=255)
                                                                                  caccTGTGCAACCCAGGTCGATGTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATGTAGAGCAGATGTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGGCCAACAGACCGGCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGcc                                                                                                                                                    <  2:387838/251‑1 (MQ=255)
                                                                                   accGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCAccc                                                                                                                                                                                                                                           >  2:580745/1‑163 (MQ=255)
                                                                                   accGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCAATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCAccc                                                                                                                                                                                                                                           <  1:580745/163‑1 (MQ=255)
                                                                                    ccGGTTCAACCCAGTGCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTGTTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGGCCAACAGACCGGCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGGGTGCAGGATAACCGGAACACCATAATGTTCAGCCAt                                                                                                                                                  <  2:655039/251‑1 (MQ=255)
                                                                                      ggTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCt                                                                                                                                                >  1:7428/1‑251 (MQ=255)
                                                                                                      aTTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATgcgc                                                                                                                                >  1:89426/1‑251 (MQ=255)
                                                                                                         tCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAcc                                                                                                                             >  1:356732/1‑251 (MQ=255)
                                                                                                           ccgagTCATGCCGCTTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTAGATGTTCTCTTGCCGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCTTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAcc                                                                                                                             >  2:76586/3‑249 (MQ=255)
                                                                                                                           taggACATGCGCTCCAGGTATTTATAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGACGATCCACGGCAGCAGTTGCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGAACGCGCCCAGGa                                                                                                            <  1:226132/248‑1 (MQ=255)
                                                                                                                            tgtacttGCGCTTCAGGTAATTTGATCAGATTGCGATGTACTCTTGCAGCGATTCTTCAGATAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGTGAAGTGGTTTTCACCCGCGTCCAACATACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATATCCGGAACCCCATAAAGTACAGCCATCGTGTGAACGTGATGCGCACCAGAGATCGCGCCCAgg                                                                                                             <  2:404011/242‑1 (MQ=255)
                                                                                                                             ggACATGCGCTCCAGGCATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCAATCTGGTGAACGTGATGCGAACAAGAGATCGCGCCCAGGaccg                                                                                                         >  2:212354/1‑248 (MQ=255)
                                                                                                                               aCATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCg                                                                                                                                                                                                                            <  1:269092/134‑1 (MQ=255)
                                                                                                                               aCATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCg                                                                                                                                                                                                                            >  2:269092/1‑134 (MQ=255)
                                                                                                                                 aTGCGCTCCAGGTATTTAGAGCAGATTTCGATGGTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATagcagc                                                                                                     <  1:625452/251‑1 (MQ=255)
                                                                                                                                   gcgcTCCAGGTATTTAGAGCAGATTGCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAtcagca                                                                                                    <  1:473912/250‑1 (MQ=255)
                                                                                                                                   gcgcTCCAGGTATTTAGAGCAGATTGCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCAc                                                                                                   <  1:447658/251‑1 (MQ=255)
                                                                                                                                     gcTCCAGGAATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGCGCGTGCTTCGAACCAGAAATCGAGCCAAGGATATCAGCAccc                                                                                                 >  2:102572/1‑251 (MQ=255)
                                                                                                                                       tCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGAACGCGCCCAGGATAGAAGCAccatg                                                                                               >  2:440602/1‑248 (MQ=255)
                                                                                                                                          aGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGCTGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCgg                                                                                            >  2:545120/1‑251 (MQ=255)
                                                                                                                                             tATTTAGTGTAGATTTCGATGTGTTCTTTCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGACTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTTCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCTGAACACCATAATGTTCAGCCAGCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCCCCCTGCGGaa                                                                                          <  2:636368/250‑1 (MQ=255)
                                                                                                                                                  agagCAAATTTCGATGTTCTCTTCCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGCACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATACACGGCAGCAGGTTCTTCGCGCCGTGGTCAGTGTGCAGGATAACCGGCACACCCTAATGTTCAGCCATCTGGGGAACGTGATGAGCACCAGAGATCGCGCCCAGGATAGAAGCACCCTGCGGAACggc                                                                                      >  2:86770/1‑247 (MQ=255)
                                                                                                                                                        gATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCa                                                                              <  1:255398/251‑1 (MQ=255)
                                                                                                                                                             cGATGTTCTCGTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGcc                                                                          <  1:61218/250‑1 (MQ=255)
                                                                                                                                                               aTGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCttt                                                                       <  1:545120/251‑1 (MQ=255)
                                                                                                                                                                         gCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGAt                                                              <  1:220546/250‑1 (MQ=255)
                                                                                                                                                                          cAGAGATTATTCAGACAGGTCGATCATGTGAGAAGAGAAAAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTAAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGATGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATa                                                             >  1:428766/1‑250 (MQ=255)
                                                                                                                                                                           agagATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATaaa                                                           <  1:35929/251‑1 (MQ=255)
                                                                                                                                                                                 tcttcAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAg                                                      <  1:548073/250‑1 (MQ=255)
                                                                                                                                                                                     aagacagggcgttcagttgAGAAGAGAAATGCTGCTTGACGGTAGCTGAGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCTGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGATATCTCGACCAGGATAGCAGCCCCCTGCGGAACGTCAGATTTCACGACTTTACCAGCGATAAAGGAAGcacca                                                 <  2:354983/234‑1 (MQ=255)
                                                                                                                                                                                         cAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCATCCCCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGt                                             <  1:319321/251‑1 (MQ=255)
                                                                                                                                                                                         cAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGCTAAAGGAAGCACCACCGt                                             >  1:379071/1‑251 (MQ=255)
                                                                                                                                                                                            gTCGAACATGTGAGAAGAGCACAGCGGATTAACGGTAGCTGAGAAGAGGATTTCACCAGAGTACAACAGACCGACGATCAACGGCAGAAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATACTGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCGGGATAGCAGCACCATGCGGAACGTCGGATTTCACGCCTTTACCAGAGATAGAGGAAGAACCACCGtggg                                          >  2:43168/1‑248 (MQ=255)
                                                                                                                                                                                              cGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTgg                                          >  1:189399/1‑249 (MQ=255)
                                                                                                                                                                                              cGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTgg                                          >  1:398938/1‑249 (MQ=255)
                                                                                                                                                                                                agcATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGGGAACGTGATGCGCACCAGAGAACGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGaa                                      <  1:102572/249‑1 (MQ=255)
                                                                                                                                                                                                aTCATGTGTGAAGAGAACAGCGGCCTTCCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGGCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGCGAACGTGATGCGCACCAGAGATCTCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGaa                                      <  2:137494/251‑1 (MQ=255)
                                                                                                                                                                                                aTCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGGTAAACGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGCGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCATTTACCAGCGCTAAAGGAAGACCCACCATTGgaga                                       >  2:595705/1‑250 (MQ=255)
                                                                                                                                                                                                  cattggagaagagaaCAGCTGCGTACCGGTAGCTGCGAAGTTGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTGCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCGTTACCAGCGATAAAGGAAGCACCACCGTTGGCGAAc                                     <  2:196410/245‑1 (MQ=255)
                                                                                                                                                                                                  cATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACt                                    <  1:44504/251‑1 (MQ=255)
                                                                                                                                                                                                          aagagaaCAGCGGCTTACCGGTAGCAGCGCAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGAGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCATACATCTGGTGAACGTGATGCGCACCAGAGCTCGAGCCCAGGATCGCAGCACCCTGCGGAACGTCGGATTTCACGCCTTAACAAGCGTTCAAGGATGCCCCACCGCTGGCGAACTGAACGAt                             >  2:449948/1‑250 (MQ=255)
                                                                                                                                                                                                             agaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACTTGATGCGCAcc                                                                                                                             <  2:443495/151‑1 (MQ=255)
                                                                                                                                                                                                             agaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAcc                                                                                                                             >  1:443495/1‑151 (MQ=255)
                                                                                                                                                                                                              gaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGCTAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCg                        >  2:162254/1‑251 (MQ=255)
                                                                                                                                                                                                               aaCAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGCACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCgg                       >  1:31420/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCTCCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACGGAACGATAACCGGCGCTTTa                <  2:669045/251‑1 (MQ=255)
                                                                                                                                                                                                                         tcccGGGAGCTGCGAGGTGGTTTTCACCGGCGTCCAACAGTCCGTCGATCCACAGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATCTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCCCCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAAc              <  2:3101/248‑1 (MQ=255)
                                                                                                                                                                                                                                    gCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGt  <  1:374968/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
GTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCGGT  >  NC_000913/3070622‑3071100

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