Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC41.2% 98.2 / 18.4 69V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (33/7);  new base C (0/28);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.98e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260906‑3261205
                                                                                                                                                       |                                                                                                                                                      
tGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCAc                                                                                                                                                       >  2:631930/1‑151 (MQ=255)
   aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCAcccc                                                                                                                                                    >  2:1249227/1‑151 (MQ=255)
    gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCccg                                                                                                                                                   <  1:1195865/151‑1 (MQ=255)
     tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCccgc                                                                                                                                                  >  1:576883/1‑151 (MQ=255)
      gAATGATATTCAGCGCGCTGATGTACCGCAACGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCccg                                                                                                                                                   <  1:504700/149‑1 (MQ=255)
        aTGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGt                                                                                                                                               >  1:1128092/1‑151 (MQ=255)
            tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                           >  1:90866/1‑151 (MQ=255)
             atTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAAtt                                                                                                                                          >  2:163710/1‑151 (MQ=255)
                   cgcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  2:1068242/1‑149 (MQ=255)
                    gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  1:987820/1‑148 (MQ=255)
                        tgttatcctcccccGCCAGCAAGTCCATTAAGTAATTGAGGTTGTCCATCACTTTGTTGTAACCCACCACGTGTTCCACCGGCGGTGTTGTTTTCGG‑CCCGTCTTAAA‑TCTTAATCTTCCCTTCCCCCCCGCCGTTAATTGCGTAGagc                                                                                                                                 <  2:1064468/137‑1 (MQ=255)
                           tttcctcccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGG‑CCCGCCCTAAA‑TCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTagag                                                                                                                                  <  2:423613/136‑1 (MQ=255)
                           tGTACCGCCCCGCCAGCCAGTCCATGAAGTGATCGGGGCTGTCCATCATTTTGTCGTAACCCAGCACGTCGTCCACCGGCGGTGTTGTTTTCGG‑CCCAACCTAAA‑TCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGagcagc                                                                                                                              <  2:334469/149‑1 (MQ=255)
                            gtccgcccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGAAATCCAGCCCGTCTTCCTCCAGCGTTGTTGTTTTCGC‑CCCGACCGGAC‑TCTTGAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTagag                                                                                                                                  <  1:907394/136‑1 (MQ=255)
                            gTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           >  1:110470/1‑151 (MQ=255)
                             tACTGCCCCGCCAGCCAGTCCATGAAGTGACCGAGGCTGCCCATCACTTTGTCTTAATCCAGCCCGTCTCCCATCAGCGGTGCT‑TTTTCGGCCCCGACCGGAA‑TCTCCAGTTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           <  2:581189/150‑1 (MQ=255)
                                  cccccccacccggTCCATGAAGTGTTCGGGGCTGTCCATCTCTTTGTCGTAAACCACCCCGTGGCCAATCAGCGGTGCTTTTTCCGG‑CCCGGCCTGAC‑TCTGCGGCTTCTGCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                                                                                      <  2:985643/145‑1 (MQ=255)
                                  ccccGCCAGCCAGCCCATGAAGTGATCGGGGCTGTCCATCTCTTTGTCGTAACCCACCCCGTCCCCCATCAGCGGTGCTTTTTGCGC‑CCCGACCTAAA‑TCTTCAGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                        <  1:674693/146‑1 (MQ=255)
                                  cACCCCCAGCCGGCCCATGAAGTGATCGAGGCTGTCCATCCCTTTGTCGTAACCCACCCCGTCGCCCATCAGGGGTTCTGTTTCGGC‑CCCGACCTGAA‑TCTTCCGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                                     <  2:414114/151‑1 (MQ=255)
                                           ccagTCCATGAGGTGATCGAGGCTGCCCATCTCTTTGTCGTAACCCAGCACGTAGTCCACCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                              <  2:751337/149‑1 (MQ=255)
                                           ccagTCAATGAAGTGATCGAGGCTGTCCTTCACTTTTTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTCCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  2:678614/150‑1 (MQ=255)
                                              gTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                          >  2:1051377/1‑150 (MQ=255)
                                                ccagaaaGTGATCGGCGCTGTCCATCTCTTTCTCGAAATCCGGCCCGTCCTCCATCAGCGTTGTTGTTTCCGG‑CCCGCCCTTAA‑TCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTg                                                                                                        <  1:120433/145‑1 (MQ=255)
                                                 cATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                                                                                       >  2:810499/1‑150 (MQ=255)
                                                     aaGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGT‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                    >  1:21331/1‑149 (MQ=255)
                                                          aTCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                                                                                              >  1:1296796/1‑150 (MQ=255)
                                                          aTCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAAc                                                                                             >  1:689454/1‑151 (MQ=255)
                                                              aGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                          >  2:386289/1‑150 (MQ=255)
                                                                 ctctctatcgttttGTGGAAAGCCAGCACTCCCTCCACCAGCGGTGTTTTTTGCGG‑CCCTACATGAA‑GCGTCGGTTTTTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                                                         <  2:627024/137‑1 (MQ=255)
                                                                         cACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                               >  2:949368/1‑150 (MQ=255)
                                                                            ttttTCGTATCCCCGCCCGTCGTCCATCAGCGTTGTTGTTTCCGC‑CCCGCCCTAAA‑TCTTCAGTTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                  <  1:919961/144‑1 (MQ=255)
                                                                             ttGTCGTAATCCAGCACGTCGTCCATCAGCGGGGCTGTTTTCGG‑CCCGACCTGAA‑TCATCAGCTTCTCGTCCACCCCGCCGTTAATTGCGCAGAGCAGCGTTATCGCCAGGTTAGCGCGGGCACCAAAGAACTGCATTTGCTTACCAAt                                                                            >  1:131199/1‑149 (MQ=255)
                                                                              tGTCTTAATCCAGCTCGTCTTCAATCAGGGGTGCTTTTTCCGG‑CCCGACCTAAT‑TCTGCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                           <  2:436389/149‑1 (MQ=255)
                                                                              tGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAc                                                                         >  2:378199/1‑151 (MQ=255)
                                                                              tGTCGTAAGCCACCACGTCGTCCATCCGGGTTTCTTTTTTCGG‑CCCGACCAAAA‑TCTTCACTTTCGCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                                                             <  2:722932/147‑1 (MQ=255)
                                                                              tGTCGTAACCCAGCCCGTCGCCCATCAGGGGTGTTTTTTTCGC‑CCCGCCCTGAT‑TCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                           <  1:886094/149‑1 (MQ=255)
                                                                               gTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                        >  1:383723/1‑151 (MQ=255)
                                                                                 cGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACCGCATTTGCTTACCAATCACCAt                                                                      >  2:587995/1‑151 (MQ=255)
                                                                                       ccaccaCTTCTTCCATCAGGGGTGCTTTTTTCGG‑CCCGCCCTGAA‑TCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 <  1:1095500/150‑1 (MQ=255)
                                                                                       ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCCCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                                >  1:1079978/1‑151 (MQ=255)
                                                                                        cAGCACGTCGCCCACCGGCGGTGTTGTTTTCGG‑CCCGCCCTAAA‑TCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 <  1:1043654/149‑1 (MQ=255)
                                                                                           cACGTCGTCCATCAGGGGTGCTGTTTCCGG‑CCCGCCCTGAT‑TCTTCAGCTTCCCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGa                                                              <  2:687138/149‑1 (MQ=255)
                                                                                           cACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGGTAGCGCGAGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGa                                                              >  2:682142/1‑149 (MQ=255)
                                                                                                  tcCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                                     >  2:551430/1‑151 (MQ=255)
                                                                                                    cccccggggTGCTTTTTTCCG‑CCCGGCCCAAATTCTTCGCCTCCGCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                                        <  1:208442/140‑1 (MQ=255)
                                                                                                     accggCGGTGTTGTTTTCGG‑CCCGACCTGAA‑TCTTCTGCTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                   <  1:631910/146‑1 (MQ=255)
                                                                                                        aGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                               >  2:337586/1‑151 (MQ=255)
                                                                                                         gCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                              >  1:964460/1‑151 (MQ=255)
                                                                                                           ggTCCTTTTTTCGG‑CCCGCCTTGAT‑TCTTAAGTTTTGCTTCCTCTCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                            <  2:648471/151‑1 (MQ=255)
                                                                                                            gTGCTGTTTCCGC‑CCCGGCCTGAA‑TCTTCACCTTCCCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                              <  1:260629/148‑1 (MQ=255)
                                                                                                              gCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTaa                                          >  2:1142131/1‑150 (MQ=255)
                                                                                                                  ttttCGC‑CCCGACCTAAA‑TTTTCAGCTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAt                                         <  1:1048748/147‑1 (MQ=255)
                                                                                                                  ttttCGC‑CCCCCCCTCAT‑TCTTCACCTTCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgt                                      <  2:452124/150‑1 (MQ=255)
                                                                                                                  tttcCGG‑CCCGACCTGAT‑TCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                                        <  1:498453/148‑1 (MQ=255)
                                                                                                                   ttccGG‑CCCGCCTTAAT‑TCTCCAGCTTCCCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                                        <  2:497571/144‑1 (MQ=255)
                                                                                                                     tCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                                  >  1:898125/1‑151 (MQ=255)
                                                                                                                          cccGACCTAAA‑TCTTCAGCTCCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                                  <  1:1274096/147‑1 (MQ=255)
                                                                                                                              aCCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAg                           >  2:1027090/1‑150 (MQ=255)
                                                                                                                               cgtaaa‑TTTTCAGCTCCGCGTCCACGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         <  1:26244/147‑1 (MQ=255)
                                                                                                                                   aa‑TCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGt                      <  2:521943/150‑1 (MQ=255)
                                                                                                                                      tcgtagccTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCTCCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    <  1:1009068/143‑1 (MQ=255)
                                                                                                                                              ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  2:5717/1‑147 (MQ=255)
                                                                                                                                              ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  1:916499/1‑147 (MQ=255)
                                                                                                                                              ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  1:1059371/1‑147 (MQ=255)
                                                                                                                                              ttCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCAt                <  2:631995/146‑1 (MQ=255)
                                                                                                                                              ttCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   <  2:990800/143‑1 (MQ=255)
                                                                                                                                                  cGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa              <  1:1151758/144‑1 (MQ=255)
                                                                                                                                                   gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc       >  1:462775/1‑150 (MQ=255)
                                                                                                                                                      cACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAttt   >  1:529197/1‑151 (MQ=255)
                                                                                                                                                       aCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtttt  >  2:138420/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGG‑CCCGACCTGAA‑TCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260906‑3261205

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