Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC45.1% 100.8 / 33.0 82V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (37/8);  new base C (0/37);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 1.13e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAT  >  NC_000913/3260907‑3261196
                                                                                                                                                     |                                                                                                                                             
gTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcc                                                                                                                                             <  2:1077319/151‑1 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCATTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAccc                                                                                                                                            >  1:1015136/1‑150 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                           <  2:568021/151‑1 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                           >  2:56561/1‑151 (MQ=255)
   gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                           <  1:19715/150‑1 (MQ=255)
      aaTGATATTCAGCGCGCTGATGTACTGCCCCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgccg                                                                                                                                       <  2:532125/151‑1 (MQ=255)
           tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                  >  2:206798/1‑151 (MQ=255)
           tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                  >  2:854747/1‑151 (MQ=255)
              tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                               >  1:733609/1‑151 (MQ=255)
                   ggggtgatttcttGCACCGCCACCCAGTCCATGAAGTGATCGAGGCTGTCCACCTCTTTGTCGTAATCCAGCACGCCGCCCACCAGCGGTGCTTTTTTCGCCCCGACCTGACTCTTCACCTTCTCCCCCCCCCCGCCGTTAAtt                                                                                                                                 <  1:64159/132‑1 (MQ=255)
                     gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                        >  2:60131/1‑151 (MQ=255)
                     gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                        >  2:251567/1‑151 (MQ=255)
                        gATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcagc                                                                                                                     >  1:156935/1‑151 (MQ=255)
                         ttgtagggccccGCCAGCCGGTCCAGGAAGTGATGGAGGCGGTCCATCTCTTTGTCGTATCCCAGCCCTTCGTCCATCAGGGGTGCTTTTTTCGGCCCGACCGAAATCTTCAGCTTCTCTTCCCCCCCGCCGTTAATTGCGTAGagca                                                                                                                       <  1:76856/138‑1 (MQ=255)
                          tttacgccccccccacccagTCCATGAAGTGATCGGGGCTGCCCATCACTTTGTCGTAATCCAGCCCGTCGTCCATCGGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCTCCCCCCCCGCCGTTAATTGCGTAGagcagc                                                                                                                     <  2:1210989/137‑1 (MQ=255)
                            tACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                 >  1:361174/1‑151 (MQ=255)
                              cgcccccGCCAGCCAGCCCATGAAGGGATGGGGGCTGTCCATCTTTTTGTCGAAACCCACCCCGTCGTCCATCGGGGGTGTTGTTTTGGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                 <  1:1140756/144‑1 (MQ=255)
                               tGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                              >  2:1007587/1‑151 (MQ=255)
                               ccccccGCCAGCCAGTCAATGAAGTGATCGAGGCTGTCCTTCATTTTGTCTTACCCCACCCCGTCGCCCACCGGCGGTGTTTTTTTGGGCCCGACCGAATTCTACATCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                               <  1:555921/146‑1 (MQ=255)
                                gCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                                                                             >  1:338291/1‑151 (MQ=255)
                                 cacccccccccagTCCATGAAGTGATCGAGGCTGCCCATCTCTTTGTCGTAACCCAGCCCGTCGCCCACCCGGGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                              <  1:472526/140‑1 (MQ=255)
                                   ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCCTCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                            >  2:1200053/1‑149 (MQ=255)
                                   ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                            >  2:602029/1‑149 (MQ=255)
                                   ccGCCAGCCAGTCAATGAAGTGATCGGGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTCAGTTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                                                                                          <  2:338300/151‑1 (MQ=255)
                                    ccccacccagCCCATGAGGTGACGGGGGCTGTCCATCCTTTTGTCGTAACCCGCCCCGTCGTCCATCAGGGGTGCTTTTTCCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                         <  1:1052601/149‑1 (MQ=255)
                                      cccgccatTCGATAAGGTGATCGCGGCTCTCCTTCGCTTTGTCGCAACCCAGCTCGTCGCCCATCAGGGTTGTTGTTTCGGCCCCACCCTAAATCTTCCGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                         <  1:93840/146‑1 (MQ=255)
                                      ccagccactcgtgaaaGTGTTGGGGGCTGCCCCTCTCTTTTTCGTAACCCAGCCCTTCTCCCATCCCGGGTGTTTTTTCCGGCCCAACAGGAATTTCAAGTTGCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCa                                                                                                           <  2:612366/133‑1 (MQ=255)
                                       cagccagTCCATGAAGTGATCGAGGCTGTCCATCTTTTTGTCGTAATCCAGCCCGTCTCCCACCGGGGGTGCTGTTTTCGGCCGTACTTGATTCTTCATCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                  <  1:913998/139‑1 (MQ=255)
                                       cagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                       >  2:124933/1‑150 (MQ=255)
                                       cacccagccctagAAGGGATCGGGCCTGTCCCCCATTTTGTCGTAACCCAGCCCGTGCCCCATCAGCGGTCCTGTTTTCGGCGCGACGTGATTCTTCGGTTCCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                  <  1:346421/127‑1 (MQ=255)
                                         gccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTTCGCCCGACCTGAATCTTCAGCTTCCCGTCCACCTCGCCGTTAATTGCGTAGAGCAGCGCTTTCGCGAGGTTAg                                                                                                     >  1:1008580/1‑150 (MQ=255)
                                               ccATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                                                                              >  1:700898/1‑151 (MQ=255)
                                                 aggaggTGATGGAGGCTGACCATCTCTTGGTCGTACCCACGCGGTTCGCCCTTGGGGGCTGCTTTTTCGGGCCCGACCAGAAGCTGCAGCTCCTGGGCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTg                                                                                               <  1:903976/143‑1 (MQ=255)
                                                    aaGTGATCGAGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGTCCATCAGGGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                          <  2:435099/150‑1 (MQ=255)
                                                    aaGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                          >  1:281027/1‑150 (MQ=255)
                                                     aGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                         >  1:81516/1‑150 (MQ=255)
                                                       tGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                         >  2:1090548/1‑148 (MQ=255)
                                                         accGAGGCTGCCCATCCTTTTGCGTTACCCCACCACGTCGTCCATCGGCGGTGTTGTTTTCGGCCCGACGTGAATCTTCTCCTTCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                          <  1:682394/143‑1 (MQ=255)
                                                           cGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                                                  >  1:274509/1‑151 (MQ=255)
                                                            gAGGCTGTCCATCACTTTGTCGTAACCCAGCCCGTCGCCCATCAGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTAAGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                                                                                     <  2:583304/147‑1 (MQ=255)
                                                             aGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                 >  1:693649/1‑150 (MQ=255)
                                                              ggCTTTCCATCATTTTGTCGTAACCCAGCCCGTCGTCCACCAGGGGTGTTTTTTTCGGCCCGCCTTAATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                                                  <  2:1179701/148‑1 (MQ=255)
                                                                 tgcccATCTCTTTTTCGTAATCCAGCCCGTCTTCCATCAGCGGTGCTTTTTTCGGCCCGACCTAAATCTTCAGCTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                                               <  2:254673/145‑1 (MQ=255)
                                                                 tGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGAGCACCAAAGAACTGCATTTg                                                                            >  2:624569/1‑151 (MQ=255)
                                                                 tGTCAATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCTGGCTGAACTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGCTTTCGCCAGGTTAGCGCGGGCACCAAAGAACTGCATTTg                                                                            >  1:199434/1‑151 (MQ=255)
                                                                   tCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCt                                                                          >  2:1086361/1‑151 (MQ=255)
                                                                    ccATCTCTTTGTCGTAACCCAGCACGTCGTCCATCGGCGGTGTTTTTTCCGCCCCGCCCGGAATCTCCGGCTTCCCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                         <  2:1082372/151‑1 (MQ=255)
                                                                    ccATCCTTTTTTCGTAACCCAGCCCGTCGCCCATCGGCGGTGCTTTTTCCGCCCCGCCCTGAATCTTCATCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                             <  1:1139079/147‑1 (MQ=255)
                                                                       tCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                      >  2:65889/1‑151 (MQ=255)
                                                                       tCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                      >  2:1003253/1‑151 (MQ=255)
                                                                           ttttTCGTAACCCAGCCCGTCGCCCTTCAGCGGTGCTTTTTTCGCCCCGCCCTAATTCTCCTTTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                     <  1:1195236/148‑1 (MQ=255)
                                                                           tttGTCGTAACCCAGCCCGCGGTCCACCAGGGGTGCTTTTTCCGCCCCGCCCTGAATCTTCGGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAAt                                                                   <  2:574414/150‑1 (MQ=255)
                                                                                cgaaccccAGCACGTCGTCCATCAGCGGGTTTTTTTTCGCCCCGACCTAAATCTCCATCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCa                                                                 <  1:639594/141‑1 (MQ=255)
                                                                                cGTAATCCAGCACGTCGTCCATCAGCGGTGGTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                             <  1:1076500/151‑1 (MQ=255)
                                                                                cGTAACCCAGCACGTCGTCCATCGGCGGTGCTGTTTCCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCa                                                              <  2:1151173/150‑1 (MQ=255)
                                                                                   aacccaccaCGTCGTCCATCGGCGGGGTTGTTTGCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                          <  2:715218/148‑1 (MQ=255)
                                                                                   aacccAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt                                                            <  2:298808/146‑1 (MQ=255)
                                                                                     ccatgcccgtcgtCCTTCAGCGGTGTTGTTTCCGCCCCGCCCTAACTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                        <  1:342701/144‑1 (MQ=255)
                                                                                      ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                         <  2:162996/149‑1 (MQ=255)
                                                                                         gcccgtcgcCCACCAGGGGGGTTTTTTTGGCCCCGACGTGATTCTTCATCTCCCCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGa                                                     <  1:544281/147‑1 (MQ=255)
                                                                                          cACGTCGTCCATCGGCGGTGCTTTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACg                                                   <  2:225974/151‑1 (MQ=255)
                                                                                             gtcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgc                                                  >  1:908779/1‑149 (MQ=255)
                                                                                                gtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGATTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                >  2:867965/1‑148 (MQ=255)
                                                                                                 tcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                            >  2:1120054/1‑151 (MQ=255)
                                                                                                  ccATCAGCGTTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGcagcag                                             <  2:1058341/149‑1 (MQ=255)
                                                                                                      cAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                       >  1:788834/1‑151 (MQ=255)
                                                                                                        gCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGTTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCaa                                         <  2:779430/147‑1 (MQ=255)
                                                                                                         cGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                    >  1:319874/1‑151 (MQ=255)
                                                                                                                 ttttCGCCCCGCCCTAAATTTTCAGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                               <  2:980952/148‑1 (MQ=255)
                                                                                                                 tttcCGGCCCGACCTGATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                               <  2:848337/148‑1 (MQ=255)
                                                                                                                   ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                          >  2:633928/1‑151 (MQ=255)
                                                                                                                      ggCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                       >  2:368738/1‑151 (MQ=255)
                                                                                                                      ggCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                       >  1:859982/1‑151 (MQ=255)
                                                                                                                      gccccGACCTGAATCTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                         <  2:1200205/147‑1 (MQ=255)
                                                                                                                        cccGACCGGAATCTTCACCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                        <  1:1116493/148‑1 (MQ=255)
                                                                                                                               cgaaaTCTTCAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCa               <  2:481300/148‑1 (MQ=255)
                                                                                                                                  aatctcccccTTCTCTCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGt             <  1:461804/141‑1 (MQ=255)
                                                                                                                                          gCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc      >  2:440020/1‑149 (MQ=255)
                                                                                                                                          gCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc      >  2:1077324/1‑149 (MQ=255)
                                                                                                                                          gCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc      >  1:401021/1‑149 (MQ=255)
                                                                                                                                            ttCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCa        <  2:281034/145‑1 (MQ=255)
                                                                                                                                             tctcGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAt  <  1:402250/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGAT  >  NC_000913/3260907‑3261196

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