Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC43.9% 118.8 / 10.9 68V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (32/5);  new base C (0/29);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.33e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAT  >  NC_000913/3260908‑3261196
                                                                                                                                                   |                                                                                                                                             
tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAccc                                                                                                                                            >  1:301765/1‑151 (MQ=255)
tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCAATCAGCGGTGCTGTTTTCGGCCCGAACTGAATCTTCAGCTTCTCGTACAccc                                                                                                                                            >  1:585307/1‑151 (MQ=255)
   tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgc                                                                                                                                         >  2:521601/1‑151 (MQ=255)
   tGAATGATATTAAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgc                                                                                                                                         >  2:694435/1‑151 (MQ=255)
    gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                        <  2:200340/151‑1 (MQ=255)
        gATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGtt                                                                                                                                     >  1:508413/1‑150 (MQ=255)
             tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                               >  2:66604/1‑151 (MQ=255)
                  gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                             >  1:243046/1‑148 (MQ=255)
                    gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                        >  2:769737/1‑151 (MQ=255)
                        tgttcctccACCGCCACCCAGCCCAAGAAGTGATCGAGGCTATCCAGCCCTTTGTCTTAATCCAGCACGTCCTCCATCAGCGGTCCTGTTTGCGGCCCGCCCTGAATCTTCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTagag                                                                                                                         <  1:673807/138‑1 (MQ=255)
                         tgtgCTCCACCGCCAGCCAGTCCAGGAGGGTATGGAGGCTGTCAATCCTTTTGTGGTAATCCAGCACGTCGTCCAGCAGCGGGGCTGTTTTCGGCCCGACTTGAATCTTCAGCTCCCCGTCCTCCCCGCCGTTAATTGCGTaga                                                                                                                          <  1:100158/144‑1 (MQ=255)
                         tGTACTCCCCCGCCAGCCAGCCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGCCCATCAGCGGTGCTGTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                    <  1:223503/150‑1 (MQ=255)
                              tGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                              >  1:458274/1‑151 (MQ=255)
                                ccccGCCAGCCAGTCCATGAAGTGATCGGGGCTGTCCATCCCTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                              <  1:34807/147‑1 (MQ=255)
                                cACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                            >  2:599578/1‑151 (MQ=255)
                                  ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                            >  1:711606/1‑149 (MQ=255)
                                         ccagccctggaggTGATCGAGGCTGTCCATCATTTTGTCGAAACCCAGCACGCCGTCCATCAGCGGTGTTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                    <  2:774643/138‑1 (MQ=255)
                                         ccagTCCATGAAGTGACCGGGGCTGTCCACCCTTTTGTCGTATTCCAGCACGTCGTCCATCGGCGGTGTTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCTTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                    <  1:490056/150‑1 (MQ=255)
                                         ccagCCCAAGAGGTGAGCGAGGCTGCCAACCTTTTTGTCGTACCCCAGCCCGGCCCCCACCAGCGGTGTTTTTTCCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                      <  2:370606/148‑1 (MQ=255)
                                           agTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                 >  1:122767/1‑151 (MQ=255)
                                            gTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                >  1:333150/1‑151 (MQ=255)
                                               caggaggggagcGAGGCTGTCCATCATTTTGTCGAAATCCAGCCCGCCGTCCATCAGGGTGGCTGTTTTGGGCCCGACCTGATTTTCCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTg                                                                                               <  2:500770/138‑1 (MQ=255)
                                                aTGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                            >  2:447084/1‑151 (MQ=255)
                                                aTGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                            >  1:6351/1‑151 (MQ=255)
                                                 tGAAGTGATCGAGGCTGTCCATCTTTTGGTCGTAATCCAGCACGTCGCCCATCAGCGGTGCTGTTTCCGGCCCGCCCTGAATCTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                            <  1:506845/150‑1 (MQ=255)
                                                 tGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                           >  1:569625/1‑151 (MQ=255)
                                                   aggTGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                            <  1:769158/146‑1 (MQ=255)
                                                   aaGTGATCGGGGCTTTCCATCTCTTTGTGGTAACCCAGCACGTCGTCCATCGGCGGTGGTGTTTTCGGCCCGACCTGAATCTCCAGTTTCCCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                    <  2:528964/140‑1 (MQ=255)
                                                       gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                      >  1:579280/1‑150 (MQ=255)
                                                          cGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                                                  >  1:887969/1‑151 (MQ=255)
                                                             ggCTGTCCATCTCTTTGTCTTAACCCAGCACGTCGTCCACCAGCGGTGCTTTTTTCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                                               <  1:374387/151‑1 (MQ=255)
                                                               cTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                             >  1:865867/1‑151 (MQ=255)
                                                                tGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                            >  1:432135/1‑151 (MQ=255)
                                                                           ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                  <  1:674546/150‑1 (MQ=255)
                                                                                 tAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTg                                                           >  1:236924/1‑151 (MQ=255)
                                                                                   aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCACGTTAGCGCCTGCACCAAAGAACTGCATTTGCTTACCAATCACCa                                                              >  2:332052/1‑146 (MQ=255)
                                                                                    tCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                        >  1:442286/1‑151 (MQ=255)
                                                                                     ccaccccgtcgtCCACCGGGGGTGCTTTTTTCGGCCCGACCTGAATCTTCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                          <  2:566700/142‑1 (MQ=255)
                                                                                      cAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTg                                                      >  1:164631/1‑151 (MQ=255)
                                                                                            gtcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgc                                                  >  2:748148/1‑149 (MQ=255)
                                                                                               gtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                               <  1:317530/149‑1 (MQ=255)
                                                                                               gtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                               >  2:487357/1‑149 (MQ=255)
                                                                                                tccaccggggttGCTGTTTGCGGCCCGCCCTGAATCTGCTTCTTCCCGACCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                 <  2:725622/135‑1 (MQ=255)
                                                                                                 ccATCGGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                            <  2:301723/150‑1 (MQ=255)
                                                                                                  caccAGCGGTGCTGTTTCCGGCCCGCCCTGAATCTTCAGCTTCTCACCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGcagcag                                             <  2:288078/145‑1 (MQ=255)
                                                                                                     cAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                       >  1:628131/1‑151 (MQ=255)
                                                                                                        cGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                    >  1:576971/1‑151 (MQ=255)
                                                                                                         gttgttGTTTGCGGCCGCTCCTGATTCTCCAGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                   <  2:260330/146‑1 (MQ=255)
                                                                                                         ggtggtGTTTTCGGCCCGCCCTGAATCTTCATTTTCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                   <  2:716932/151‑1 (MQ=255)
                                                                                                         ggTGCTGTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                   <  2:905220/151‑1 (MQ=255)
                                                                                                          gtggttttttGGGCCGGGCCTAAATCTTAAGCTTCTCGTCCTCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                  <  1:521693/145‑1 (MQ=255)
                                                                                                                ttttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                            >  2:925162/1‑151 (MQ=255)
                                                                                                                 tttCGGCCCGCCCTCAGTCGGCGTCTTGGCCTCCCCTCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                              <  2:915450/148‑1 (MQ=255)
                                                                                                                 tttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                           >  1:132011/1‑151 (MQ=255)
                                                                                                                  tccGGCCCGACCTGAATCTTCAGCTCCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                           <  1:195941/148‑1 (MQ=255)
                                                                                                                     ggCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                          <  1:284262/148‑1 (MQ=255)
                                                                                                                       ccggcccTGATTCTCCGGCTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                       <  2:53985/144‑1 (MQ=255)
                                                                                                                        ccGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                    <  1:59481/151‑1 (MQ=255)
                                                                                                                            cgtGATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                <  2:777887/149‑1 (MQ=255)
                                                                                                                            ccTAAATCTTCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                <  1:758932/151‑1 (MQ=255)
                                                                                                                                attCTTCAGTTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTa            <  2:348922/149‑1 (MQ=255)
                                                                                                                                aaTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                <  2:247996/147‑1 (MQ=255)
                                                                                                                                     tCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCa        <  1:138815/150‑1 (MQ=255)
                                                                                                                                      cAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCa        <  2:927435/149‑1 (MQ=255)
                                                                                                                                      cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc      >  2:844472/1‑151 (MQ=255)
                                                                                                                                          ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGCTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc      >  1:434507/1‑147 (MQ=255)
                                                                                                                                            ctcGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAt  <  1:353609/149‑1 (MQ=255)
                                                                                                                                            ctcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg    <  2:23395/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAT  >  NC_000913/3260908‑3261196

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