Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,955,828 G→C 5.2% intergenic (+54/+150) 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,8280GC 5.2% 251.4 / 3.8 77intergenic (+54/+150)elbB/PP_5221putative hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base C (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCT  >  NC_002947/5955687‑5955976
                                                                                                                                             |                                                                                                                                                    
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)
                                                                                                                                             |                                                                                                                                                    
GGCGTACATGGAGGCCAAGTCGATCAGCGAGGCGGCGGGCGGGATCTACAAGCTGGTGGACCGAGTGCTGGAGCTGACGCACGAGTAACTTACTGCGTTGGCTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCGCAGGCTTCAAGGCCACGCGGTACCTGTGGGGCAACTGCGCGATAGCCTGTGGGGCAACTGCGCGATCGCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGCCAGCCTGGGCTGTCGCCTACAGGAATCT  >  NC_002947/5955687‑5955976

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: