Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,244,248 G→C 5.7% R78R (CGC→CGG hmgA ← homogentisate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,244,2480GC 5.7% 230.1 / 3.6 70R78R (CGC→CGGhmgAhomogentisate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (27/39);  new base C (2/2);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACA  >  NC_002947/5244103‑5244395
                                                                                                                                                 |                                                                                                                                                   
ccGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAATCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGccgc                                                                                                                                                  >  2:192349/1‑149 (MQ=255)
  ggTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAATCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGCCTGGCGCGCCAGCCGCTc                                                                                                                                                <  1:192349/149‑1 (MQ=255)
          cGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTAGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGt                                                                                                                                         <  3:51236/148‑1 (MQ=255)
           gCTCCGTCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTgg                                                                                                                                       <  4:31378/149‑1 (MQ=255)
           gCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTgg                                                                                                                                       <  5:299825/149‑1 (MQ=255)
             tCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGt                                                                                                                                     >  7:211692/1‑149 (MQ=255)
               gaaggTTGGCTGCCCATGGCTGCCATCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGgcgc                                                                                                                                                  <  8:137249/130‑1 (MQ=255)
               gaaggTTGGCTGCCCATGGCTGCCATCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGgcgc                                                                                                                                                  >  7:137249/5‑134 (MQ=255)
               cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                                                                                                                                             <  6:123167/139‑1 (MQ=255)
               cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                                                                                                                                             >  5:123167/1‑139 (MQ=255)
               cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTg                                                                                                                                    <  6:180287/148‑1 (MQ=255)
                 gCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCaa                                                                                                                                 <  8:17300/149‑1 (MQ=255)
                 gCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCaa                                                                                                                                 <  3:154569/149‑1 (MQ=255)
                   gTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                               >  7:254683/1‑149 (MQ=255)
                      ggCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                               >  3:96216/1‑146 (MQ=255)
                        cAGCCATGGGCAGGCAACCTTCGATGAAATCGGTGGGTTCAGCAGGAATCGGGTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGGGACCCTGTGGCCGGCGGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGTAGaa                                                                                                                          <  8:82515/149‑1 (MQ=255)
                           ccATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                               <  6:318073/141‑1 (MQ=255)
                           ccATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                               >  5:318073/1‑141 (MQ=255)
                            cATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGc                                                                                                                      <  1:49739/149‑1 (MQ=255)
                                ggCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCTGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGaa                                                                                                                   <  8:254683/148‑1 (MQ=255)
                                ggCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAAt                                                                                                                  <  5:263718/149‑1 (MQ=255)
                                   aGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCg                                                                                                               <  8:227374/149‑1 (MQ=255)
                                       aaCCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTa                                                                                                            >  6:332000/1‑148 (MQ=255)
                                          cTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGAGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGc                                                                                                        >  6:192510/1‑149 (MQ=255)
                                               aTGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAgg                                                                                                    >  6:82148/1‑148 (MQ=255)
                                                  aaaTCGGTCGGTTCAGCAGGAATCGGCCGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACg                                                                                                >  6:14849/1‑149 (MQ=255)
                                                     tcggtcggTTCCGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGgcg                                                                                             <  5:88879/149‑1 (MQ=255)
                                                            gtTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCa                                                                                           <  2:331650/144‑1 (MQ=255)
                                                            gtTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCa                                                                                           >  1:331650/1‑144 (MQ=255)
                                                               cagcagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCt                                                                                    <  5:332000/148‑1 (MQ=255)
                                                                cgcagGCTTCGGCTGCGGGCTCCCGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGc                                                                                  <  8:211692/148‑1 (MQ=255)
                                                                 gcagGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg                                                                                                                               >  2:95339/1‑103 (MQ=255)
                                                                      aaTCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCa                                                                            <  5:341639/149‑1 (MQ=255)
                                                                        tCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATg                                                                          <  6:134020/149‑1 (MQ=255)
                                                                          ggCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATgg                                                                         <  4:96216/148‑1 (MQ=255)
                                                                          ggCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                        <  2:212315/149‑1 (MQ=255)
                                                                           gCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGt                                                                        >  2:202850/1‑148 (MQ=255)
                                                                                gggCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGc                                                                  >  2:14644/1‑149 (MQ=255)
                                                                                   cTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                               <  5:174659/149‑1 (MQ=255)
                                                                                    tCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTg                                                              >  5:59918/1‑149 (MQ=255)
                                                                                      cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATg                                                                                                                 <  1:201605/96‑1 (MQ=255)
                                                                                      cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATg                                                                                                                 >  2:201605/1‑96 (MQ=255)
                                                                                      cAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGcc                                                            >  8:69642/1‑149 (MQ=255)
                                                                                       aGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGccc                                                           <  7:208457/149‑1 (MQ=255)
                                                                                          gcAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAc                                                        <  6:59918/149‑1 (MQ=255)
                                                                                            aGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAcag                                                      >  3:328859/1‑149 (MQ=255)
                                                                                                ggTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACagcagc                                                  <  5:14849/149‑1 (MQ=255)
                                                                                                 gTTGGGGTTGATGCCAGCCAGTGGCCGGCCGTGCGGCTGGGGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGGAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCTTGGTGAACGCCGTGCCCGACAGCAGCt                                                 <  4:15392/149‑1 (MQ=255)
                                                                                                 gTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCt                                                 >  5:336690/1‑149 (MQ=255)
                                                                                                  ttGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGTGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTc                                                >  4:161508/1‑149 (MQ=255)
                                                                                                     gggTGGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGc                                             <  5:342808/149‑1 (MQ=255)
                                                                                                      ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                            >  3:178696/1‑149 (MQ=255)
                                                                                                      ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                            >  3:48925/1‑149 (MQ=255)
                                                                                                      ggTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCa                                            >  4:75721/1‑149 (MQ=255)
                                                                                                           aTGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAg                                       >  3:263393/1‑149 (MQ=255)
                                                                                                               cACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCa                                   >  1:214320/1‑149 (MQ=255)
                                                                                                                    aGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACgg                              <  5:157710/149‑1 (MQ=255)
                                                                                                                        gCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGcc                          >  4:201312/1‑149 (MQ=255)
                                                                                                                            ttcgAGCGCCTGGCGCGCCAGCCGCTCGGAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCt                      >  2:22386/3‑149 (MQ=255)
                                                                                                                              ggAGCGGCTGGCGCGCCAGGCGCTCGTAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCGCGCTGCGGGCCATTGTGAACGCCGTGCCAGACATCAGCTCGGCATACAGGCCATACGGGGCCTTCTGc                    <  5:192510/148‑1 (MQ=255)
                                                                                                                                  cGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAcagc                                                     >  7:134974/1‑112 (MQ=255)
                                                                                                                                  cGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGAcagc                                                     <  8:134974/112‑1 (MQ=255)
                                                                                                                                   ggCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGa               <  1:14644/149‑1 (MQ=255)
                                                                                                                                   ggCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGa               <  4:154260/149‑1 (MQ=255)
                                                                                                                                        gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGTGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGtt            <  3:161508/147‑1 (MQ=255)
                                                                                                                                        gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                               <  7:189025/96‑1 (MQ=255)
                                                                                                                                        gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGt                                                               >  8:189025/1‑96 (MQ=255)
                                                                                                                                        gcgcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCt          <  4:263393/149‑1 (MQ=255)
                                                                                                                                          gcgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGc        <  5:82148/149‑1 (MQ=255)
                                                                                                                                          gcgcCAGCCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCgg                  >  5:63646/1‑139 (MQ=255)
                                                                                                                                          gcgcCAGCCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCgg                  <  6:63646/139‑1 (MQ=255)
                                                                                                                                           cgcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGcc       <  4:328859/149‑1 (MQ=255)
                                                                                                                                            gcCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGccc      <  4:178696/149‑1 (MQ=255)
                                                                                                                                                ggCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACa  <  4:48925/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGCGAATGCGGTACAGCCAGGTACGGCGCAGCTCGCTGCGGGCCATGGTGAACGCCGTGCCCGACAGCAGCTCGGCATACAGGCCATACGGGGCCTTCTGCGGGGAGTTCTGCCCAACA  >  NC_002947/5244103‑5244395

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