Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,101,298 G→A 6.6% A110V (GCA→GTA)  phhR ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,101,2980GA 6.6% 176.9 / 2.1 60A110V (GCA→GTA) phhRsigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

AACACATCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCGGCAATG  >  NC_002947/5101151‑5101441
                                                                                                                                                   |                                                                                                                                                    
aaCACATCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTg                                                                                                                                                      >  1:335956/1‑148 (MQ=255)
 acacATCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCa                                                                                                                                                    <  5:127563/149‑1 (MQ=255)
   acaTCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACg                                                                                                                                                  >  6:95190/1‑149 (MQ=255)
      tCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGtt                                                                                                                                               >  8:87574/1‑149 (MQ=255)
         ccTTTGACTTTGATGCGCAAGCCTTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGc                                                                                                                                            <  6:49203/149‑1 (MQ=255)
         ccTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGc                                                                                                                                            <  3:280820/149‑1 (MQ=255)
          ctttgactttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTgg                                                                                                                                             <  5:224244/147‑1 (MQ=255)
          ctttgactttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTgg                                                                                                                                             >  6:224244/1‑147 (MQ=255)
          ctttgactttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGcc                                                                                                                                           <  2:335956/149‑1 (MQ=255)
           tttgactttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCg                                                                                                                                          >  4:303700/1‑149 (MQ=255)
             tgactttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGc                                                                                                                                         <  6:251453/148‑1 (MQ=255)
                 tttgaTGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGa                                                                                                                                    >  1:77318/1‑149 (MQ=255)
                     aTGCGCAAGCCGTTGATGCGCGTCTTGTTGGCCGGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgc                                                                                                                                 <  7:325697/148‑1 (MQ=255)
                      tGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgctg                                                                                                                               >  5:327331/1‑149 (MQ=255)
                                gTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGt                                                                                                                     <  5:95190/149‑1 (MQ=255)
                                 ttGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGtc                                                                                                                    >  5:115342/1‑149 (MQ=255)
                                    aTGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGtcga                                                                                                                  <  5:260853/148‑1 (MQ=255)
                                        gcgACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGt                                                                                                                     >  2:50629/1‑141 (MQ=255)
                                        gcgACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGt                                                                                                                     <  1:50629/141‑1 (MQ=255)
                                          gACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgctg                                                                                                                               >  1:295962/1‑129 (MQ=255)
                                          gACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgctg                                                                                                                               <  2:295962/129‑1 (MQ=255)
                                           aCTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTACACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGc                                                                                                          >  3:228069/1‑149 (MQ=255)
                                           aCTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTACACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGc                                                                                                          >  3:301464/1‑149 (MQ=255)
                                               gttgGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCAcc                                                                                                       <  3:303700/148‑1 (MQ=255)
                                               gttgGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCAc                                                                                                        <  5:87965/147‑1 (MQ=255)
                                                  ggCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCggg                                                                                                    <  7:87574/148‑1 (MQ=255)
                                                    cccGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAg                                                                                                 >  2:57058/1‑149 (MQ=255)
                                                    cccGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAg                                                                                                 <  6:115342/149‑1 (MQ=255)
                                                     ccGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGt                                                                                                >  2:238717/1‑149 (MQ=255)
                                                               tCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGaa                                                                                                  >  6:58924/1‑137 (MQ=255)
                                                               tCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGaa                                                                                                  <  5:58924/137‑1 (MQ=255)
                                                                 cGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGccccc                                                                                     <  1:238717/148‑1 (MQ=255)
                                                                         cGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgc                                                                                                                                 <  4:12715/96‑1 (MQ=255)
                                                                         cGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGAtgc                                                                                                                                 >  3:12715/1‑96 (MQ=255)
                                                                         cGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAg                                                                                   >  7:298608/1‑142 (MQ=255)
                                                                         cGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAg                                                                                   <  8:298608/142‑1 (MQ=255)
                                                                          gAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑At                                                                      >  7:221067/1‑149 (MQ=255)
                                                                           aaGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑Att                                                                     >  3:69865/1‑149 (MQ=255)
                                                                                     cGTAGCGCGCCAGCGGCATGCCCGCCACCTCCTCACCCCGCACGCCAACCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGCTCGACAGCACCGGGAAGTCCAGCGCCCCCACCAA‑‑CG‑‑C‑ATTCAGCTCCa                                                             <  6:327331/147‑1 (MQ=255)
                                                                                       tAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCa                                                                   <  1:126179/139‑1 (MQ=255)
                                                                                       tAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCa                                                                   >  2:126179/1‑139 (MQ=255)
                                                                                        aGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGc                                                        >  4:167137/1‑149 (MQ=255)
                                                                                          cgcgcCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCg                                                      <  2:77318/149‑1 (MQ=255)
                                                                                                          ccGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCAt                                       >  4:7016/1‑148 (MQ=255)
                                                                                                            ggCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAg                                    <  8:221067/149‑1 (MQ=255)
                                                                                                             gCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑CAATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAg                                    <  1:57058/149‑1 (MQ=255)
                                                                                                               aCCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGcc                                 <  3:7016/149‑1 (MQ=255)
                                                                                                                 cTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑Att                                                                     <  1:174119/111‑1 (MQ=255)
                                                                                                                 cTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑Att                                                                     >  2:174119/1‑111 (MQ=255)
                                                                                                                 cTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCa                                        >  2:264100/1‑140 (MQ=255)
                                                                                                                 cTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCa                                        <  1:264100/140‑1 (MQ=255)
                                                                                                                     tcaACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACg                           >  2:217662/1‑149 (MQ=255)
                                                                                                                       aaCCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAATTGGTAG‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCcaca                             >  1:287513/1‑149 (MQ=255)
                                                                                                                        aCCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCtt                        <  4:69865/149‑1 (MQ=255)
                                                                                                                             cACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGc                   >  3:73442/1‑149 (MQ=255)
                                                                                                                             cACGCCAAGCAACTGCGCGGCTGCACCGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGc                   >  8:30112/1‑149 (MQ=255)
                                                                                                                                  cAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGaa               >  1:210008/1‑148 (MQ=255)
                                                                                                                                         cTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCgg       <  5:292438/149‑1 (MQ=255)
                                                                                                                                              cgGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCGGCAATg  <  4:73442/149‑1 (MQ=255)
                                                                                                                                              cgGCTACACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCGGCAATg  <  4:228069/149‑1 (MQ=255)
                                                                                                                                              cgGCTACACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCGGCAATg  <  4:301464/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACACATCACCTTTGACTTTGATGCGCAAGCCGTTGATGCGCGACTTGTTGGCCCGCACCAACTCCGGAAGATCGAAGTCTTCTACGTAGCGCGCCAGCGGCATGCCCGGCACCTCATCAACCCGCACGCCAAGCAACTGCGCGGCTGCACGGTTGGCCGCAACGATGCTGCCGCCCATGTCGATCGACAGCACCGGGAAGTCCAGCGCCCCCAGCAA‑‑CG‑‑C‑ATTCAGCTCCATGTGCCGGCGTTCGCTGGGCATCAGGCCCACACGCTTAACGCCGAATACCCCGGCAATG  >  NC_002947/5101151‑5101441

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: ‑