Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT40.3% 70.2 / 32.0 63S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (13/24);  new base T (25/0);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 7.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604194‑4604494
                                                                                                                                                        |                                                                                                                                                      
cACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttc                                                                                                                                                        <  1:1016399/151‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                    <  1:626063/151‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                                     <  1:465581/150‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                                      <  1:173420/149‑1 (MQ=255)
      tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAATCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCAGCGGTGATtttttcttttg                                                                                                                                                   <  2:976728/150‑2 (MQ=255)
      tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCACGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttt                                                                                                                                                   >  1:315786/1‑150 (MQ=255)
       aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                 <  1:475375/151‑1 (MQ=255)
        gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                                                                                <  1:902230/151‑1 (MQ=255)
               gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                                          <  2:727259/150‑1 (MQ=255)
                cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                                 >  2:526396/1‑140 (MQ=255)
                cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                        <  2:728350/151‑1 (MQ=255)
                   cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGaa                                                                                                                                       >  1:639493/1‑149 (MQ=255)
                     tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGTg                                                                                                                                    >  1:464807/1‑150 (MQ=255)
                     tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCACGTTTGCGGCGGTttcttttttctttt                                                                                                                                                    >  2:1016343/1‑135 (MQ=255)
                       agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTAAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCagcagacgcggg                                                                                                                                  >  2:841142/1‑139 (MQ=255)
                          agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                               <  1:177986/150‑1 (MQ=255)
                          agtaTCGGTATTGAGAGTTTATTCAAAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                              <  2:200776/151‑1 (MQ=255)
                            taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAgag                                                                                                                                >  2:771516/1‑145 (MQ=255)
                            taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCAATTTTTCTTTTTTTTCCTTGAGAAgatgg                                                                                                                                  >  2:990275/1‑141 (MQ=255)
                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCGTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                            <  2:672998/150‑1 (MQ=255)
                                  tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGCATTTTTTTTTTTTTTCCATGGGAAGGGAGCGCAGAGAggg                                                                                                                       >  2:203958/1‑150 (MQ=255)
                                    ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGGAGAGAgg                                                                                                                        >  2:111846/1‑147 (MQ=255)
                                    ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGGGAGGGAtt                                                                                                                    >  1:680084/1‑151 (MQ=255)
                                      gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGAGGGAtt                                                                                                                    >  1:831595/1‑149 (MQ=255)
                                       agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                 <  2:69847/151‑1 (MQ=255)
                                           tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                               >  1:1010682/1‑149 (MQ=255)
                                            ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGGGAGGGCCGAGAGGGGTTATCtt                                                                                                               >  1:927466/1‑148 (MQ=255)
                                             tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACACTTTTTTTTTTTCTGCCATGAGAAGGGAGTGTAGAGAgggggtatttt                                                                                                               >  2:568456/1‑139 (MQ=255)
                                              aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGGGGGATTAtttt                                                                                                               >  2:680179/1‑146 (MQ=255)
                                              aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                          <  2:596468/151‑1 (MQ=255)
                                                tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCAGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                        >  1:808565/1‑151 (MQ=255)
                                                    aaaGTTTGCGGTTCACCCTTATAAGCTCCATACCTATACCCGTCAGGATTCATTTCAGGATGCCATGTCGCGCATCTCTTTTGCGGCGGTCATTTTTACTTATTCTGACATTAGAAGTAAGCGCAGAGAGGGCTTATCTTGCctgactga                                                                                                     <  2:179717/150‑1 (MQ=255)
                                                     aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCACTTTTTTTTTTTCTGCCATGAGAAGTGGGCGCAGAGAGGGATTTTCTTGCctgactg                                                                                                      >  1:805661/1‑148 (MQ=255)
                                                      aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                                  <  2:1040383/151‑1 (MQ=255)
                                                            cGGGTAAACATGATAAGCGCCATACCTATACCAGTAAGGAGTTATTTCAGGATGCCATGTCGGGGTTCACGTTTACGGCGGTCATTTTTTCGTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCGTGCCTGACTGCGCTgg                                                                                                <  2:1044234/147‑1 (MQ=255)
                                                                        ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgcg                                                                                <  1:318534/151‑1 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTTTGGCATGGGGAGTGAGCGCAGAGAGGGAGTATTTTGCTTGACTTAACGGGTGATTAAAtttttgcg                                                                                >  2:746612/1‑144 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGt                                                                               <  1:253131/150‑1 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTa                                                                              <  1:136551/151‑1 (MQ=255)
                                                                            gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                            <  1:112081/151‑1 (MQ=255)
                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTAccc                                                                           >  1:70243/1‑150 (MQ=255)
                                                                                 aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                       <  1:219646/151‑1 (MQ=255)
                                                                                  tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                         >  1:343161/1‑148 (MQ=255)
                                                                                  tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCcgcg                                                                       >  1:753903/1‑150 (MQ=255)
                                                                                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGTGGGCGCCGAGGGGGGTTATCTTGCCTGACTGGACTGGCGATTTAGGTTTCGCGTCCCCGGCGttt                                                                    >  2:317121/1‑148 (MQ=255)
                                                                                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                  <  1:533854/150‑1 (MQ=255)
                                                                                         aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTCTTTTTCTGCCATGAGGAGTGAGCGCTGAGAGGGAATTTCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTTGtt                                                                >  2:294412/1‑150 (MQ=255)
                                                                                          ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                              <  1:679373/151‑1 (MQ=255)
                                                                                          ccAGTCAGGAGTCATCTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGtt                                                                >  1:848815/1‑149 (MQ=255)
                                                                                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGCTTATCTTGCCTGTCTGAACTGGCGATTAAGTTTTCGCGGACCCCGCGTTTTGTTAtt                                                             >  2:427833/1‑150 (MQ=255)
                                                                                            aGTCAGGAGGCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCACTCATTAGAGGAGAGCGCAGAGGGGGGTTTTCTTTTCTGGCTGCAGAGGCGATTTAATTTACCCGTTACCGGAGttttggttttt                                                            >  2:409660/1‑143 (MQ=255)
                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGGGGGCGCAGAGGGGGGTTATCTTTGCTGACTTAACTGGCGATTTAGGTTCCGCGTACCCGGGGGTTTGGTAtt                                                             >  2:273768/1‑149 (MQ=255)
                                                                                                aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGTCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGGGTTTTGTTTTTGCgg                                                         >  2:149365/1‑150 (MQ=255)
                                                                                                   aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGCGCAGAGGGGGGTTTTCTTTCCTGACTGAACTGGCGATTTAGGTTCCGCGTACCCCGGGGTTTGTTTTTTCgggt                                                       >  2:737739/1‑147 (MQ=255)
                                                                                                    gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTTCCGCGTACCCGGCGTTTTGTTTTTGCgggt                                                       >  1:877020/1‑146 (MQ=255)
                                                                                                     tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGgatg                                                      >  1:208256/1‑148 (MQ=255)
                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGTGAAGTGAGGGCCGGGAGGGAGTTTCTTTCCTCACTTAACTTGCGATTTAGTTTCCGCGTACCCGGGGTTTAGTTTTTGTGggtg                                                      >  2:702587/1‑146 (MQ=255)
                                                                                                        tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGggtgatg                                                   >  1:168408/1‑148 (MQ=255)
                                                                                                           cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                              <  2:300124/150‑1 (MQ=255)
                                                                                                            aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                            <  1:47943/151‑1 (MQ=255)
                                                                                                              gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                          <  1:526885/151‑1 (MQ=255)
                                                                                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                     <  1:3577/151‑1 (MQ=255)
                                                                                                                         gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                               <  1:345855/151‑1 (MQ=255)
                                                                                                                           ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  1:745917/151‑1 (MQ=255)
                                                                                                                                tcGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                        >  1:97294/1‑151 (MQ=255)
                                                                                                                                 cGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTCTTTTCTGACTTAAATTGCGATTAAGTTTTCGCGGACCCGGCGTTTTGTTATTTCGGGTGGTGATATTTAAAATCGGGTGAtttggttatt                       >  2:997181/1‑143 (MQ=255)
                                                                                                                                    ttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                     >  1:666276/1‑150 (MQ=255)
                                                                                                                                        ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                 >  1:594383/1‑150 (MQ=255)
                                                                                                                                         gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt               >  1:67594/1‑151 (MQ=255)
                                                                                                                                           ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              >  1:938117/1‑150 (MQ=255)
                                                                                                                                              tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc           >  2:272128/1‑151 (MQ=255)
                                                                                                                                                 aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc         >  2:998951/1‑151 (MQ=255)
                                                                                                                                                      ttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc    >  2:267713/1‑151 (MQ=255)
                                                                                                                                                        ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg  >  1:672834/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604194‑4604494

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