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

TAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604200‑4604490
                                                                                                                                                 |                                                                                                                                                   
tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAt‑ttttttttttt                                                                                                                                                >  1:416411/1‑150 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAt‑ttttttttttttt                                                                                                                                              >  2:327794/1‑149 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑tttttctttttc                                                                                                                                               <  1:246252/150‑1 (MQ=255)
 aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑Tttttctttttct                                                                                                                                              <  1:1149593/151‑1 (MQ=255)
    aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAt‑ttttttttttttt                                                                                                                                              >  2:596086/1‑146 (MQ=255)
    aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGcc                                                                                                                                           <  2:924915/151‑1 (MQ=255)
     aaaggcggtaTGAGTATTAACGGTATTGAGATTTTATTAAGAAAGTTTGCAGGTAACCTTAATAGGCTCAATGACTAAACGGGTTTGGTTTCTTTTCACGATACCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTGTTTCTGcgg                                                                                                                                          <  2:1097496/142‑3 (MQ=255)
     aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTt‑ttttttttttttt                                                                                                                                              >  2:440400/1‑145 (MQ=255)
      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCAt                                                                                                                                         <  2:704204/151‑1 (MQ=255)
      aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCAt                                                                                                                                         <  2:1019184/151‑1 (MQ=255)
       aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCT‑TTTTTCTTTTTTTGCCAt                                                                                                                                         >  1:190193/1‑150 (MQ=255)
              aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                  <  2:1140220/150‑1 (MQ=255)
                gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                               <  2:1014303/151‑1 (MQ=255)
                   tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAA‑ttttttttttttttcc                                                                                                                                           >  2:982321/1‑131 (MQ=255)
                    agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTATTTTTTTTCCATGAGAAGTGAgcg                                                                                                                             >  1:1173802/1‑149 (MQ=255)
                      taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                          <  2:786096/150‑1 (MQ=255)
                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagag                                                                                                                        <  1:437737/151‑1 (MQ=255)
                       aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                         <  2:924801/150‑1 (MQ=255)
                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTCTTTTTTTTCCGTGGAAAAGGAGCGCggag                                                                                                                        >  1:1134287/1‑149 (MQ=255)
                         cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTTCCTGAGAAGTGAGCGCag                                                                                                                          >  2:58292/1‑147 (MQ=255)
                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGAT‑TTTTTTTTTTTTTTCCATGAGGAGGGTGCggagggag                                                                                                                      >  2:23823/1‑150 (MQ=255)
                          ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTTTTCCATGAGAAGTGAGCGCAGGGAgg                                                                                                                     >  1:1190526/1‑151 (MQ=255)
                                gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                               <  1:583177/151‑1 (MQ=255)
                                       tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTTTCATTAGGAGTGGAGGCGGAGAgggggtatttt                                                                                                            >  2:513173/1‑139 (MQ=255)
                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCct                                                                                                       <  2:806826/150‑1 (MQ=255)
                                        aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCct                                                                                                       <  1:875695/150‑1 (MQ=255)
                                            aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTATTTTTAATTTTTTGCCTTCAG‑AGGGAGCGCGGAGGGGGATTGTC‑TTGGctgcacg                                                                                                  >  1:1277435/1‑147 (MQ=255)
                                             gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGATTATC‑ttttcttact                                                                                                   >  1:397854/1‑141 (MQ=255)
                                              aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCctgactga                                                                                                 <  2:1257967/150‑1 (MQ=255)
                                               aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATCAGAAGAGAGCGGCAAGAGGGGCTATC‑TTTCcttactg                                                                                                  >  2:194603/1‑148 (MQ=255)
                                                aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCAGGGAGGGGTTATC‑TTTCCTGACTGGACt                                                                                              >  2:1022422/1‑151 (MQ=255)
                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTATT‑TTGCctgact                                                                                                   >  1:213572/1‑145 (MQ=255)
                                                 gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATT‑TTGCCTGACTGaa                                                                                                >  1:770292/1‑148 (MQ=255)
                                                       gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGAtt                                                                                       <  1:709921/151‑1 (MQ=255)
                                                         gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCAT‑TTTTTTTTTTTTTGCCCTTGGAAGGGAGGGGCCGGGGGGATTATCATTCCCTG‑CTGAACTGGCgggt                                                                                       >  2:677557/1‑146 (MQ=255)
                                                             cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGttt                                                                                 <  2:941064/151‑1 (MQ=255)
                                                              ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCG‑TTTTTTTTTTTCTGGCCAGAGAAGGGAGGGCACAGGGGGATTATG‑TTGTCTTTATTTACTGGCGATTAAgct                                                                                  >  2:2227/1‑147 (MQ=255)
                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                               <  1:233742/149‑1 (MQ=255)
                                                                 tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAAATGGCGGTTAAGttttc                                                                               >  1:376714/1‑147 (MQ=255)
                                                                   taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCCGAGAGGGGTTTTT‑TTTCCTGACTGAACTGGCGATTAAGTTTCcg                                                                              >  1:639572/1‑148 (MQ=255)
                                                                      gCTCCATACCTGTACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGTTATT‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGGAcc                                                                        >  1:26184/1‑151 (MQ=255)
                                                                      gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                        <  1:829233/151‑1 (MQ=255)
                                                                          cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGGGGGGTTATT‑TTTCCTGACTGAACTGGGGGTTAAGTTTTcgcggaacc                                                                       >  2:930972/1‑143 (MQ=255)
                                                                              ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt                                                                <  1:268028/151‑1 (MQ=255)
                                                                                   aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                           <  2:822629/151‑1 (MQ=255)
                                                                                    ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                          <  1:488651/151‑1 (MQ=255)
                                                                                      aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                        <  1:639250/151‑1 (MQ=255)
                                                                                       gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                       <  1:408404/151‑1 (MQ=255)
                                                                                        tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTTCCATGCGGAGTGAGCGCCGAGAGGGGTTATT‑TTTCCTGACTGAACCTGCGGTTAAGTTTCCGCCTACCCCGCGTTTTGTTTTTGCg                                                      >  1:663196/1‑151 (MQ=255)
                                                                                            gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAAATGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTATTGCggg                                                    >  1:620451/1‑148 (MQ=255)
                                                                                              gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGTGAAGTGAGGGCGGAGGGGGGTTATC‑TTGCCTGACTTAACTGGCGAGTAAGTTTTCGCGTTCCCCGCGCTTTGTTTTTGCgggt                                                   >  1:55531/1‑146 (MQ=255)
                                                                                                 aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCCGCGGTTAGTTATTTCGggtgatg                                               >  1:683905/1‑149 (MQ=255)
                                                                                                    tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCAGAGAGGGATTATC‑TTGTCTGACTGAACTGGGGATTAAGTTTTCCCGTACCCGGGGCTTTGTTATTTCGGGTGATGATAtt                                          >  1:1159192/1‑151 (MQ=255)
                                                                                                    acAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTAT‑TTTTTTTTTTTTTGCCCTTAGGAGGGGGGGCCTGGAGGGATTTTC‑TTTGCTGACCTAAATGGCGGTTAAGTTTTCGCGTCACCGGCGGTTTGTTATTTCGGATGGTGAAAtt                                          >  1:865573/2‑151 (MQ=255)
                                                                                                             cATGTCGCGGATCTCGTTTGCGGCGGTCCT‑TTTTTTTTTTTCTGCCATGGGAAGGGTGTGCTGGGAGGGATTTTC‑TTGCCTGACTGAATTGGCGGTTAAGTTTCCGCGTACCAGGGGTTTCGGTATTTTGGATGTTGATTTTGTAGTTc                                   >  1:135433/1‑149 (MQ=255)
                                                                                                              aTGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                <  1:338614/151‑1 (MQ=255)
                                                                                                                gTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                              <  1:567635/151‑1 (MQ=255)
                                                                                                                   gcgGATCTCGTTTGCGGCGGTCAT‑TTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAGGGGTTTTC‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGTTTTGTTTTTTCGGGTGATGATATTTAAGCTCGGCTggt                            >  2:923749/1‑148 (MQ=255)
                                                                                                                         ctcGTTTGCGGCGGTCAT‑TTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCCGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                       >  1:1107106/1‑149 (MQ=255)
                                                                                                                           cGTTTGCGGCGGTCAT‑TTTTTCTTTTTTTGCCATGAGGAGAGAGCGGGCGGAGGGGTTATC‑TTTGCTTAATGTACTGGCGGTTTAGTTTCCGCGGAACCCGCGTCTTGTTATTGCGGGTGATGTTATTGAAGATCgcgcggttggtgtctt                   >  2:437912/1‑136 (MQ=255)
                                                                                                                           cGTTTGCGGCCGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCGGAGAGGGATTATC‑TTGCCTGAGTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTTATTGGTTCa                     <  2:1151825/149‑1 (MQ=255)
                                                                                                                            gTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg                  <  1:244919/151‑1 (MQ=255)
                                                                                                                               tGCGGCGGTCAT‑TTTTTTTTTTTTTTCCCTGCGAAGTGAGCGCAGAGGGGGGTTTTC‑TTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTGCGGGTGATGATATTTAAGCTCGGCTGGTTGGGTTATTGt                 >  2:804299/1‑149 (MQ=255)
                                                                                                                                gcggcgGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg               <  2:971657/150‑1 (MQ=255)
                                                                                                                                     ggTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc          >  1:296951/1‑150 (MQ=255)
                                                                                                                                        taT‑TTTTTCTTTTTCTGCGATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCTTGTCTGAACTGGCGGTTAAGTTTCCGCGTCCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      >  2:292691/2‑151 (MQ=255)
                                                                                                                                        cAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg      >  1:607343/1‑151 (MQ=255)
                                                                                                                                        cAT‑TTTTTCTTTTTCTGCCATGAGAAGAGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc       >  2:755754/1‑150 (MQ=255)
                                                                                                                                          t‑tttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCTTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCTTTGTCGCCATCAc        >  2:493020/1‑147 (MQ=255)
                                                                                                                                               ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCGTCACCGCTgg  >  2:528737/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAT‑TTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATC‑TTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  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: ‑