Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT41.6% 144.0 / 94.9 137S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (14/66);  new base T (57/0);  total (71/66)
Fisher's exact test for biased strand distribution p-value = 1.48e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 161 total aligned reads displayed.

ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA  >  NC_000913/4604207‑4604483
                                                                                                                                         |                                                                                                                                            
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                            >  1:931157/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                            >  1:266914/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                            >  2:2273825/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCAATGTCGCGGATCTCGTTTGCGGCGGTTATtttttttt                                                                                                                                            >  2:563363/1‑140 (MQ=255)
  gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttctttt                                                                                                                                          >  1:1025728/1‑140 (MQ=255)
  gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                          >  1:1657927/1‑140 (MQ=255)
  gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                          >  2:619079/1‑140 (MQ=255)
  gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                          >  1:680135/1‑140 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                       <  2:1126412/140‑1 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                       <  2:1559034/140‑1 (MQ=255)
               taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCACACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                             >  1:1363452/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGTGAg                                                                                                                        >  2:166913/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAg                                                                                                                        >  2:1929261/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAg                                                                                                                        >  1:122508/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                        >  1:1586592/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAg                                                                                                                        >  1:1002060/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTAAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAg                                                                                                                        >  1:298596/1‑140 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCC‑‑ACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgcg                                                                                                                      >  1:2229956/1‑140 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgcg                                                                                                                      >  1:617055/1‑140 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                      >  1:91377/1‑140 (MQ=255)
                        tGAGAGTTTATTCTGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                    <  1:598611/140‑1 (MQ=255)
                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                    <  2:2048463/140‑1 (MQ=255)
                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                    <  1:108354/140‑1 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCag                                                                                                                   >  2:1371576/1‑140 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCag                                                                                                                   >  2:84159/1‑140 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                  <  2:872591/140‑1 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                  <  1:535273/140‑1 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                 <  2:760066/140‑1 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                 <  1:619079/140‑1 (MQ=255)
                             gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCagagag                                                                                                               >  1:1607074/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTTTAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAggg                                                                                                             >  2:2163852/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCAGAGAggg                                                                                                             >  2:1398891/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                             >  1:1233111/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                             <  2:1302217/140‑1 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                             <  1:2273825/140‑1 (MQ=255)
                                 aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCATTGGAAATGAGCGCAGAGAggggt                                                                                                           >  2:2024641/1‑138 (MQ=255)
                                   tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTTTGAAGTGAGCGCAGAGAgggggtt                                                                                                         >  2:390984/1‑136 (MQ=255)
                                     aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                       <  1:2024641/140‑1 (MQ=255)
                                       aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGAGAGGGCAGAGAGGGGTTATCtt                                                                                                     >  2:1170626/1‑140 (MQ=255)
                                       aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                     >  1:1837732/1‑140 (MQ=255)
                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                    <  2:679074/140‑1 (MQ=255)
                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                    <  2:2216106/140‑1 (MQ=255)
                                          gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCCGGGAGGGTTATCTTGCct                                                                                                 >  2:605744/1‑140 (MQ=255)
                                             tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCctga                                                                                               >  2:187504/1‑140 (MQ=255)
                                                   tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                                                                        <  2:266914/140‑1 (MQ=255)
                                                     aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                      >  1:1382005/1‑140 (MQ=255)
                                                         ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGa                                                                                  <  1:563363/140‑1 (MQ=255)
                                                           ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGGTATCTTGCCTGACTGAACTGGCGGtt                                                                                >  2:2174835/1‑140 (MQ=255)
                                                               gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGt                                                                            <  1:390984/140‑1 (MQ=255)
                                                                cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                           >  2:1608572/1‑140 (MQ=255)
                                                                cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                           <  2:1482059/140‑1 (MQ=255)
                                                                 tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGAGGGGAGGGAGGGGTTATTTTTCCTGGCTGAAATGGCGGTTTAGttt                                                                          >  2:299854/1‑140 (MQ=255)
                                                                     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGCGCAGAGAGGGGTTATCTTGGCTGACTGAACTGGCGATTAAGGTTTcgc                                                                      >  2:618613/1‑140 (MQ=255)
                                                                     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGGGCGCAGAGAGGGGTTATCTTTTCTGACTGAACTGGGGATTAAGttttccc                                                                      >  2:805008/1‑136 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCATTGGAAGGGGGCGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGATTAAGTTTTCCCCTAcc                                                                 >  2:1271523/1‑140 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTAcc                                                                 >  2:1096301/1‑140 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                 <  2:122508/140‑1 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                                <  1:266841/140‑1 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCATGTTGAGTGAGACGAGAGGGGGGTTAATTTGCCTTAATCTACCGGGGATTAAAGTTTCTCGTAcac                                                                >  2:1484419/1‑138 (MQ=255)
                                                                             ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGCCGGGAGGGAGTTTCTTTCCTGACTGAACTGGCGATTAAGGTTTCGCGGACCCgg                                                              >  2:666560/1‑140 (MQ=255)
                                                                               aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                            <  2:1657927/140‑1 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGGACCCGGCgg                                                           >  2:1150441/1‑139 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGCGCGGAGAGGGGGTTTTTTTCCTGACTGAACTGGGGATTAAGTTTTCGCGTACCCGGCgg                                                           >  2:509817/1‑139 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCgg                                                           >  2:777194/1‑139 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTGCATGAGAAGAGAGCGCAGGGGGGGGTTAATTTGCCTGCCTGAAATGGCGATTAAATTTTCCCGTACCCcgcgg                                                           >  2:2129174/1‑139 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGGGGTTAAGTTTCCGCGTACCCcgcgg                                                           >  2:381227/1‑139 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTCGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCggt                                                          >  2:1042637/1‑138 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGGGGGGGTTTGTTTCCTCACCTAAATTGGGGATAAGGTTTCCCGTAccccggggtt                                                         >  2:1551829/1‑133 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGAGCGCAGAGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGttt                                                         >  2:1796173/1‑140 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAATGGGGAGTTAGTTTTCGCGCACCCCGCGgtt                                                         >  2:1598116/1‑140 (MQ=255)
                                                                                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGGCGAGGGGGGTTATTTTGCCTGACTGAAATGGGGATTTAGTTTTCCCGTAACCGGCGTTTTGt                                                      >  2:220044/1‑140 (MQ=255)
                                                                                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                 <  1:2254442/140‑1 (MQ=255)
                                                                                           tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAAATTGCGGTTAAGTTTTCGCGGACCCGGCGGTTAGTTTTTGc                                                >  2:1827948/1‑140 (MQ=255)
                                                                                           tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTTAACTTGCGATTAAGTTTTCGCGGACCCCGCGTTTTGTTAtttt                                                >  2:2254442/1‑138 (MQ=255)
                                                                                            ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGCACCCGGCGTTTAGTTATTGCg                                               <  2:91377/140‑1 (MQ=255)
                                                                                              cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                             <  1:147556/140‑1 (MQ=255)
                                                                                               aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                            <  1:1486926/140‑1 (MQ=255)
                                                                                               aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                            <  2:680135/140‑1 (MQ=255)
                                                                                               aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                            <  1:74524/140‑1 (MQ=255)
                                                                                                 gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                          <  2:668077/140‑1 (MQ=255)
                                                                                                  aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                         <  2:2250271/140‑1 (MQ=255)
                                                                                                  aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                         <  1:1621962/140‑1 (MQ=255)
                                                                                                   tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                        <  2:48275/140‑1 (MQ=255)
                                                                                                    gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                       <  2:2072392/140‑1 (MQ=255)
                                                                                                    gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatga                                       <  1:2163852/140‑1 (MQ=255)
                                                                                                        tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGGGAAGGGAGCGCAGAGGGGGGTTTTCTTGCCTGACTGGACTGGCGATTAAGTTTTCGCGTACCCGGGGGTTTGTTATTTCGGGTGATGATAtt                                   >  2:67435/1‑140 (MQ=255)
                                                                                                        tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                   <  1:668755/140‑1 (MQ=255)
                                                                                                            gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                               <  1:1140858/140‑1 (MQ=255)
                                                                                                              ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                             <  1:2181386/140‑1 (MQ=255)
                                                                                                                    cGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGTGAAGTGAGAGCAGAGAGGGAATATCTTGCCTTAACGAACTGGGGGATAAGATTTCGCGTACCCCGCGTTTTGTTTTTGGCGATGGTGATATTGAAGcgccggtg                       >  2:357992/1‑133 (MQ=255)
                                                                                                                     gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                      <  2:1363452/140‑1 (MQ=255)
                                                                                                                      tttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAt                     <  1:362013/140‑1 (MQ=255)
                                                                                                                          cggcggTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAATTTTCGCGTACCCCGCGGTTAGTTTTTTCGGGTGATGATATTTAAGCTCGGCTGATTTGt                 >  2:1562445/1‑140 (MQ=255)
                                                                                                                            gcggTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc               >  2:1779818/1‑140 (MQ=255)
                                                                                                                             cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa              <  1:1263856/140‑1 (MQ=255)
                                                                                                                               gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTCCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt            <  2:488261/140‑1 (MQ=255)
                                                                                                                                   ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg        <  1:1154370/140‑1 (MQ=255)
                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt    <  1:1648530/140‑1 (MQ=255)
                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt    <  2:379945/140‑1 (MQ=255)
                                                                                                                                         ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCa  <  1:1796173/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCA  >  NC_000913/4604207‑4604483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: