Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,731,814 Δ1 bp 5.5% coding (157/2574 nt) opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,731,8140C. 5.5% 862.7 / 42.3 256coding (157/2574 nt)opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base C (105/136);  new base . (11/3);  total (117/139)
Fisher's exact test for biased strand distribution p-value = 1.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 260 total aligned reads displayed.

CGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAGT  >  NC_002947/5731666‑5731963
                                                                                                                                                     |                                                                                                                                                      
   cATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGcc                                                                                                                                                    >  2:3536358/1‑150 (MQ=255)
    aTGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCt                                                                                                                                                   >  2:3899901/1‑150 (MQ=255)
     tGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGTGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                                  <  2:3043451/150‑1 (MQ=255)
     tGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                                  <  2:4660395/150‑1 (MQ=255)
     tGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                                  <  2:1375279/150‑1 (MQ=255)
       gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                                <  2:4437720/150‑1 (MQ=255)
       gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                                <  2:1143894/150‑1 (MQ=255)
       gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                                <  1:5011879/150‑1 (MQ=255)
       gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                                <  1:2788792/150‑1 (MQ=255)
        cTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCGTGCAGGCCTCggt                                                                                                                                               <  1:5241495/150‑2 (MQ=255)
         tCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCa                                                                                                                                              >  2:2071904/1‑150 (MQ=255)
          cGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCggtgg                                                                                                                                             <  1:1057926/150‑4 (MQ=255)
            ggCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCagag                                                                                                                                           <  2:4529685/150‑1 (MQ=255)
             gCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGc                                                                                                                                          >  2:50029/1‑150 (MQ=255)
               aCGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGCGGCCTGCACGGCCTCGGCAGAGCTg                                                                                                                                        >  1:338106/1‑150 (MQ=255)
               aCGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTg                                                                                                                                        >  2:4043617/1‑150 (MQ=255)
               aCGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTg                                                                                                                                        >  1:1043493/1‑150 (MQ=255)
               aCGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTg                                                                                                                                        <  1:6373599/150‑1 (MQ=255)
                 gACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTggc                                                                                                                                      >  2:4305069/1‑150 (MQ=255)
                 gACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTggc                                                                                                                                      <  1:5382769/150‑1 (MQ=255)
                   ccTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTggcgg                                                                                                                                    >  2:2850235/1‑150 (MQ=255)
                    cTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTggcggc                                                                                                                                   <  1:2186335/150‑1 (MQ=255)
                    cTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCGCGGCAGAGCAggcggc                                                                                                                                   >  1:3499605/1‑150 (MQ=255)
                      tGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCg                                                                                                                                 <  1:6054428/150‑1 (MQ=255)
                       ggTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCgg                                                                                                                                >  2:159563/1‑150 (MQ=255)
                          gCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGtt                                                                                                                             >  2:659212/1‑150 (MQ=255)
                          gCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGtt                                                                                                                             <  2:5825977/150‑1 (MQ=255)
                          gCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGtt                                                                                                                             <  1:221962/150‑1 (MQ=255)
                            gTTTGCTCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTgg                                                                                                                           <  2:5114045/150‑1 (MQ=255)
                             tttGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGgc                                                                                                                          >  1:4080825/1‑150 (MQ=255)
                                        gggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGtt                                                                                                               <  1:4734917/150‑1 (MQ=255)
                                        gggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCGTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGtt                                                                                                               <  2:5514486/150‑1 (MQ=255)
                                         ggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTg                                                                                                              <  1:3931913/150‑1 (MQ=255)
                                         ggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTg                                                                                                              <  1:4886752/150‑1 (MQ=255)
                                           cGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgt                                                                                                            <  1:1569757/150‑1 (MQ=255)
                                           cGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgt                                                                                                            <  1:674601/150‑1 (MQ=255)
                                           cGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgt                                                                                                            <  1:3142123/150‑1 (MQ=255)
                                            gATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgtg                                                                                                           <  2:2158819/150‑1 (MQ=255)
                                              tCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTa                                                                                                         <  1:1936453/150‑1 (MQ=255)
                                                ttCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGc                                                                                                       >  2:785457/1‑150 (MQ=255)
                                                ttCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGc                                                                                                       >  2:6068466/1‑150 (MQ=255)
                                                  cAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTTGCCGCCAAGCGTTGGTGTAGCTc                                                                                                     >  1:3350747/1‑150 (MQ=255)
                                                   aGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCa                                                                                                    >  1:3868770/1‑150 (MQ=255)
                                                      cacaGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTg                                                                                                 >  2:3590836/1‑150 (MQ=255)
                                                            cggcCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAg                                                                                           <  1:5716004/150‑1 (MQ=255)
                                                            cggcCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAg                                                                                           <  2:2999517/150‑1 (MQ=255)
                                                             ggcCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAgg                                                                                          <  2:311912/150‑1 (MQ=255)
                                                             ggcCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAgg                                                                                          <  1:2033233/150‑1 (MQ=255)
                                                                    gCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCt                                                                                   <  2:445647/150‑1 (MQ=255)
                                                                       gTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGc                                                                                <  1:4974945/150‑1 (MQ=255)
                                                                         cGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCaa                                                                             <  2:3586077/150‑1 (MQ=255)
                                                                          gACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAg                                                                            <  1:3786453/150‑1 (MQ=255)
                                                                             gCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGGTGGTGTAGCTCACAGAACGAGGTGCAGCTGGCAAGCTc                                                                         >  2:2468816/1‑150 (MQ=255)
                                                                              ccGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCa                                                                        <  1:881716/150‑1 (MQ=255)
                                                                              ccGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCAGGCGGCCGGGTTGGCCGCCAACCGTTGGTGTAGCTCACTGAAAGAGGTGCAGCTGGCCAGcaaa                                                                        >  1:324211/1‑147 (MQ=255)
                                                                               cGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAg                                                                       <  1:1478284/150‑1 (MQ=255)
                                                                               cGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAg                                                                       <  2:4201013/150‑1 (MQ=255)
                                                                                gAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGc                                                                      >  2:5331089/1‑150 (MQ=255)
                                                                                gAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCAGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGc                                                                      >  2:5298016/1‑150 (MQ=255)
                                                                                      ttCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCAGCACGGCCGCGGCAGAGCTGGCGGCCGGGGTGGCCGCCAAGCGGTGGGGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCgctg                                                                >  1:1873640/1‑150 (MQ=255)
                                                                                         ggCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCtc                                                             <  2:4811803/150‑1 (MQ=255)
                                                                                         ggCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCtc                                                             <  1:3543751/150‑1 (MQ=255)
                                                                                            ttcttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCtcgtc                                                          <  2:6213930/150‑1 (MQ=255)
                                                                                               ttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCt                                                       >  2:1534689/1‑150 (MQ=255)
                                                                                               ttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCt                                                       >  2:5051745/1‑150 (MQ=255)
                                                                                                tcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTc                                                      <  1:1420282/150‑1 (MQ=255)
                                                                                                 ccAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCa                                                     <  1:3287787/150‑1 (MQ=255)
                                                                                                 ccAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCa                                                     <  1:5281642/150‑1 (MQ=255)
                                                                                                   aGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGt                                                   <  1:3667956/150‑1 (MQ=255)
                                                                                                   aGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGt                                                   <  2:6491217/150‑1 (MQ=255)
                                                                                                   aGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGt                                                   <  2:3345583/150‑1 (MQ=255)
                                                                                                   aGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGGTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCACGCTCAGCCCGCTGCGCGTCGCTCAgg                                                   >  2:1422220/1‑149 (MQ=255)
                                                                                                     cTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTgg                                                 >  2:5995559/1‑150 (MQ=255)
                                                                                                       cggcggCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAGCGACGAGCAGCGGGCTGAGCTTGCCAGCTGCACCTCGTTCAGTgagc                                               <  1:4142275/150‑4 (MQ=21)
                                                                                                          cggcgCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGt                                            <  2:6664879/150‑1 (MQ=255)
                                                                                                              gcTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGc                                        <  1:2342137/150‑1 (MQ=255)
                                                                                                              gcTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGc                                        <  1:2996473/150‑1 (MQ=255)
                                                                                                                tGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCa                                      >  2:3913656/1‑150 (MQ=255)
                                                                                                                 gCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCATAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAg                                     <  1:4574704/150‑1 (MQ=255)
                                                                                                                 gCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAg                                     <  1:1357472/150‑1 (MQ=255)
                                                                                                                 gCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAg                                     <  1:785457/150‑1 (MQ=255)
                                                                                                                  cccACGGTCAGGCACGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAgg                                    >  1:3317274/1‑150 (MQ=255)
                                                                                                                   ccACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGGTGGGGGAGCTCACTGAACGAGGTGCAGCTGGCAAGCACAGCCCGCTGCTCGTCGCTCCGTGGTCGGTGGGCCAGGt                                   >  2:3713852/1‑150 (MQ=255)
                                                                                                                        gTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTTGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCg                              >  2:67977/1‑150 (MQ=255)
                                                                                                                           aGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGGTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCa                           >  2:4021837/1‑150 (MQ=255)
                                                                                                                               gcgGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCCCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCCCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCg                       >  2:1170893/1‑150 (MQ=255)
                                                                                                                                     ccACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg                 >  1:4870709/1‑150 (MQ=255)
                                                                                                                                     ccACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg                 >  2:5717855/1‑150 (MQ=255)
                                                                                                                                      tgccGAGGCCGTGCAGGCCTCGGTAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCgg                >  1:6157915/3‑150 (MQ=255)
                                                                                                                                      tgccGAGGCCGTGCAGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCtcgtcg                                                         <  1:5975049/107‑1 (MQ=255)
                                                                                                                                      cACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCgg                <  2:3536804/150‑1 (MQ=255)
                                                                                                                                       aCCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGc               <  1:6435188/150‑1 (MQ=255)
                                                                                                                                        ccGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGcc              >  2:5892051/1‑150 (MQ=255)
                                                                                                                                          gAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCtt            <  1:1625148/150‑1 (MQ=255)
                                                                                                                                               cTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCt       >  2:5525618/1‑150 (MQ=255)
                                                                                                                                                   cACGGCGTAGGCAGAGCTGGCGGCGGGTGTGGCCGCCAAGCGTTGGTGTATCTCACTGAACGAGGTGCAGCTGGAAAGCTCAGCACTCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTTCTCGCCAAGCGATTGCGGCCTTGCGCTTGa    <  1:5421456/150‑1 (MQ=255)
                                                                                                                                                    aCGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAg   >  1:5693576/1‑150 (MQ=255)
                                                                                                                                                     cGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTATCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAGt  <  1:2677450/150‑1 (MQ=255)
                                                                                                                                                     cGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTATCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAGt  <  2:4152958/150‑1 (MQ=255)
                                                                                                                                                     cGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAGt  <  2:1738038/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTTGCGCTTGAGT  >  NC_002947/5731666‑5731963

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