Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,955,832 C→G 6.1% intergenic (+58/+146) elbB → / ← PP_5221 putative hydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,955,8320CG 6.1% 254.6 / 4.1 79intergenic (+58/+146)elbB/PP_5221putative hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base G (3/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

GGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTG  >  NC_002947/5955687‑5955977
                                                                                                                                                 |                                                                                                                                                 
ggCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAg                                                                                                                                                >  2:228832/1‑149 (MQ=255)
     aCATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAc                                                                                                                                           >  7:32906/1‑149 (MQ=255)
       aTGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCAcgc                                                                                                                                         >  1:386096/1‑149 (MQ=255)
         ggAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCgg                                                                                                                                       >  4:78431/1‑149 (MQ=255)
          gAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCgg                                                                                                                                       >  8:83642/1‑148 (MQ=255)
            ggCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTAc                                                                                                                                    >  5:131728/1‑149 (MQ=255)
            ggCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTAc                                                                                                                                    <  3:353949/149‑1 (MQ=255)
               cAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCtg                                                                                                                                 >  1:50524/1‑149 (MQ=255)
                  gTCGATCAGCGAGGCGGCGGGCGGTATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTgg                                                                                                                              <  5:234300/149‑1 (MQ=255)
                  gTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTgg                                                                                                                              <  3:95600/149‑1 (MQ=37)
                    cGAGCAGCGAGGCGGCGGGCGGGATCTACAAGCGGGGGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTggga                                                                                                                            <  2:342147/149‑2 (MQ=25)
                     gATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTgggg                                                                                                                            <  7:251419/148‑1 (MQ=255)
                            gatgcggcgGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTAACTGTGGGGCAACTgcg                                                                                                                    >  5:363666/4‑149 (MQ=255)
                                 ggcgggcgGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATa                                                                                                               <  5:7144/149‑1 (MQ=255)
                                  gcgggcggGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAg                                                                                                              >  6:345927/1‑149 (MQ=255)
                                        gggATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCtgtg                                                                                                        >  7:240962/1‑149 (MQ=255)
                                          gATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAATGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTggg                                                                                                      >  2:353688/1‑149 (MQ=255)
                                           aTCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGGGGGGCAACTGAGCGATAGCCGGTgggg                                                                                                     >  2:87394/1‑149 (MQ=255)
                                              tACAAGCTGGTGGACCGAGTGCTGGAGCAGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCACGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCaa                                                                                                  >  8:130239/1‑149 (MQ=255)
                                               aCAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAAc                                                                                                 >  8:296972/1‑149 (MQ=255)
                                               aCAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAAc                                                                                                 >  6:221640/1‑149 (MQ=255)
                                                cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACt                                                                                                <  3:78431/149‑1 (MQ=255)
                                                cAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACt                                                                                                >  8:172153/1‑149 (MQ=255)
                                                     tggtggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGcgcg                                                                                           <  2:386096/149‑1 (MQ=255)
                                                       gtggACGGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCGTGAAGGGCCCGCGGTGCCGGTGGGGCCACTGCGCGATCGCCTGTGGGGCAACGGCGCGAt                                                                                         >  4:101854/1‑149 (MQ=37)
                                                        tggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                        <  5:278106/149‑1 (MQ=255)
                                                        tggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                        <  8:222365/149‑1 (MQ=255)
                                                        tggACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATc                                                                                        <  7:111758/149‑1 (MQ=255)
                                                              gAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTg                                                                                   >  5:88879/1‑148 (MQ=255)
                                                              gAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGt                                                                                  <  1:154599/149‑1 (MQ=255)
                                                               aGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTg                                                                                   >  7:71861/1‑147 (MQ=255)
                                                                gTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAg                                                                                <  6:131728/149‑1 (MQ=255)
                                                                 tGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAg                                                                                >  3:148154/1‑148 (MQ=255)
                                                                     ggAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCg                                                                           >  2:120287/1‑149 (MQ=255)
                                                                      gAGCTGACGCACGAGTAACTTACTGCGTTTGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCg                                                                           <  8:344504/148‑1 (MQ=255)
                                                                      gAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCgg                                                                          <  2:50524/149‑1 (MQ=255)
                                                                      gAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCgg                                                                          <  4:148154/149‑1 (MQ=255)
                                                                      gAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCgg                                                                          <  7:172153/149‑1 (MQ=255)
                                                                      gAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTCGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCgg                                                                          >  1:319125/1‑149 (MQ=255)
                                                                              gCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCGGGGGGGCAACTGCGCGATCGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTAccc                                                                  >  7:207707/1‑149 (MQ=255)
                                                                                aCGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCCGCGCAGGCTTCAAGGCCACGCGGTAGCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTGTACCcgc                                                                >  6:313196/1‑149 (MQ=255)
                                                                                aCGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAAATGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCcgc                                                                >  4:365299/1‑149 (MQ=255)
                                                                                  gAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCcgcg                                                               >  4:352309/1‑148 (MQ=255)
                                                                                       acttacttCGTTGGCTTCTTCGCGGGGATGCCCGGTCCCACAGGTCCAGCGCATGCTTCAAGGGCACGCGGTACCTTTGGGGCAACTGCGCGATATCCTTTGGGGCAACTGCGCGATCGCTTGTAGGAg                                                                             <  6:308343/129‑1 (MQ=255)
                                                                                       aCTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAg                                                                             >  5:308343/1‑129 (MQ=255)
                                                                                       aCTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAacac                                                          <  8:32906/148‑1 (MQ=255)
                                                                                         ttACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCgg                                                       >  5:363460/1‑149 (MQ=255)
                                                                                                gTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGcc                                                <  5:345927/149‑1 (MQ=255)
                                                                                                    gCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCt                                                                                                           <  7:70060/86‑1 (MQ=255)
                                                                                                    gCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCt                                                                                                           >  8:70060/1‑86 (MQ=255)
                                                                                                     cttcttCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGc                                           <  4:141314/149‑1 (MQ=255)
                                                                                                       tcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCTCTGCCGGTgtt                                          <  8:240962/148‑3 (MQ=255)
                                                                                                       tcttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCa                                         >  2:40918/1‑149 (MQ=255)
                                                                                                         ttcGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATc                                       >  1:15774/1‑149 (MQ=255)
                                                                                                           cgcgGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTTTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATc                                       <  6:363666/147‑1 (MQ=255)
                                                                                                           cgcgGGCATGCCCGCGCCCACAGGTCCAGCGCAGGCGTCAAGGCCACGCGGTACCTGGGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGAGCGCTTGTAGGCGCGGGTTTACCCGCGAACACCGGCTCTGCCGGTgttcgc                                       <  3:304038/147‑6 (MQ=255)
                                                                                                            gcgGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCAc                                    >  3:192542/1‑149 (MQ=255)
                                                                                                               ggCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACcgc                                 <  1:228832/149‑1 (MQ=255)
                                                                                                                gCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACcgcg                                >  1:149481/1‑149 (MQ=255)
                                                                                                                     cccGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGc                           <  1:120287/149‑1 (MQ=255)
                                                                                                                        gCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTg                        >  7:239483/1‑149 (MQ=255)
                                                                                                                            ccACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTggg                      <  2:319125/147‑1 (MQ=255)
                                                                                                                            ccACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCt                    >  8:105764/1‑149 (MQ=255)
                                                                                                                            ccACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCt                    >  2:237820/1‑149 (MQ=255)
                                                                                                                            ccACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCt                    >  2:232357/1‑149 (MQ=255)
                                                                                                                            ccACAGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCt                    >  7:317953/1‑149 (MQ=25)
                                                                                                                             cacaGGTACCGCGTGGCCTTGAAGCCTGCGCTGGACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTg                   >  3:99771/1‑149 (MQ=25)
                                                                                                                              acaGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGt                  >  4:174601/1‑149 (MQ=255)
                                                                                                                               caGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTc                 <  2:15774/149‑1 (MQ=255)
                                                                                                                                aGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCg                <  7:83642/149‑1 (MQ=255)
                                                                                                                                 gggccAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGc               <  5:221640/146‑1 (MQ=255)
                                                                                                                                  gTCCAGCGCAGGCTTCAATGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGACAGCCCAGGCTGgcgcg              <  1:353688/149‑2 (MQ=38)
                                                                                                                                     cAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTa            <  8:230460/148‑1 (MQ=255)
                                                                                                                                      aGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACa          <  3:352309/149‑1 (MQ=255)
                                                                                                                                           aGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCtt                                                                                    >  1:162330/1‑70 (MQ=255)
                                                                                                                                           aGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCtt                                                                                    <  2:162330/70‑1 (MQ=255)
                                                                                                                                           aGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAAt     <  6:363460/149‑1 (MQ=255)
                                                                                                                                             gCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCt   >  1:318442/1‑149 (MQ=255)
                                                                                                                                              cTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTg  <  7:296972/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCTG  >  NC_002947/5955687‑5955977

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: