Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,726,860:1 +G 6.4% coding (162/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,726,8601.G 6.4% 322.9 / 14.8 78coding (162/2574 nt)opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base . (36/37);  new base G (3/2);  total (39/39)
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.09e-01

GCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTT  >  NC_002947/5726716‑5727004
                                                                                                                                                 |                                                                                                                                                
gcgcCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGccgtgca                                                                                                                                              <  1:1555819/150‑6 (MQ=255)
   ccATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGATCAGGCGCGGGCCCACCGAGGCCTGCACGGc                                                                                                                                          <  2:738668/150‑1 (MQ=255)
      tGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                       <  2:1090330/150‑1 (MQ=255)
        gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCCGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                     >  1:21628/1‑150 (MQ=255)
        gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                     <  1:931576/150‑1 (MQ=255)
        gCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCgg                                                                                                                                     >  2:1248614/1‑150 (MQ=255)
          tCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCa                                                                                                                                   <  2:654995/150‑1 (MQ=255)
           cGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCGCCGAGGCCGTGCAGGCCTCggtgg                                                                                                                                  <  1:1354714/150‑4 (MQ=255)
            gggCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCaga                                                                                                                                 <  1:1719828/150‑1 (MQ=255)
                  gACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTggc                                                                                                                           >  2:1024966/1‑150 (MQ=255)
                         gTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCggg                                                                                                                    >  2:406359/1‑150 (MQ=255)
                             gTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTgg                                                                                                                >  2:324224/1‑150 (MQ=255)
                              tttGATCGGCGGGGCGATCTTCATGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGgc                                                                                                               >  2:1267714/1‑150 (MQ=255)
                              tttGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGgc                                                                                                               >  1:538742/1‑150 (MQ=255)
                              tttGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGgc                                                                                                               >  2:1890143/1‑150 (MQ=255)
                                  aTCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGgccgcc                                                                                                           >  2:1831045/1‑150 (MQ=255)
                                         gggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGtt                                                                                                    <  2:1216794/150‑1 (MQ=255)
                                          ggCGATCTTCATGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTg                                                                                                   <  1:1267714/150‑1 (MQ=255)
                                          ggCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCACACCGAGGCCTTCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTg                                                                                                   >  1:1887219/1‑150 (MQ=255)
                                          ggCGATCTTCAGGCACAGGCGGCCGCTGCCGCCGACGCCGAACATTTCGTCTTCGTCCAGCTCGGCGGCGCTGCCCACGGACAGGCGCGCGCCCACCGAGACCTACACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTg                                                                                                   <  2:544282/150‑1 (MQ=255)
                                           gCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                       >  2:1000414/1‑113 (MQ=255)
                                           gCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTc                                                                                                                                       <  1:1000414/113‑1 (MQ=255)
                                             gATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCGGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgtg                                                                                                >  1:1185342/1‑150 (MQ=255)
                                              aTCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgtgt                                                                                               <  1:1890143/150‑1 (MQ=255)
                                               tCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTa                                                                                              <  1:944530/150‑1 (MQ=255)
                                                 ttCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGc                                                                                            <  1:1363709/150‑1 (MQ=255)
                                                 ttCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGc                                                                                            >  2:1363709/1‑150 (MQ=255)
                                                      gCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCGGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACt                                                                                       <  2:1185342/150‑1 (MQ=255)
                                                        acaGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGa                                                                                     <  1:1024966/150‑1 (MQ=255)
                                                        acaGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGa                                                                                     <  2:1344076/150‑1 (MQ=255)
                                                                    tGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgt                                                                                                 >  2:1307433/1‑126 (MQ=255)
                                                                    tGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGgt                                                                                                 <  1:1307433/126‑1 (MQ=255)
                                                                    tGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGc                                                                         <  1:324224/150‑1 (MQ=255)
                                                                     gCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCt                                                                        <  2:1673470/150‑1 (MQ=255)
                                                                        gTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGc                                                                     >  1:1445047/1‑150 (MQ=255)
                                                                        gTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGc                                                                     <  2:78601/150‑1 (MQ=255)
                                                                          cGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGATGTGCAGCTGGCaa                                                                   <  2:1734318/150‑1 (MQ=255)
                                                                            aCGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGc                                                                 >  2:1252539/1‑150 (MQ=255)
                                                                                   gCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGccc                                                          <  2:1299179/150‑1 (MQ=255)
                                                                                    cATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCt                                                                        <  1:1864357/135‑1 (MQ=255)
                                                                                    cATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCt                                                                        >  2:1864357/1‑135 (MQ=255)
                                                                                       ttCGGCTTCTTCCAGCTGGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAg                                                                                >  2:996312/1‑124 (MQ=255)
                                                                                       ttCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAg                                                                                <  1:996312/124‑1 (MQ=255)
                                                                                       ttCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCCGGCAAGCTCAGCCCgctg                                                      >  1:1236048/1‑150 (MQ=255)
                                                                                        tCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCgctgc                                                     <  2:1409804/150‑1 (MQ=255)
                                                                                                ttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACt                                                                                       <  1:1837546/108‑1 (MQ=255)
                                                                                                ttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACt                                                                                       >  2:1837546/1‑108 (MQ=255)
                                                                                                ttcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCt                                             >  2:1369429/1‑150 (MQ=255)
                                                                                                atcCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTTAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCt                                             >  2:1327461/2‑150 (MQ=255)
                                                                                                  ccAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCa                                           >  2:327210/1‑150 (MQ=255)
                                                                                                   cAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAg                                          >  1:478945/1‑150 (MQ=255)
                                                                                                      cTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGCACAGTAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCGGAGTgg                                       >  1:1476603/1‑150 (MQ=255)
                                                                                                      cTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTgg                                       <  1:1518401/150‑1 (MQ=255)
                                                                                                       tCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGt                                      >  2:1593956/1‑150 (MQ=255)
                                                                                                       tCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCCCGTCGCTCAGTGGt                                      >  1:881654/1‑150 (MQ=255)
                                                                                                           cggcgCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTgg                                                                                                                <  1:1758838/72‑1 (MQ=255)
                                                                                                           cggcgCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTgg                                                                                                                >  2:1758838/1‑72 (MQ=255)
                                                                                                               gcTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGc                              >  2:365190/1‑150 (MQ=255)
                                                                                                               gcTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGc                              >  1:177596/1‑150 (MQ=255)
                                                                                                                cTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGcc                             <  2:1445047/150‑1 (MQ=255)
                                                                                                                 tGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCa                            >  2:546137/1‑150 (MQ=255)
                                                                                                                     cACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCa                                                                                                                                   >  2:925296/1‑43 (MQ=255)
                                                                                                                     cACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCa                                                                                                                                   <  1:925296/43‑1 (MQ=255)
                                                                                                                      aCGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTAc                       <  2:1377570/150‑1 (MQ=255)
                                                                                                                      aCGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTAc                       <  1:1369429/150‑1 (MQ=255)
                                                                                                                           aaGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGcc                  <  1:1145418/149‑1 (MQ=255)
                                                                                                                              gcgcGGGCCCACCGCGGCCTGCAAGGCATCGGCCGGGCTGGCGGCCGCACTGTCCTCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCGAg               <  2:1476603/150‑1 (MQ=255)
                                                                                                                                  agctctgccGAGGCCGTGCAGGCCTCGGTAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGATGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGAt           >  1:900968/8‑150 (MQ=255)
                                                                                                                                    gCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTc                     >  1:1073520/1‑138 (MQ=255)
                                                                                                                                    gCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTc                     <  2:1073520/138‑1 (MQ=255)
                                                                                                                                      ctgccGAGGCCGTGCAGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg       >  1:1302320/4‑150 (MQ=255)
                                                                                                                                      ccACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg       >  1:1396015/1‑150 (MQ=255)
                                                                                                                                      ccACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg       >  2:304942/1‑150 (MQ=255)
                                                                                                                                      ccACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCGCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCg       >  1:437107/1‑150 (MQ=255)
                                                                                                                                       tgccGAGGCCGTGCAGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGa                                                                                 >  2:548977/3‑75 (MQ=38)
                                                                                                                                       tgccGAGGCCGTGCAGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGa                                                                                 <  1:548977/73‑1 (MQ=38)
                                                                                                                                         ccGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGcc    <  1:1593956/150‑1 (MQ=255)
                                                                                                                                          cGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCCGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCt   <  2:1236048/150‑1 (MQ=255)
                                                                                                                                           gAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCtt  <  2:21628/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGCCATGGCTCGGGCACGACCTTGGTGCGTTTGATCGGCGGGGCGATCTTCAGGCACAGGCGGCCGCTGCCGTCGACGCCGAGCATTTCGGCTTCTTCCAGCTCGGCGGCGCTGCCCACGGTCAGGCGCGGGCCCACCGAGGCCTGCACGGCCTCGGCAGAGCTGGCGGCCGGGTTGGCCGCCAAGCGTTGGTGTAGCTCACTGAACGAGGTGCAGCTGGCAAGCTCAGCCCGCTGCTCGTCGCTCAGTGGTAGGTGGGCCAGGTACTCGCCAAGCGATTCCGGCCTT  >  NC_002947/5726716‑5727004

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