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

CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  NC_000913/4604194‑4604494
                                                                                                                                                        |                                                                                                                                                      
cACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATtttttc                                                                                                                                                        <  1:1016399/151‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATtttttctttt                                                                                                                                                    <  1:626063/151‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATtttttcttt                                                                                                                                                     <  1:465581/150‑1 (MQ=255)
    cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATtttttctt                                                                                                                                                      <  1:173420/149‑1 (MQ=255)
      tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAATCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCAGCGG‑T‑GATtttttcttttg                                                                                                                                                   <  2:976728/150‑2 (MQ=255)
      tAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCACGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CAtttttttttttt                                                                                                                                                   >  1:315786/1‑150 (MQ=255)
       aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTttttctttttct                                                                                                                                                 <  1:475375/151‑1 (MQ=255)
        gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTg                                                                                                                                                <  1:902230/151‑1 (MQ=255)
               gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATga                                                                                                                                          <  2:727259/150‑1 (MQ=255)
                cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CAtttttttttttttt                                                                                                                                                 >  2:526396/1‑140 (MQ=255)
                cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATgaga                                                                                                                                        <  2:728350/151‑1 (MQ=255)
                   cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTCTTCCATGAGaa                                                                                                                                       >  1:639493/1‑149 (MQ=255)
                     tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CTTTTTTTTTTTTTTTGCCATGAGAAGTg                                                                                                                                    >  1:464807/1‑150 (MQ=255)
                     tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCACGTTTGCGGCGGTt‑tcttttttctttt                                                                                                                                                    >  2:1016343/1‑135 (MQ=255)
                       agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTAAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CCTTTTTTTTTTTTCTTCCagcagacgcggg                                                                                                                                  >  2:841142/1‑139 (MQ=255)
                          agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                               <  1:177986/150‑1 (MQ=255)
                          agtaTCGGTATTGAGAGTTTATTCAAAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCa                                                                                                                              <  2:200776/151‑1 (MQ=255)
                            taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTCTGCCATGAGAAGTGAgag                                                                                                                                >  2:771516/1‑145 (MQ=255)
                            taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑C‑CAATTTTTCTTTTTTTTCCTTGAGAAgatgg                                                                                                                                  >  2:990275/1‑141 (MQ=255)
                             aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCGTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                                                                                            <  2:672998/150‑1 (MQ=255)
                                  tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TGCATTTTTT‑TTTTTTTTCCATGGGAAGGGAGCGCAGAGAggg                                                                                                                       >  2:203958/1‑150 (MQ=255)
                                    ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGGAGAGAgg                                                                                                                        >  2:111846/1‑147 (MQ=255)
                                    ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGGGAGGGAtt                                                                                                                    >  1:680084/1‑151 (MQ=255)
                                      gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGAGGGAtt                                                                                                                    >  1:831595/1‑149 (MQ=255)
                                       agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                 <  2:69847/151‑1 (MQ=255)
                                           tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                               >  1:1010682/1‑149 (MQ=255)
                                            ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTTTTCCATGAGAAGGGAGGGCCGAGAGGGGTTATCtt                                                                                                               >  1:927466/1‑148 (MQ=255)
                                             tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑A‑CACTTTTTTTTTTTCTGCCATGAGAAGGGAGTGTAGAGAgggggtatttt                                                                                                               >  2:568456/1‑139 (MQ=255)
                                              aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGGGGGATTAtttt                                                                                                               >  2:680179/1‑146 (MQ=255)
                                              aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                          <  2:596468/151‑1 (MQ=255)
                                                tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCAG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                        >  1:808565/1‑151 (MQ=255)
                                                    aaaGTTTGCGGTTCACCCTTATAAGCTCCATACCTATACCCGTCAGGATTCATTTCAGGATGCCATGTCGCGCATCTCTTTTGCGGCGG‑T‑CATTTTTACTTATTCTGACATTAGAAGTAAGCGCAGAGAGGGCTTATCTTGCctgactga                                                                                                     <  2:179717/150‑1 (MQ=255)
                                                     aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CACTTTTTTTTTTTCTGCCATGAGAAGTGGGCGCAGAGAGGGATTTTCTTGCctgactg                                                                                                      >  1:805661/1‑148 (MQ=255)
                                                      aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                                                                                  <  2:1040383/151‑1 (MQ=255)
                                                            cGGGTAAACATGATAAGCGCCATACCTATACCAGTAAGGAGTTATTTCAGGATGCCATGTCGGGGTTCACGTTTACGGCGG‑T‑CATTTTTTCGTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCGTGCCTGACTGCGCTgg                                                                                                <  2:1044234/147‑1 (MQ=255)
                                                                        ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgcg                                                                                <  1:318534/151‑1 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑GATTTTTTTTTTTTTTGGCATGGGGAGTGAGCGCAGAGAGGGAGTATTTTGCTTGACTTAACGGGTGATTAAAtttttgcg                                                                                >  2:746612/1‑144 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGt                                                                               <  1:253131/150‑1 (MQ=255)
                                                                          aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTa                                                                              <  1:136551/151‑1 (MQ=255)
                                                                            gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                                            <  1:112081/151‑1 (MQ=255)
                                                                              tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTAccc                                                                           >  1:70243/1‑150 (MQ=255)
                                                                                 aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                       <  1:219646/151‑1 (MQ=255)
                                                                                  tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                         >  1:343161/1‑148 (MQ=255)
                                                                                  tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGGACCCcgcg                                                                       >  1:753903/1‑150 (MQ=255)
                                                                                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTCTTCCCTGAGAAGTGGGCGCCGAGGGGGGTTATCTTGCCTGACTGGACTGGCGATTTAGGTTTCGCGTCCCCGGCGttt                                                                    >  2:317121/1‑148 (MQ=255)
                                                                                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                  <  1:533854/150‑1 (MQ=255)
                                                                                         aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑AATTTTTTCTTTTTCTGCCATGAGGAGTGAGCGCTGAGAGGGAATTTCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTTGtt                                                                >  2:294412/1‑150 (MQ=255)
                                                                                          ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                              <  1:679373/151‑1 (MQ=255)
                                                                                          ccAGTCAGGAGTCATCTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTTGtt                                                                >  1:848815/1‑149 (MQ=255)
                                                                                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑G‑CATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGCTTATCTTGCCTGTCTGAACTGGCGATTAAGTTTTCGCGGACCCCGCGTTTTGTTAtt                                                             >  2:427833/1‑150 (MQ=255)
                                                                                            aGTCAGGAGGCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑G‑CATTTTTTTTTTTTCACTCATTAGAGGAGAGCGCAGAGGGGGGTTTTCTTTTCTGGCTGCAGAGGCGATTTAATTTACCCGTTACCGGAGttttggttttt                                                            >  2:409660/1‑143 (MQ=255)
                                                                                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCCTGGGAAGGGGGCGCAGAGGGGGGTTATCTTTGCTGACTTAACTGGCGATTTAGGTTCCGCGTACCCGGGGGTTTGGTAtt                                                             >  2:273768/1‑149 (MQ=255)
                                                                                                aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGTCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGGGTTTTGTTTTTGCgg                                                         >  2:149365/1‑150 (MQ=255)
                                                                                                   aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTTCCATGGGAAGTGAGCGCAGAGGGGGGTTTTCTTTCCTGACTGAACTGGCGATTTAGGTTCCGCGTACCCCGGGGTTTGTTTTTTCgggt                                                       >  2:737739/1‑147 (MQ=255)
                                                                                                    gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTTCCGCGTACCCGGCGTTTTGTTTTTGCgggt                                                       >  1:877020/1‑146 (MQ=255)
                                                                                                     tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGgatg                                                      >  1:208256/1‑148 (MQ=255)
                                                                                                       aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CCTTTTTTTTTTTTTTTCCCTGTGAAGTGAGGGCCGGGAGGGAGTTTCTTTCCTCACTTAACTTGCGATTTAGTTTCCGCGTACCCGGGGTTTAGTTTTTGTGggtg                                                      >  2:702587/1‑146 (MQ=255)
                                                                                                        tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTTGTTATTGCGggtgatg                                                   >  1:168408/1‑148 (MQ=255)
                                                                                                           cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                              <  2:300124/150‑1 (MQ=255)
                                                                                                            aGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGa                                            <  1:47943/151‑1 (MQ=255)
                                                                                                              gATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                          <  1:526885/151‑1 (MQ=255)
                                                                                                                   cATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                     <  1:3577/151‑1 (MQ=255)
                                                                                                                         gcgGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                               <  1:345855/151‑1 (MQ=255)
                                                                                                                           ggATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTgg                             <  1:745917/151‑1 (MQ=255)
                                                                                                                                tcGTTTGCGGCGG‑T‑CATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                        >  1:97294/1‑151 (MQ=255)
                                                                                                                                 cGTTTGCGGCGG‑T‑CATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGGGGGGTTTTCTTTTCTGACTTAAATTGCGATTAAGTTTTCGCGGACCCGGCGTTTTGTTATTTCGGGTGGTGATATTTAAAATCGGGTGAtttggttatt                       >  2:997181/1‑143 (MQ=255)
                                                                                                                                    ttGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                     >  1:666276/1‑150 (MQ=255)
                                                                                                                                        ggcgg‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc                 >  1:594383/1‑150 (MQ=255)
                                                                                                                                         gcgg‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt               >  1:67594/1‑151 (MQ=255)
                                                                                                                                           gg‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc              >  1:938117/1‑150 (MQ=255)
                                                                                                                                              t‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc           >  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)
                                                                                                                                                        |                                                                                                                                                      
CACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑T‑CATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGG  >  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: ‑