Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT43.8% 118.3 / 53.0 96S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (9/45);  new base T (42/0);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.51e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604198‑4604494
                                                                                                                                                   |                                                                                                                                                       
cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttttttt                                                                                                                                                     >  1:1400222/1‑151 (MQ=255)
cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttcttt                                                                                                                                                      <  1:506279/150‑1 (MQ=255)
 aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttcttttt                                                                                                                                                    <  1:1632760/151‑1 (MQ=255)
  tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttcttttt                                                                                                                                                    >  2:304373/1‑150 (MQ=255)
   aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCtttttttttttttt                                                                                                                                                  >  2:679707/1‑149 (MQ=255)
   aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttctttttc                                                                                                                                                   <  1:1092796/150‑1 (MQ=255)
    gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CAtttttttttttttt                                                                                                                                                  >  1:441505/1‑148 (MQ=255)
    gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTttttttttttct                                                                                                                                                  >  2:1370856/1‑150 (MQ=255)
      aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGcc                                                                                                                                               >  2:1125295/1‑151 (MQ=255)
      aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGG‑ACTTTTTTTTTTTTTTgtc                                                                                                                                               >  1:1609155/1‑149 (MQ=255)
      aaaaGGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CTtttttttttttttttc                                                                                                                                                >  2:1383941/1‑146 (MQ=255)
       aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CAtttttttttttttttcc                                                                                                                                               >  1:302348/1‑145 (MQ=255)
        aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCAt                                                                                                                                             <  1:1137034/151‑1 (MQ=255)
           gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATgag                                                                                                                                          <  1:221457/151‑1 (MQ=255)
           gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATga                                                                                                                                           <  1:1112438/150‑1 (MQ=255)
            cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTTTTTTTTTTTTTTgtc                                                                                                                                               >  1:1302528/1‑143 (MQ=255)
            cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGcc                                                                                                                                               >  2:666820/1‑145 (MQ=255)
               cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                       <  2:1432504/150‑1 (MQ=255)
                aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                     <  2:841920/151‑1 (MQ=255)
                aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGG‑CTttttttttttttttcaccaaagca                                                                                                                                        >  2:1437798/1‑136 (MQ=255)
                 tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGTg                                                                                                                                     >  1:90837/1‑150 (MQ=255)
                   agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                                  <  1:1343672/151‑1 (MQ=255)
                           cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagagag                                                                                                                          <  2:1360049/151‑1 (MQ=255)
                               aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAt                                                                                                                      <  1:1284458/151‑1 (MQ=255)
                                 tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGAtt                                                                                                                     <  1:1311274/150‑1 (MQ=255)
                                  gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CTTTTTTTTTTTTTTTTCCATGAGGAGGGAGGGCAGAGAGGGGtt                                                                                                                     >  2:1052101/1‑149 (MQ=255)
                                        ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCCTGAGAAGGGAGCGCGGAGGGGGGTTATCTTgt                                                                                                              >  1:204522/1‑149 (MQ=255)
                                         tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGGGGGGGGAGAGAggggatatttt                                                                                                                >  1:1631298/1‑139 (MQ=255)
                                         tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc                                                                                                             <  1:981351/150‑1 (MQ=255)
                                           ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGGGAGGGGTTTTCtt                                                                                                                >  2:1463674/1‑145 (MQ=255)
                                           ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGGGAAGGGGGGGCAGAGAGGGGTTTTCTTTCct                                                                                                            >  2:1636328/1‑149 (MQ=255)
                                           ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCAGAGAGGGGTTATCtt                                                                                                                >  2:713239/1‑145 (MQ=255)
                                           ttCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTTTCTTTCctg                                                                                                           >  2:1220054/1‑150 (MQ=255)
                                            tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctga                                                                                                          <  1:1369391/150‑1 (MQ=255)
                                                 aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGATTATTTTGTCTGACTGaa                                                                                                     >  1:215466/1‑150 (MQ=255)
                                                  aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTTTCTTGCCTGACTGAAc                                                                                                    >  1:142813/1‑150 (MQ=255)
                                                  aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                                   <  1:67175/151‑1 (MQ=255)
                                                   gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGaa                                                                                                     >  2:809852/1‑148 (MQ=255)
                                                      tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGCGGGGTTATCTTGCCTGACTTAACTGGCg                                                                                               >  1:512889/1‑151 (MQ=255)
                                                       gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCg                                                                                               <  2:208591/150‑1 (MQ=255)
                                                         gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGG‑tt                                                                                           >  2:796550/1‑151 (MQ=255)
                                                          ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCGGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCAAGAGAAGAGGGCGCAGAGAGGGGTTATCTTGCCTGATCGAAATGGGGaat                                                                                            >  2:363048/1‑148 (MQ=255)
                                                            tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTTCCCTTAGAAGGGAGCGCAGAGGGGGGTTTTTTTTCCTGGCTGAACTGGCGA‑tt                                                                                           >  1:1614824/1‑149 (MQ=255)
                                                            tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGGAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTa                                                                                          >  2:907049/1‑149 (MQ=255)
                                                              aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGGAGGGGGAGCAGGGGGGGATTATCTTGCCTTACTGAAATgggggtt                                                                                            >  2:978720/1‑141 (MQ=255)
                                                              aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGGGGGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGG‑TTaa                                                                                         >  2:1503890/1‑148 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCTTGTCGAGTTTTTGGTTTGCGGCGGT‑CAGTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTcc                                                                                   <  2:1633354/149‑1 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTcc                                                                                   <  2:302355/149‑1 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCcg                                                                                  <  1:1197221/150‑1 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCcg                                                                                  <  2:804106/150‑1 (MQ=255)
                                                                     taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCCGAGAGGGGTTATTTTTCCTGACTGAACTGGCGA‑TTTAGTTTTcgcg                                                                                >  1:682808/1‑150 (MQ=255)
                                                                      aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTa                                                                              <  1:477075/151‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGG‑TATTTTTTTTTTTTTTTCCCTGAGGAAAGAGGGCAGAGGGGGATTATTTTGCCTTACTGTCATGGCGA‑TTAATGTTTCGCGGAccc                                                                           >  1:175928/1‑151 (MQ=255)
                                                                           ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTAccc                                                                           <  1:918699/149‑1 (MQ=255)
                                                                             aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTGCCTGGGAAGGGAGGGCCGGGAGGGGTTTTCTTGCCTGGCTGAAATGGCGG‑TTAAATTTCCCCGGAccc                                                                           >  2:1142218/1‑147 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTTCCCTTAGAAGTGAGCGCCGAGAGGGATTATCTTGCCTGACTGAAATGGGGG‑TTTAGTTTCCGCGTACCCcgcg                                                                       >  2:418306/1‑150 (MQ=255)
                                                                                 cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGGAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTTAGTTTCCGCGGACCCGGCGgtt                                                                    >  1:1073365/1‑150 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGgtt                                                                    >  2:1456854/1‑148 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                 <  1:1125253/151‑1 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                 <  1:1647017/151‑1 (MQ=255)
                                                                                    taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                 <  1:318512/150‑1 (MQ=255)
                                                                                     aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                <  1:958082/150‑1 (MQ=255)
                                                                                     aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                               <  2:1671479/151‑1 (MQ=255)
                                                                                        aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                            <  1:733533/151‑1 (MQ=255)
                                                                                         gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                           <  1:1071303/151‑1 (MQ=255)
                                                                                           cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGTGGAGAGGGGGCCGAGGGGGGGTAACTTGCCTTGCTGAAATGGGGA‑TTAAGTTTTCGCGTACCCCGCGTTTAGGTATTGCgg                                                         >  1:1070122/1‑151 (MQ=255)
                                                                                            aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                        <  2:1538873/151‑1 (MQ=255)
                                                                                                gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGC‑CATTTTTTTTTTTTTTGCCCTGGGGAGGGGGCGGACAGGGGGATTTTTTTGCCTGGCTGAAAAGGCGG‑TTAAGTTTTTGCGCGCTCGGCGGTTGGCTAttttgggtg                                                      >  1:434887/1‑142 (MQ=255)
                                                                                                   aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGGAGGGGGGGCAGGGAGGGATTATTTTGCCTGACTGAACTGGGGG‑TTAAGTTTTCGCGTACCCGGCGTTTTGTTATTTCCggtgatgat                                                 >  2:627347/1‑151 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTTCCCTGAGGAGGGAGCGCAGAGAGGGGTTTTTTTTCCTGACTGAACTGGGGG‑TTAAGATTTCGCGGACCCGGCGTTTTGTTTTTGCGGATGGTGatat                                               >  1:692200/1‑151 (MQ=255)
                                                                                                     ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGata                                                <  1:1672038/150‑1 (MQ=255)
                                                                                                      tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAAGGAGCGCAGGGGGGGGTTATCTTGCCTGACTGAAATGGGGG‑TTTAGTTTTCGCGGACCCCGCGTTTTGTTATTTCGGGGGGTGATAtt                                              >  1:1188540/1‑151 (MQ=255)
                                                                                                         ggATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                           <  2:194403/151‑1 (MQ=255)
                                                                                                          gATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                          <  1:415282/151‑1 (MQ=255)
                                                                                                             gCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                       <  2:1282902/151‑1 (MQ=255)
                                                                                                              ccATGTCGCGGATCTCGTTTGCGGCGGT‑TAATTTTTTTTTTTCTTGCACGGGAAAGGAGCGCAGAGGGGGGTTTTTTTTGCTTTTTTCACTAGCGATTTAAGTTT‑TTCGTACCCGGCGGTTTTTTTTTGTGGGAGGTGGTATTGAAggccg                                      >  1:238066/1‑147 (MQ=255)
                                                                                                              ccATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                      <  1:513038/151‑1 (MQ=255)
                                                                                                               cATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                     <  2:588555/151‑1 (MQ=255)
                                                                                                               cATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                     <  1:561294/151‑1 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGT‑CCTTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGAGGGGTTATCTTTCCTGACTGAAATGGCGA‑TTAAGTTTTCGCGTACCCCGCGGTTTGTTATTGCGGGTGGTGATATTTAAGCTCGGCTg                                  >  2:667065/1‑151 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTGCCTGGCTGAACTTGCGG‑TTAAGTTTCCGCGTACCCGGGGGTTTGTTATTGCGGGTGATGATATTTAAGCTCgg                                     >  2:1550575/1‑148 (MQ=255)
                                                                                                                    cgcgGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGG‑TTAAGTTTCCGCGTACCCGGCGGTTTGTTATTGCGGGTGATGATATTGAAGCTCGGCTg                                  >  2:1025287/1‑149 (MQ=255)
                                                                                                                     gcgGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                                  <  2:294749/148‑1 (MQ=255)
                                                                                                                             cGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGGAGAGAGGGGTTATTTTTCCTGACTTAACTGGCGA‑TTTAGTTTTCGCGGACCCGGCGTTTTGTTATTTCGGGGGATGATATTGAAGCTCGGCTGGTTGGTTTAtt                       >  1:1365024/1‑151 (MQ=255)
                                                                                                                             cGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                         <  1:666797/149‑1 (MQ=255)
                                                                                                                             cGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                        <  1:825235/150‑1 (MQ=255)
                                                                                                                                   cggcggT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                 <  2:746881/151‑1 (MQ=255)
                                                                                                                                   cggcggT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                 >  2:1387532/1‑151 (MQ=255)
                                                                                                                                      cggT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              >  2:965449/1‑151 (MQ=255)
                                                                                                                                      cggT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              <  2:682832/151‑1 (MQ=255)
                                                                                                                                      cggT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              >  2:232041/1‑151 (MQ=255)
                                                                                                                                         t‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc           >  2:1657496/1‑151 (MQ=255)
                                                                                                                                         t‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc            >  2:1110318/1‑150 (MQ=255)
                                                                                                                                           cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg          <  2:1657703/151‑1 (MQ=255)
                                                                                                                                            aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc         >  1:1483372/1‑151 (MQ=255)
                                                                                                                                               ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg      >  2:1329922/1‑151 (MQ=255)
                                                                                                                                               ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg       <  2:792770/150‑1 (MQ=255)
                                                                                                                                               ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg       <  1:614971/150‑1 (MQ=255)
                                                                                                                                                  tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACg   >  2:357254/1‑151 (MQ=255)
                                                                                                                                                   ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  2:1108053/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGA‑TTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604198‑4604494

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