Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,984,447 G→A 16.4% intergenic (+43/+46) PP_5242 → / ← ycaC‑II GAF domain‑containing sensor histidine kinase/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,984,4470GA16.4% 132.4 / 18.2 61intergenic (+43/+46)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (3/7);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

TCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCG  >  NC_002947/5984301‑5984595
                                                                                                                                                   |                                                                                                                                                     
tCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                                     >  2:37500/1‑149 (MQ=255)
    tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGTCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGAGTAAACCCGCTCCTACGGTCGGCTGTa                                                                                                                                                  <  5:217455/148‑1 (MQ=255)
    tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGCGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAg                                                                                                                                                 <  1:340087/149‑1 (MQ=255)
       ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAgggg                                                                                                                                              <  5:189232/149‑1 (MQ=255)
       ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAgggg                                                                                                                                              >  8:112439/1‑149 (MQ=255)
         cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCTGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTTTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                             <  6:213619/148‑1 (MQ=37)
          ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                           <  2:32943/149‑1 (MQ=38)
          ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                           <  2:27229/149‑1 (MQ=38)
          ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCg                                                                                                                                            <  3:235331/148‑1 (MQ=38)
           gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                                          <  4:179433/149‑1 (MQ=38)
                tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTTCGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                                     >  8:340021/1‑149 (MQ=255)
                tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                                     >  4:27850/1‑149 (MQ=255)
                  aGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                                   <  1:37500/149‑1 (MQ=255)
                  aGGCCCACGGCGGGCAGCTGCAGGCGCTCGTGACGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCGGCCTTTGATTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTATTCCCGCTGCTACGGTCGGGGGTAGGGGCGGGTTTAccc                                                                                                                                   >  1:46756/1‑149 (MQ=255)
                      ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGGTCCTACGGGCGGCGGGAGGGGCGGGTTTTCCcgcg                                                                                                                               >  3:43443/1‑148 (MQ=39)
                      ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcgcg                                                                                                                               >  2:104658/1‑148 (MQ=255)
                         cggcggGCAGCGGCCGGCGCTGGTGTCGCAGGTGGGGGGAGGCCGGATTCGTGTCAGGCTGCCTTTGAAGCACGCATAGCCTGTGCCGGCCTCTGCGCGGGTAAACCCGCTCCTCCGGGCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                           <  3:27850/149‑1 (MQ=37)
                                        ggCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGttt                                                                                                            <  7:340021/149‑1 (MQ=255)
                                                  tCGCAGTTGGGTGGATGCCGGATTCGGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAGGTCGGGTGTAAGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTAAATAAgcg                                                                                                   >  1:197245/1‑149 (MQ=37)
                                                          gggTGCATGCCGGATTCGTGTCAGTCCGCCTTTGAATCACGCATAGCCTGTGCCTGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCAGTAGGGGCGGGTTTACCCGAGAAGAGGCCAGTACAGGTTTAAATCAGCCCGAGGCCac                                                                                          >  1:41948/1‑148 (MQ=38)
                                                                    cgtaTTGGTGTCTGTCTGCATTTGAATCCCGAATTGCCGGTGCGGGCGTCTTCGCGGGTAAACCCGCTCCTAGGGGGGGTTTTTGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAg                                                                                <  2:197245/146‑1 (MQ=25)
                                                                    ccgATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCATCGCGAGGCCAGCGCTGCGTAg                                                                                >  7:174704/3‑149 (MQ=255)
                                                                    cGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGttt                                                                                                                                       >  1:235905/1‑95 (MQ=255)
                                                                    cGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGttt                                                                                                                                       <  2:235905/95‑1 (MQ=255)
                                                                                 aGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAgg                                                                    <  6:105205/149‑1 (MQ=255)
                                                                                      cccttgtaaTCACGCATCCCTTGGGCGGGCTTCTTCGCGGGTAACCCCGCTCCTAGGGGCGGCTGTTGGGGGGGTTTAACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGc                                                                <  8:74125/141‑1 (MQ=255)
                                                                                        cTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCgg                                                              >  5:353578/1‑148 (MQ=255)
                                                                                         tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCGGCGACGGTGGGCGGGAGGGGCGGGGTTTCCCGCGAAGAGGGCAGTACAGGGTTACATAAGCGCGAGGCCCGCGCGGGGGTGCTGTTGATGAGGGTggggga                                                            >  8:341138/1‑147 (MQ=255)
                                                                                              aTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                       <  7:112439/149‑1 (MQ=255)
                                                                                                   gCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTGCAACGGTCGGCGGTAGGGGCGGGGTTTCCCGCGAAGACGCCAGGACAGGTTTACATGAGCGCGAGGCCACCGCGGCGTAGCTGTTCAAGAGGTTGGGGTAGTTGGGGATa                                                  >  7:100090/1‑149 (MQ=37)
                                                                                                   gCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                                      <  2:114749/129‑1 (MQ=255)
                                                                                                   gCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                                      >  1:114749/1‑129 (MQ=255)
                                                                                                    cATAGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATAc                                                 <  3:301713/149‑1 (MQ=255)
                                                                                                        gCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCt                                              <  8:264005/148‑1 (MQ=255)
                                                                                                         ccTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCACTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                            >  3:287941/1‑149 (MQ=255)
                                                                                                           tgtcCCGGCCGCTTGGCGGGTAACCCCGCTCCTACGCGGGGGGGTAGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGTTGCGa                                          <  8:100090/149‑1 (MQ=38)
                                                                                                                cGGCCTCTTCGGGGGTATACCCGCTCCTACGGTGGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagca                                     <  8:89814/149‑1 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTCCCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGCTACGCTGCGACAGCAGGTTAcc                              >  3:234853/1‑149 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCTGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGGGTAGCTGTTCAGCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAc                               >  1:89551/1‑149 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc                         <  7:341138/149‑1 (MQ=255)
                                                                                                                            gggTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc                         >  2:45133/1‑149 (MQ=25)
                                                                                                                              gTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGATGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCt                       >  6:253981/1‑149 (MQ=255)
                                                                                                                               tAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                      <  1:45133/149‑1 (MQ=37)
                                                                                                                                   gccGCTCGTACGGTCGGCTGTAGGGGCGGGTTTACCCTCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCTTCAGGTTGGGGTAGTTGGGGATACGGTGCGACAGCAGGTTACCCAGCCCTTCGAt                  <  2:46756/148‑1 (MQ=255)
                                                                                                                                   cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGGTACCCAGCCCTGCggt                  >  8:261283/1‑147 (MQ=255)
                                                                                                                                   cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCACTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt                  <  4:287941/149‑1 (MQ=255)
                                                                                                                                   accGCTCCTACGTTCGGGTTTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTTCGTAGCTGTTCATCAGGTTGCGGTAGTTGTGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt                  <  2:41948/148‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACGGTAGGCTGTAGCGGCGGCTTTACCCGCGAAGAGGCCAGTACAGTTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATAAGCTGAGACAGCAGGTTACCCCGCCCTTCGATg                 >  6:147029/1‑149 (MQ=255)
                                                                                                                                    ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCg                    >  4:279667/1‑146 (MQ=38)
                                                                                                                                    ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCg                    <  3:279667/146‑1 (MQ=38)
                                                                                                                                     cGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt                >  1:188823/1‑149 (MQ=255)
                                                                                                                                     cGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt                >  3:47051/1‑149 (MQ=38)
                                                                                                                                      gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc               <  1:104658/149‑1 (MQ=255)
                                                                                                                                         ccTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGtt            >  1:18330/1‑149 (MQ=255)
                                                                                                                                            aCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgcg         >  7:162863/1‑149 (MQ=255)
                                                                                                                                             cGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgcgc        <  2:188823/149‑1 (MQ=255)
                                                                                                                                             cGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgcgc        <  3:162582/149‑1 (MQ=255)
                                                                                                                                             cGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc          <  4:43443/147‑1 (MQ=255)
                                                                                                                                             cGGGCGGCGGTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc          <  6:353578/147‑1 (MQ=255)
                                                                                                                                               gTCGGCTGTAGGGTCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCa      <  5:253981/149‑1 (MQ=255)
                                                                                                                                                tCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAg     >  4:32458/1‑149 (MQ=255)
                                                                                                                                                   gCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTcg  >  8:202862/1‑149 (MQ=255)
                                                                                                                                                   gCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTcg  >  8:208377/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCAGTCG  >  NC_002947/5984301‑5984595

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