Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT42.2% 67.6 / 41.2 71S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (6/35);  new base T (30/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604196‑4604494
                                                                                                                                                    |                                                                                                                                                        
cTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttttt                                                                                                                                                        >  2:1193597/1‑151 (MQ=255)
 tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGttttttttt                                                                                                                                                       >  2:1257898/1‑151 (MQ=255)
  cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttttttt                                                                                                                                                      >  1:678175/1‑151 (MQ=255)
  cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                                        <  1:752595/149‑1 (MQ=255)
   aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                     <  2:588817/151‑1 (MQ=255)
     aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                                   >  1:437664/1‑149 (MQ=255)
     aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                                   >  2:781745/1‑149 (MQ=255)
        aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTtttttttttttttt                                                                                                                                                   >  1:143658/1‑146 (MQ=255)
         aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTcc                                                                                                                                                >  1:1161927/1‑150 (MQ=255)
           aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTtttttttttttttttcctt                                                                                                                                              >  2:499866/1‑143 (MQ=255)
           aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                                                                             <  2:683475/151‑1 (MQ=255)
               gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                         <  2:963453/151‑1 (MQ=255)
                 cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGaa                                                                                                                                         >  1:295906/1‑149 (MQ=255)
                  aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCGtttttttttttttcggc                                                                                                                                                >  1:1285414/1‑136 (MQ=255)
                    gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                     <  2:1103494/150‑1 (MQ=255)
                     agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACACTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaggag                                                                                                                                        >  1:943486/1‑146 (MQ=255)
                        agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCAGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTCTCCTTAGAAGTg                                                                                                                                      >  1:1143296/1‑145 (MQ=255)
                          taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCag                                                                                                                               >  2:226392/1‑150 (MQ=255)
                           aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCTTTTTTTCTTTTTTTGCCATGTCAAGAGGgcgc                                                                                                                                 >  1:1094081/1‑147 (MQ=255)
                             cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                           <  2:195199/151‑1 (MQ=255)
                              ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAgg                                                                                                                          >  1:152267/1‑151 (MQ=255)
                                 aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                       <  2:103474/151‑1 (MQ=255)
                                     agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                   <  2:712689/151‑1 (MQ=255)
                                      gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                                                                                                  <  2:1015967/151‑1 (MQ=255)
                                       agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                 <  2:1296375/151‑1 (MQ=255)
                                       agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                 <  2:27770/151‑1 (MQ=255)
                                            aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGGGGGTGGCGGGAGGGGGTATCTTGcc                                                                                                              >  2:806269/1‑149 (MQ=255)
                                              tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAggggattttttttctta                                                                                                           >  1:423977/1‑136 (MQ=255)
                                               cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac‑t                                                                                                        <  2:633995/151‑1 (MQ=255)
                                                        tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCCTGTGGAGTGAGCGCCGGGAGGGATTATCTTTCCTGAC‑TGGACt                                                                                                   >  1:692801/1‑147 (MQ=255)
                                                           gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGAC‑TGAACTGGC‑GAtt                                                                                           >  2:1253412/1‑151 (MQ=255)
                                                                  ccgtataAGCTCCATACCTATACCGGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCATTTTTTCTTTTTCTGCCATGATAAGTGAGCGCAGATAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTc                                                                                    <  1:1252471/148‑1 (MQ=255)
                                                                  ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCCTGAC‑TGAACTGGC‑GATTAAGttt                                                                                     >  2:561498/1‑150 (MQ=255)
                                                                     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTGAC‑TGAAATGGC‑GGTTAAGTTTcc                                                                                   >  2:202880/1‑149 (MQ=255)
                                                                     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCcg                                                                                  <  1:658561/150‑1 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCcgc                                                                                 <  2:346635/150‑1 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTTCCACGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGAC‑TGAACTGGC‑GATTAAATTTTcgcg                                                                                >  1:743833/1‑151 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTTTGCCCTGGGAGGTGAGGGCAGGGGGGGGGTTTGTTGCCCGGC‑TGAACTTGG‑GATTAAGgtttc                                                                                   >  2:530081/1‑146 (MQ=255)
                                                                       taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGAC‑TTAACTGGC‑GATTAAGATTCCGCGt                                                                               >  2:1052624/1‑151 (MQ=255)
                                                                        aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTa                                                                              <  1:403123/151‑1 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCATTTTTCTTTTTCTGCCATTAGAAGAGTGGGCAAAGGGGGGTTATCTTTCCTGACTTAAATTGGG‑GGTTAATTTTCcgctacc                                                                             >  2:1152520/1‑146 (MQ=255)
                                                                             ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCgg                                                                         <  2:1164370/151‑1 (MQ=255)
                                                                              cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGc                                                                        <  2:1179413/151‑1 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATTTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTTCGCGTACCCGGCg                                                                       >  2:764225/1‑149 (MQ=255)
                                                                                  ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTCAGAAGCGAGCGCAGAGGGGGATTATGTTGCCTGCC‑CGTAATGGG‑GAATAAAGTTTCTCGTGCCCgg                                                                         >  2:769795/1‑146 (MQ=255)
                                                                                  ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGtt                                                                     <  1:777730/150‑1 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGGGCAGAGAGGGATTTTTTTTCCTGAC‑TGTACTTGC‑GATTAAGTTTTCGCGTACCCGGCGGTTTGt                                                                 >  2:797305/1‑151 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTa                                                                   <  1:465704/149‑1 (MQ=255)
                                                                                      taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGAGGGCGCAGGGAGGGGTTATCTTGCCCGAC‑TGAACTGGC‑GACTAAATTTCCGCGTACCCCGCGTTTAGtt                                                                >  1:933362/1‑151 (MQ=255)
                                                                                       aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAAGGAGGGCCGGGGGGGAGTATCTTGCCTGAC‑TGAAATGGC‑GATTAAGCTTCCCCGTTCCCGGCGGTTAGtt                                                                >  2:1216779/1‑150 (MQ=255)
                                                                                        ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                              <  1:202876/151‑1 (MQ=255)
                                                                                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTTTGCCCTTGGAAGGGGGTGGCGGGGGGGGTTGTCTTGTCTTGC‑TGCACTGGCGGGTTAAATTT‑CTCGTACCCAGCGtttttttttt                                                             >  2:501942/1‑144 (MQ=255)
                                                                                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTTCGCGTACCCGGCGGTTTGTTAtt                                                             >  2:370856/1‑148 (MQ=255)
                                                                                             cAGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg                                                         <  1:923545/151‑1 (MQ=255)
                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAAGGAGCGCACCGGGGGATTTTATTTTCTGAT‑TGGACCGGC‑GATTAAGATTTCGCCGACCCCGCGTTTAGGTATTGCgg                                                         >  1:377858/1‑151 (MQ=255)
                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGGGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGAC‑TGAACTGGC‑GATTAAGTTTTCGCGTACCCCGCGTTTAGTTATTGCgg                                                         >  2:1168737/1‑151 (MQ=255)
                                                                                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGAC‑TGAACTGGC‑GATTAAGTTTTCGCGTACCCCGCGGTTAGTTATTGTgg                                                         >  2:605357/1‑149 (MQ=255)
                                                                                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                       <  2:829381/151‑1 (MQ=255)
                                                                                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTTTTTTTTTTCCCTGGGAAGGGAGAGCAGAGGGGGGATATTTTTCCTGAC‑TTGACTTGC‑GGTTAAGTTTTCGCGTTGCCCGGGTTTTGTTTTTGCgg                                                         >  1:1108211/1‑143 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGCGCAGAGAGGGGTTTTCTTGCCTGAC‑TGAACTGGC‑GAATAAATTTCCGCGTACCCGGCGGTTAGTTATTGCGGATGAAGATAtt                                              >  2:5640/1‑150 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                              <  1:297809/150‑1 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                             <  2:321786/151‑1 (MQ=255)
                                                                                                          aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCTTTTAGTTATTGCGGATGATGATATTGa                                            <  2:332535/151‑1 (MQ=255)
                                                                                                          aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                            <  2:530186/151‑1 (MQ=255)
                                                                                                           ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGAC‑TGAACTGGC‑GGTTAAGGTTTCGCGTACCCCGCGTTTTGTTTTTGCGGGTGATGATAtt                                              >  1:889917/1‑148 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTCTGCCCAGAGGAGTGAGGGCGGAGAGGGGTTATTTTTCCTCAC‑TGAAACGGG‑GTTTAAGTTTTCGCGGCCCCCGCGCTTTGTTATTGCGGGTGAAGATATTTAAg                                          >  2:78359/1‑150 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                         <  2:1292842/151‑1 (MQ=255)
                                                                                                                ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                      <  1:13409/151‑1 (MQ=255)
                                                                                                                 cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATTTTTCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTTTTTCGGATGATGATATTGAAGCTCgg                                     >  2:265416/1‑151 (MQ=255)
                                                                                                                  aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                    <  2:1143331/151‑1 (MQ=255)
                                                                                                                       gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                               <  2:266389/151‑1 (MQ=255)
                                                                                                                         ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  1:1011286/151‑1 (MQ=255)
                                                                                                                         ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  2:386714/151‑1 (MQ=255)
                                                                                                                           aTCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGAGAGGGCAGAGAGGGGTTATTTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTTCGCGTACCCGGCGTTTAGTTATTGAGGATGATGATATTTGAACTCGGCCGATTgggt                           >  2:985985/1‑149 (MQ=255)
                                                                                                                               cGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGAC‑TGAACTGGC‑GGTTAAGTTTTCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTGATTGGTTTAtt                       >  2:4396/1‑151 (MQ=255)
                                                                                                                               cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                           <  1:1285322/147‑1 (MQ=255)
                                                                                                                                gTTTGCGGCGGTCCTTTTTTCTTTTTCTTCCCTGAGAAGTGAGCGCAGAGGGGGATTATTTTTCCTGAC‑TGTAATTGC‑GATTTAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTGAAACTCGGCTGGTTTGTTCAtt                       >  1:253495/1‑150 (MQ=255)
                                                                                                                                   tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                    <  2:1014666/150‑1 (MQ=255)
                                                                                                                                       gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt               >  1:33041/1‑151 (MQ=255)
                                                                                                                                        cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              <  2:316904/151‑1 (MQ=255)
                                                                                                                                         ggTCATTTTTGCTTTTTCTGACATGAGGAGTGATCGCAGAGAGGGACTATCTGGCCTGAC‑TGAACTGTG‑GAGTAAGTTTCCGCGTGCCCGGCGTTTAGTTAATGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCTTTGTCGCCATCa             <  1:249020/151‑1 (MQ=255)
                                                                                                                                               tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt        <  1:428572/150‑1 (MQ=255)
                                                                                                                                               tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt        <  2:147007/150‑1 (MQ=255)
                                                                                                                                                 tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      >  1:148489/1‑150 (MQ=255)
                                                                                                                                                  ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGa     >  1:437795/1‑150 (MQ=255)
                                                                                                                                                    ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  2:767822/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAC‑TGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604196‑4604494

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