Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC40.9% 209.7 / 20.9 116V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (52/16);  new base C (0/47);  total (52/64)
Fisher's exact test for biased strand distribution p-value = 8.70e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

TGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTC  >  NC_000913/3260906‑3261200
                                                                                                                                                     |                                                                                                                                                 
   aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                               >  1:1026875/1‑151 (MQ=255)
     tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgc                                                                                                                                             >  2:1904845/1‑151 (MQ=255)
      gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                            >  1:1013459/1‑151 (MQ=255)
       aaTGATATTCAGCGCGCTGATTTCCTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                            <  1:1434775/150‑1 (MQ=255)
       aaTGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgccg                                                                                                                                           <  1:898639/151‑1 (MQ=255)
             atTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAtt                                                                                                                                     >  2:906974/1‑151 (MQ=255)
              ttCGGCGCGCTGTTGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCATTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAAt                                                                                                                                      <  2:1038046/149‑1 (MQ=255)
              tccagcggtcgggtGTCCTGCCCCGCCAGCGAGTCCATGAAGTGTGCGAGGCTGTCCTTCCCTTTGTCGTAACCCGGCACGCCGTCCAACGGCGGTGCTGTTTGGGCCCCGCCCTAAATCTCAAGCTTTTCGTCCACCCCGCCGTTAAt                                                                                                                                      <  1:1352209/136‑1 (MQ=255)
                cAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                   >  2:141767/1‑150 (MQ=255)
                      gcTGATGTACTGCCCCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTCCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                            <  1:120618/151‑1 (MQ=255)
                      gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                            >  1:816566/1‑151 (MQ=255)
                      gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                            >  2:234128/1‑151 (MQ=255)
                         gATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcagc                                                                                                                         >  1:1370309/1‑151 (MQ=255)
                         gATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcagc                                                                                                                         >  2:1746014/1‑151 (MQ=255)
                         gATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAATTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcagc                                                                                                                         >  1:86978/1‑151 (MQ=255)
                          aTGTACTGCCCCCCCACCCAGTCAATGAAGTGAGCGAGGCTGTCCATCACTTGGTCGTAATCCAGCCCGTCGCCCACCAGCGGTGCTTTTTTGGGCCCGCCCGGAATCTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGagca                                                                                                                           <  1:1745749/148‑1 (MQ=255)
                          aTGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                        >  1:1797202/1‑151 (MQ=255)
                           tttACTGCCCCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCGCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGagcag                                                                                                                          <  2:1539263/146‑1 (MQ=255)
                           tGTACTGCCCCGCCAGCCAGCCCATGAAGTGATCGGGGCTGTCCATCCCTTTGTCGTAACCCAGCACGTCGTCCATCGGCGGTGTTGTTTTCGCCCCGACCTGAATCTTCAGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                        <  1:900863/150‑1 (MQ=255)
                            ttcttGCCCCCCCGGCCGGTCCAGAAAGTGATCGAGGCTGTCCACCACTTTTGCGAAATCAAGCCCGTCGCCCAACAGGGGTGCTGTTTTCGGCCCCACCTGATTCTTCTGCTTCTCGTCCGCCCCGCCGTTAATTGCGTAGAGCAGCGt                                                                                                                       <  2:658023/146‑1 (MQ=255)
                             tccggccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGagcag                                                                                                                          <  2:1502872/139‑1 (MQ=39)
                             tACTCCCCCGCCAGCCAGTCCAGGAAGTGATCGAGGCTGCCCATCATTTTGTCGTAATCCAGCACGCCGTCCATCAGCGGGGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                      <  2:908532/150‑1 (MQ=255)
                               cTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                   >  2:1294619/1‑151 (MQ=255)
                                tGCCCCCCCAGCCAGTCCATGAAGTGACCGAGGCTGTCCATCCCTTTGTCGTAATCCACCCCGTCGCCCATCAGCGGTGCTGTTTCCGCCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGagcagc                                                                                                                         <  2:414286/144‑1 (MQ=255)
                                  ccccGCCAGCCAGTCCAGGAAGTGATCGGGGCTGTCCATCTTTTTGTCGTAACCCAGCGCGTCGCCCACCAGGGGTGCTGTTTTGGGCCCGCCCAGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                                  <  1:604941/147‑1 (MQ=255)
                                  cACCGCCCGCCAGTCCATGAAGTGATCGAGGCTGTCCATCCCTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGTTGCTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                                                                                 <  1:1762964/150‑1 (MQ=255)
                                    cccccggccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTACTCCAGCCCGTCTTCCACCGGCGGTGCTGTTTTCGCCCCGACCTGAATCTTCACCTTCTCTTCCACCCCGCCGTTAATTGCGTAGAGCAGCGt                                                                                                                       <  2:460272/136‑1 (MQ=255)
                                    ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGCCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCa                                                                                                               <  2:481567/150‑1 (MQ=255)
                                    ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAACCCACCACGTCGCCCACCAGCGGTGCTGTTTTCGCCCCGACCTGAATTTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                   <  2:1410067/146‑1 (MQ=38)
                                    ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGTAACCCACCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                   <  1:814034/146‑1 (MQ=255)
                                     cGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                             <  2:606627/151‑1 (MQ=255)
                                     cGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                             >  1:144554/1‑151 (MQ=255)
                                     cGCCAGCCAGCCCATGAAGTGATCGAGGCTGTCCATCATTTTGTCGTAACCCAGCACGTCGTCCATCGGCGGTGCTGTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                             <  2:1458083/151‑1 (MQ=255)
                                         acccagTCCGTGAGGTGCTCGAGCCTTCCCACCGTTTTGTTGTAATCCAGCACGTGGTCCACCAGCGGTTCTTTTTTCGGCCCGCCCGGATTCTCCACCTCCGCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                          <  1:467063/148‑1 (MQ=255)
                                           ccagTCCATGAGGTGATCGGGGTTTTCCTTCACTTTGTCGTAATCCAGCCCGTCGCCCATCAGCGGTGTTGTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                         <  1:1937508/149‑1 (MQ=255)
                                           ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                       >  2:469445/1‑151 (MQ=255)
                                           ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                       <  2:868498/151‑1 (MQ=255)
                                            cagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcgc                                                                                                      >  2:399996/1‑151 (MQ=255)
                                             agTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                     >  1:1700968/1‑151 (MQ=255)
                                              gcccATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCACTTCGCCCATCAGCGGTGCTTTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                       <  2:1234499/146‑1 (MQ=255)
                                                ccATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                                                                                  >  2:662541/1‑151 (MQ=255)
                                                    gAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                               >  1:461426/1‑150 (MQ=255)
                                                     aaGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                              >  2:1794374/1‑150 (MQ=255)
                                                     aaGTGAGCGGGGCTGCCCATCACTTTGTCGTAATCCAGCCCGTCCTCCACCGGCGGTGTTTTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                             <  1:815419/151‑1 (MQ=255)
                                                       gggATCGAGGCTGTCCATCATTTTGTCGTAACCCACCCCGTCGCCCACCGGCGGTGCTGTTTTCGGCCCGACCAGAATCTTCAGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                                                                                            <  1:1294412/148‑1 (MQ=255)
                                                         gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGCGCACCAAAg                                                                                           >  1:1259732/1‑149 (MQ=255)
                                                          accGAGGCTGTCCATCCTTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                          <  2:114122/147‑1 (MQ=255)
                                                          aTCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAAc                                                                                        >  2:869272/1‑151 (MQ=255)
                                                              aGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                     >  2:707039/1‑150 (MQ=255)
                                                               ggCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                                                   >  1:512848/1‑151 (MQ=255)
                                                                gCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                                                  >  2:1345095/1‑151 (MQ=255)
                                                                    tCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCt                                                                              >  2:1703134/1‑151 (MQ=255)
                                                                     ccctccttttGTGGTAATCCAGCCCGTCTCCAATCAGGGGTGCTGTTTGCGGCCCGCCCGAAATCTTCGCCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                             <  2:1160210/144‑1 (MQ=255)
                                                                     ccATCTTTTTGTCTTAATCCAGCACGTCGTCCACCAGCGGTGCTGTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGc                                                                               <  2:1258414/149‑1 (MQ=255)
                                                                     ccATCCTTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGTTGTTTGCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                                 <  2:118460/147‑1 (MQ=255)
                                                                     ccATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                             >  1:387824/1‑151 (MQ=255)
                                                                      catctttttGTCGAAACCCAGCACTTCGTCCATCAGGGGTGTTGTTTTCGCCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                <  1:906812/141‑1 (MQ=255)
                                                                            ttggTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTGCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                <  1:1461963/138‑1 (MQ=255)
                                                                             ttGTCGTAACCCAGCACGTCGTCCATCAGGGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                      <  2:608965/150‑1 (MQ=255)
                                                                             ttGCCGTAATCCACCCCGTCGCCCATCAGCGTTGTTTTTTTCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                      <  1:1701672/150‑1 (MQ=255)
                                                                              tGTCGTAACCCAGCACGTCTTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCa                                                                     <  2:1365511/150‑1 (MQ=255)
                                                                                 cGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                 >  2:881271/1‑151 (MQ=255)
                                                                                    aaTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                              >  2:48469/1‑151 (MQ=255)
                                                                                       ccgccccgtcccccaccaGCGTTGTTTTTTTCGGCCCGCCTTGAATTTTCAGCTCCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTg                                                               <  1:1463176/135‑1 (MQ=255)
                                                                                       ccaccaCGTCGCCCATCAGCGGTCTTGTTTTCGGCCCGCCCTGATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                           <  1:1745772/151‑1 (MQ=255)
                                                                                       ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                           <  2:1762797/151‑1 (MQ=255)
                                                                                            aCGTCGTCCATCAGCGGGGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGCTTTCGCCAGGGTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgc                                                      >  2:1902387/1‑151 (MQ=255)
                                                                                              gtcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgc                                                      >  2:912195/1‑149 (MQ=255)
                                                                                                cgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                     <  2:365993/148‑1 (MQ=255)
                                                                                                 gtcCATCAGGGGGGCTTTTTTCGGCCCGCCCTGAATCTTCACCTTCCCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagca                                                  <  2:1646458/150‑1 (MQ=255)
                                                                                                   ccaccggCGGTGTTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACg                                                       <  1:1301395/137‑1 (MQ=255)
                                                                                                         gCGGTGCTTTTTCCGCCCCGCCCTGATTCTCCATCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                         <  2:1273607/151‑1 (MQ=255)
                                                                                                         gCGGTGCTGTTTTCGCCCCGACCTAAATCTTCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                         <  2:683510/151‑1 (MQ=255)
                                                                                                          cGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                        >  1:1009163/1‑151 (MQ=255)
                                                                                                           ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                       >  1:27353/1‑151 (MQ=255)
                                                                                                           ggTGCTGTTTGCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                          <  1:878875/148‑1 (MQ=255)
                                                                                                           ggTCCTGTTTTCGGCCGGCCCGGAATCTTCATCTTCCCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTCTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                        <  2:1316955/150‑1 (MQ=255)
                                                                                                            ggtcttttttCGCCCCGACCTGAATCTTCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                           <  1:439169/140‑1 (MQ=255)
                                                                                                            gTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                       >  1:1142592/1‑150 (MQ=255)
                                                                                                              gCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAt                                    >  1:1158990/1‑151 (MQ=255)
                                                                                                                tGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                                  >  2:1525006/1‑151 (MQ=255)
                                                                                                                  tttgCGGCCCGCCCTGAATCTTCAGCTCCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                                   <  1:802434/148‑1 (MQ=255)
                                                                                                                    ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                              >  1:1359591/1‑151 (MQ=255)
                                                                                                                       ggCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                           >  2:1700239/1‑151 (MQ=255)
                                                                                                                            gACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                        >  2:1469817/1‑149 (MQ=255)
                                                                                                                            gACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                        >  1:799343/1‑149 (MQ=255)
                                                                                                                              cctgtttctcCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAg                      <  1:1652051/143‑1 (MQ=255)
                                                                                                                                    ttttCAGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc               <  1:143855/148‑1 (MQ=255)
                                                                                                                                    tCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg              >  2:1589251/1‑151 (MQ=255)
                                                                                                                                      ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCa                   <  1:408669/144‑1 (MQ=255)
                                                                                                                                       tCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc               <  1:1919821/147‑1 (MQ=255)
                                                                                                                                        cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc          >  1:1570719/1‑151 (MQ=255)
                                                                                                                                        cAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc          <  1:1788390/151‑1 (MQ=255)
                                                                                                                                         aGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCa            <  1:746281/148‑1 (MQ=255)
                                                                                                                                            ttCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg        <  2:1813121/149‑1 (MQ=255)
                                                                                                                                             tctcGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa         <  1:1037881/147‑1 (MQ=255)
                                                                                                                                              ctcGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAt      <  2:1243468/149‑1 (MQ=255)
                                                                                                                                                cGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc  >  1:1866300/1‑151 (MQ=255)
                                                                                                                                                 gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc  >  2:1734131/1‑150 (MQ=255)
                                                                                                                                                 gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc  >  2:1921703/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTC  >  NC_000913/3260906‑3261200

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