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

CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  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)
                                               cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                        <  2:633995/151‑1 (MQ=255)
                                                        tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCCTGTGGAGTGAGCGCCGGGAGGGATTATCTTTCCTGACTGGACt                                                                                                   >  1:692801/1‑147 (MQ=255)
                                                           gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGAtt                                                                                           >  2:1253412/1‑151 (MQ=255)
                                                                  ccgtataAGCTCCATACCTATACCGGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCTGTCATTTTTTCTTTTTCTGCCATGATAAGTGAGCGCAGATAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTc                                                                                    <  1:1252471/148‑1 (MQ=255)
                                                                  ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                     >  2:561498/1‑150 (MQ=255)
                                                                     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAAATGGCGGTTAAGTTTcc                                                                                   >  2:202880/1‑149 (MQ=255)
                                                                     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcg                                                                                  <  1:658561/150‑1 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                                                                                 <  2:346635/150‑1 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTTCCACGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAATTTTcgcg                                                                                >  1:743833/1‑151 (MQ=255)
                                                                      ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTTTGCCCTGGGAGGTGAGGGCAGGGGGGGGGTTTGTTGCCCGGCTGAACTTGGGATTAAGgtttc                                                                                   >  2:530081/1‑146 (MQ=255)
                                                                       taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTTAACTGGCGATTAAGATTCCGCGt                                                                               >  2:1052624/1‑151 (MQ=255)
                                                                        aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTa                                                                              <  1:403123/151‑1 (MQ=255)
                                                                           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCATTTTTCTTTTTCTGCCATTAGAAGAGTGGGCAAAGGGGGGTTATCTTTCCTGACTTAAATTGGGGGTTAATTTTCcgctacc                                                                             >  2:1152520/1‑146 (MQ=255)
                                                                             ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                         <  2:1164370/151‑1 (MQ=255)
                                                                              cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc                                                                        <  2:1179413/151‑1 (MQ=255)
                                                                                 aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATTTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCg                                                                       >  2:764225/1‑149 (MQ=255)
                                                                                  ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTCAGAAGCGAGCGCAGAGGGGGATTATGTTGCCTGCCCGTAATGGGGAATAAAGTTTCTCGTGCCCgg                                                                         >  2:769795/1‑146 (MQ=255)
                                                                                  ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt                                                                     <  1:777730/150‑1 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGGGCAGAGAGGGATTTTTTTTCCTGACTGTACTTGCGATTAAGTTTTCGCGTACCCGGCGGTTTGt                                                                 >  2:797305/1‑151 (MQ=255)
                                                                                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                                   <  1:465704/149‑1 (MQ=255)
                                                                                      taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGAGGGCGCAGGGAGGGGTTATCTTGCCCGACTGAACTGGCGACTAAATTTCCGCGTACCCCGCGTTTAGtt                                                                >  1:933362/1‑151 (MQ=255)
                                                                                       aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAAGGAGGGCCGGGGGGGAGTATCTTGCCTGACTGAAATGGCGATTAAGCTTCCCCGTTCCCGGCGGTTAGtt                                                                >  2:1216779/1‑150 (MQ=255)
                                                                                        ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                              <  1:202876/151‑1 (MQ=255)
                                                                                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTTTGCCCTTGGAAGGGGGTGGCGGGGGGGGTTGTCTTGTCTTGCTGCACTGGCGGGTTAAATTTCTCGTACCCAGCGtttttttttt                                                             >  2:501942/1‑144 (MQ=255)
                                                                                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGTTAtt                                                             >  2:370856/1‑148 (MQ=255)
                                                                                             cAGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCgg                                                         <  1:923545/151‑1 (MQ=255)
                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAAGGAGCGCACCGGGGGATTTTATTTTCTGATTGGACCGGCGATTAAGATTTCGCCGACCCCGCGTTTAGGTATTGCgg                                                         >  1:377858/1‑151 (MQ=255)
                                                                                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGGGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGTTTAGTTATTGCgg                                                         >  2:1168737/1‑151 (MQ=255)
                                                                                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGGTTAGTTATTGTgg                                                         >  2:605357/1‑149 (MQ=255)
                                                                                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                       <  2:829381/151‑1 (MQ=255)
                                                                                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTTTTTTTTTTCCCTGGGAAGGGAGAGCAGAGGGGGGATATTTTTCCTGACTTGACTTGCGGTTAAGTTTTCGCGTTGCCCGGGTTTTGTTTTTGCgg                                                         >  1:1108211/1‑143 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGCGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGAATAAATTTCCGCGTACCCGGCGGTTAGTTATTGCGGATGAAGATAtt                                              >  2:5640/1‑150 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                              <  1:297809/150‑1 (MQ=255)
                                                                                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                             <  2:321786/151‑1 (MQ=255)
                                                                                                          aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCTTTTAGTTATTGCGGATGATGATATTGa                                            <  2:332535/151‑1 (MQ=255)
                                                                                                          aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                            <  2:530186/151‑1 (MQ=255)
                                                                                                           ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTTTCGCGTACCCCGCGTTTTGTTTTTGCGGGTGATGATAtt                                              >  1:889917/1‑148 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTCTGCCCAGAGGAGTGAGGGCGGAGAGGGGTTATTTTTCCTCACTGAAACGGGGTTTAAGTTTTCGCGGCCCCCGCGCTTTGTTATTGCGGGTGAAGATATTTAAg                                          >  2:78359/1‑150 (MQ=255)
                                                                                                             aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc                                         <  2:1292842/151‑1 (MQ=255)
                                                                                                                ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                      <  1:13409/151‑1 (MQ=255)
                                                                                                                 cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATTTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTTTTTCGGATGATGATATTGAAGCTCgg                                     >  2:265416/1‑151 (MQ=255)
                                                                                                                  aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                    <  2:1143331/151‑1 (MQ=255)
                                                                                                                       gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                               <  2:266389/151‑1 (MQ=255)
                                                                                                                         ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  1:1011286/151‑1 (MQ=255)
                                                                                                                         ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  2:386714/151‑1 (MQ=255)
                                                                                                                           aTCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGAGAGGGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGTTTAGTTATTGAGGATGATGATATTTGAACTCGGCCGATTgggt                           >  2:985985/1‑149 (MQ=255)
                                                                                                                               cGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTGATTGGTTTAtt                       >  2:4396/1‑151 (MQ=255)
                                                                                                                               cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                           <  1:1285322/147‑1 (MQ=255)
                                                                                                                                gTTTGCGGCGGTCCTTTTTTCTTTTTCTTCCCTGAGAAGTGAGCGCAGAGGGGGATTATTTTTCCTGACTGTAATTGCGATTTAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTGAAACTCGGCTGGTTTGTTCAtt                       >  1:253495/1‑150 (MQ=255)
                                                                                                                                   tGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTc                    <  2:1014666/150‑1 (MQ=255)
                                                                                                                                       gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt               >  1:33041/1‑151 (MQ=255)
                                                                                                                                        cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              <  2:316904/151‑1 (MQ=255)
                                                                                                                                         ggTCATTTTTGCTTTTTCTGACATGAGGAGTGATCGCAGAGAGGGACTATCTGGCCTGACTGAACTGTGGAGTAAGTTTCCGCGTGCCCGGCGTTTAGTTAATGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCTTTGTCGCCATCa             <  1:249020/151‑1 (MQ=255)
                                                                                                                                               tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt        <  1:428572/150‑1 (MQ=255)
                                                                                                                                               tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt        <  2:147007/150‑1 (MQ=255)
                                                                                                                                                 tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      >  1:148489/1‑150 (MQ=255)
                                                                                                                                                  ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGa     >  1:437795/1‑150 (MQ=255)
                                                                                                                                                    ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  2:767822/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
CTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  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: