Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT40.0% 69.4 / 29.8 62S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (9/27);  new base T (24/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604200‑4604490
                                                                                                                                                 |                                                                                                                                                   
tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttt                                                                                                                                                >  1:416411/1‑150 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                              >  2:327794/1‑149 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                                                                               <  1:246252/150‑1 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                              <  1:1149593/151‑1 (MQ=255)
    aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                              >  2:596086/1‑146 (MQ=255)
    aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                                                                           <  2:924915/151‑1 (MQ=255)
     aaaggcggtaTGAGTATTAACGGTATTGAGATTTTATTAAGAAAGTTTGCAGGTAACCTTAATAGGCTCAATGACTAAACGGGTTTGGTTTCTTTTCACGATACCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTGTTTCTGcgg                                                                                                                                          <  2:1097496/142‑3 (MQ=255)
     aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTtttttttttttttt                                                                                                                                              >  2:440400/1‑145 (MQ=255)
      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt                                                                                                                                         <  2:704204/151‑1 (MQ=255)
      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt                                                                                                                                         <  2:1019184/151‑1 (MQ=255)
       aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTCTTTTTTTGCCAt                                                                                                                                         >  1:190193/1‑150 (MQ=255)
              aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                  <  2:1140220/150‑1 (MQ=255)
                gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                               <  2:1014303/151‑1 (MQ=255)
                   tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAAttttttttttttttcc                                                                                                                                           >  2:982321/1‑131 (MQ=255)
                    agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTATTTTTTTTCCATGAGAAGTGAgcg                                                                                                                             >  1:1173802/1‑149 (MQ=255)
                      taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                          <  2:786096/150‑1 (MQ=255)
                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                        <  1:437737/151‑1 (MQ=255)
                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                         <  2:924801/150‑1 (MQ=255)
                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCGTGGAAAAGGAGCGCggag                                                                                                                        >  1:1134287/1‑149 (MQ=255)
                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCCTGAGAAGTGAGCGCag                                                                                                                          >  2:58292/1‑147 (MQ=255)
                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTTTTCCATGAGGAGGGTGCggagggag                                                                                                                      >  2:23823/1‑150 (MQ=255)
                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGAGCGCAGGGAgg                                                                                                                     >  1:1190526/1‑151 (MQ=255)
                                gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                               <  1:583177/151‑1 (MQ=255)
                                       tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTTTCATTAGGAGTGGAGGCGGAGAgggggtatttt                                                                                                            >  2:513173/1‑139 (MQ=255)
                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                       <  2:806826/150‑1 (MQ=255)
                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct                                                                                                       <  1:875695/150‑1 (MQ=255)
                                            aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTATTTTTAATTTTTTGCCTTCAGAGGGAGCGCGGAGGGGGATTGTCTTGGctgcacg                                                                                                  >  1:1277435/1‑147 (MQ=255)
                                             gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGATTATCttttcttact                                                                                                   >  1:397854/1‑141 (MQ=255)
                                              aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                 <  2:1257967/150‑1 (MQ=255)
                                               aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATCAGAAGAGAGCGGCAAGAGGGGCTATCTTTCcttactg                                                                                                  >  2:194603/1‑148 (MQ=255)
                                                aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCAGGGAGGGGTTATCTTTCCTGACTGGACt                                                                                              >  2:1022422/1‑151 (MQ=255)
                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTATTTTGCctgact                                                                                                   >  1:213572/1‑145 (MQ=255)
                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGaa                                                                                                >  1:770292/1‑148 (MQ=255)
                                                       gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAtt                                                                                       <  1:709921/151‑1 (MQ=255)
                                                         gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCCTTGGAAGGGAGGGGCCGGGGGGATTATCATTCCCTGCTGAACTGGCgggt                                                                                       >  2:677557/1‑146 (MQ=255)
                                                             cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                 <  2:941064/151‑1 (MQ=255)
                                                              ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCGTTTTTTTTTTTCTGGCCAGAGAAGGGAGGGCACAGGGGGATTATGTTGTCTTTATTTACTGGCGATTAAgct                                                                                  >  2:2227/1‑147 (MQ=255)
                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                               <  1:233742/149‑1 (MQ=255)
                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCGGTTAAGttttc                                                                               >  1:376714/1‑147 (MQ=255)
                                                                   taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCCGAGAGGGGTTTTTTTTCCTGACTGAACTGGCGATTAAGTTTCcg                                                                              >  1:639572/1‑148 (MQ=255)
                                                                      gCTCCATACCTGTACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGTTATTTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGGAcc                                                                        >  1:26184/1‑151 (MQ=255)
                                                                      gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                        <  1:829233/151‑1 (MQ=255)
                                                                          cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGGGGGGTTATTTTTCCTGACTGAACTGGGGGTTAAGTTTTcgcggaacc                                                                       >  2:930972/1‑143 (MQ=255)
                                                                              ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                <  1:268028/151‑1 (MQ=255)
                                                                                   aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                           <  2:822629/151‑1 (MQ=255)
                                                                                    ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                          <  1:488651/151‑1 (MQ=255)
                                                                                      aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                        <  1:639250/151‑1 (MQ=255)
                                                                                       gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                       <  1:408404/151‑1 (MQ=255)
                                                                                        tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGCGGAGTGAGCGCCGAGAGGGGTTATTTTTCCTGACTGAACCTGCGGTTAAGTTTCCGCCTACCCCGCGTTTTGTTTTTGCg                                                      >  1:663196/1‑151 (MQ=255)
                                                                                            gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTATTGCggg                                                    >  1:620451/1‑148 (MQ=255)
                                                                                              gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGTGAAGTGAGGGCGGAGGGGGGTTATCTTGCCTGACTTAACTGGCGAGTAAGTTTTCGCGTTCCCCGCGCTTTGTTTTTGCgggt                                                   >  1:55531/1‑146 (MQ=255)
                                                                                                 aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCCGCGGTTAGTTATTTCGggtgatg                                               >  1:683905/1‑149 (MQ=255)
                                                                                                    tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCAGAGAGGGATTATCTTGTCTGACTGAACTGGGGATTAAGTTTTCCCGTACCCGGGGCTTTGTTATTTCGGGTGATGATAtt                                          >  1:1159192/1‑151 (MQ=255)
                                                                                                    acAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTTTTTTTTTTGCCCTTAGGAGGGGGGGCCTGGAGGGATTTTCTTTGCTGACCTAAATGGCGGTTAAGTTTTCGCGTCACCGGCGGTTTGTTATTTCGGATGGTGAAAtt                                          >  1:865573/2‑151 (MQ=255)
                                                                                                             cATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGGGAAGGGTGTGCTGGGAGGGATTTTCTTGCCTGACTGAATTGGCGGTTAAGTTTCCGCGTACCAGGGGTTTCGGTATTTTGGATGTTGATTTTGTAGTTc                                   >  1:135433/1‑149 (MQ=255)
                                                                                                              aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                <  1:338614/151‑1 (MQ=255)
                                                                                                                gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                              <  1:567635/151‑1 (MQ=255)
                                                                                                                   gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTTTTTCGGGTGATGATATTTAAGCTCGGCTggt                            >  2:923749/1‑148 (MQ=255)
                                                                                                                         ctcGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCCGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                       >  1:1107106/1‑149 (MQ=255)
                                                                                                                           cGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGGAGAGAGCGGGCGGAGGGGTTATCTTTGCTTAATGTACTGGCGGTTTAGTTTCCGCGGAACCCGCGTCTTGTTATTGCGGGTGATGTTATTGAAGATCgcgcggttggtgtctt                   >  2:437912/1‑136 (MQ=255)
                                                                                                                           cGTTTGCGGCCGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCGGAGAGGGATTATCTTGCCTGAGTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTTATTGGTTCa                     <  2:1151825/149‑1 (MQ=255)
                                                                                                                            gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                  <  1:244919/151‑1 (MQ=255)
                                                                                                                               tGCGGCGGTCATTTTTTTTTTTTTTTCCCTGCGAAGTGAGCGCAGAGGGGGGTTTTCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTGGTTGGGTTATTGt                 >  2:804299/1‑149 (MQ=255)
                                                                                                                                gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg               <  2:971657/150‑1 (MQ=255)
                                                                                                                                     ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc          >  1:296951/1‑150 (MQ=255)
                                                                                                                                        taTTTTTTCTTTTTCTGCGATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCTTGTCTGAACTGGCGGTTAAGTTTCCGCGTCCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      >  2:292691/2‑151 (MQ=255)
                                                                                                                                        cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      >  1:607343/1‑151 (MQ=255)
                                                                                                                                        cATTTTTTCTTTTTCTGCCATGAGAAGAGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc       >  2:755754/1‑150 (MQ=255)
                                                                                                                                          ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCTTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCTTTGTCGCCATCAc        >  2:493020/1‑147 (MQ=255)
                                                                                                                                               ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCGTCACCGCTgg  >  2:528737/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604200‑4604490

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