Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,984,498 G→A 9.0% intergenic (+38/+51) 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,4980GA 9.0% 249.9 / 9.5 88intergenic (+38/+51)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (45/35);  new base A (3/5);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCC  >  NC_002947/5984349‑5984646
                                                                                                                                                     |                                                                                                                                                    
cTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACg                                                                                                                                                      <  1:539564/150‑1 (MQ=255)
        tCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTg                                                                                                                                              <  1:2711586/150‑1 (MQ=255)
           gTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAg                                                                                                                                           <  2:1305104/150‑1 (MQ=255)
              gTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAgggg                                                                                                                                        >  1:540137/1‑150 (MQ=255)
               ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGc                                                                                                                                       >  2:159522/1‑150 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGGTGCAGGCGCTGGTGTCGCAGGTGGGTGGATGCCGGATTCGTCTCGGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  2:2629440/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     >  2:1047573/1‑150 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  1:153039/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  1:1047573/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGACTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  1:2372038/150‑1 (MQ=255)
                  ggTGGCTCAGGCCCACGGCGTGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                                    <  1:159522/150‑1 (MQ=255)
                  ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                                    <  1:2655600/150‑1 (MQ=255)
                    tGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGtt                                                                                                                                  <  2:540137/150‑1 (MQ=255)
                       cTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTAc                                                                                                                               <  2:496098/150‑1 (MQ=37)
                                    cggGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGa                                                                                                                         >  1:347856/1‑143 (MQ=255)
                                    cggGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGa                                                                                                                         <  2:347856/143‑1 (MQ=255)
                                           cTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                     <  1:117844/124‑1 (MQ=25)
                                              cAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                     <  1:2400666/121‑1 (MQ=25)
                                                   gcTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACAt                                                                                                   <  2:2663657/150‑1 (MQ=255)
                                                     tGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  2:2662124/114‑1 (MQ=255)
                                                     tGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     >  1:2662124/1‑114 (MQ=255)
                                                          tCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                               <  1:3039927/99‑1 (MQ=255)
                                                          tCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                               >  2:3039927/1‑99 (MQ=255)
                                                                       gATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAg                                                                                          <  1:1760533/139‑1 (MQ=255)
                                                                       gATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAg                                                                                          >  2:1760533/1‑139 (MQ=255)
                                                                         tGCCGGATTCGTGTCAGTCTGCCTTTGAATCCCGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     <  2:2719598/94‑1 (MQ=255)
                                                                         tGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                     >  1:2719598/1‑94 (MQ=255)
                                                                         tGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                             <  1:313058/102‑1 (MQ=255)
                                                                         tGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                             >  2:313058/1‑102 (MQ=255)
                                                                          gCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                   >  2:1897120/1‑95 (MQ=255)
                                                                          gCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                   <  1:1897120/95‑1 (MQ=255)
                                                                               aTTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTc                                                                                                                                                   <  1:522068/74‑1 (MQ=255)
                                                                               aTTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTc                                                                                                                                                   >  2:522068/1‑74 (MQ=255)
                                                                                    tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                               <  1:2598141/89‑1 (MQ=255)
                                                                                    tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                               >  2:2598141/1‑89 (MQ=255)
                                                                                     gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGaa                                                                                                                        >  1:2238220/1‑95 (MQ=255)
                                                                                     gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGaa                                                                                                                        <  2:2238220/95‑1 (MQ=255)
                                                                                         gTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                      <  1:1352906/93‑1 (MQ=255)
                                                                                         gTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                      >  2:1352906/1‑93 (MQ=255)
                                                                                             gCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCgg                                                         <  1:1163715/150‑1 (MQ=255)
                                                                                                tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAg                                                      >  2:1185791/1‑150 (MQ=255)
                                                                                                      tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTAc                                                                                                     <  2:2491127/97‑1 (MQ=255)
                                                                                                      tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGTGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  >  2:1628679/1‑148 (MQ=255)
                                                                                                      tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  <  1:1628679/148‑1 (MQ=255)
                                                                                                      tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACACGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGtttcc                                                                                                     >  1:2491127/1‑95 (MQ=255)
                                                                                                      tCACGCATAGCCTGTGCCGGCCTCTTCGCCGGCAAACCCCCTCCTACGGTCGGCTGTAGGGGCTGGTGTACCCGCGACGAGGCCAGTACAGGTTTACGTCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCCTCAGGTTGCGGCAGCTgggg                                                >  2:1389383/1‑150 (MQ=255)
                                                                                                                catgtgCCGGCCTCTTCGCGGGTGAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCg                                      >  2:2348219/3‑150 (MQ=255)
                                                                                                                           ctctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgc                                                                                    >  2:2214745/1‑93 (MQ=255)
                                                                                                                           ctctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgc                                                                                    <  1:2214745/93‑1 (MQ=255)
                                                                                                                            tctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAc                          >  2:1207055/1‑150 (MQ=255)
                                                                                                                            tctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAc                          >  2:814757/1‑150 (MQ=255)
                                                                                                                                 gcgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                     >  2:679171/1‑150 (MQ=255)
                                                                                                                                  cgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATAAGCGTGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc                    >  1:371300/1‑150 (MQ=255)
                                                                                                                                   gggTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc                   >  1:227673/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt                 >  1:1077820/1‑150 (MQ=255)
                                                                                                                                     gTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCg                                                          <  2:722006/109‑1 (MQ=18)
                                                                                                                                     gTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCg                                                          >  1:722006/1‑109 (MQ=21)
                                                                                                                                      tAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc                <  1:2348219/150‑1 (MQ=255)
                                                                                                                                      tAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc                >  2:228006/1‑150 (MQ=255)
                                                                                                                                      tAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc                >  2:228616/1‑150 (MQ=255)
                                                                                                                                         aCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                                   >  2:394454/1‑128 (MQ=255)
                                                                                                                                         aCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                                   <  1:394454/128‑1 (MQ=255)
                                                                                                                                         aCCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgggg                                                <  2:815445/115‑1 (MQ=25)
                                                                                                                                         aCCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgggg                                                >  1:815445/1‑115 (MQ=25)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgc                                                                                    >  1:2603934/1‑78 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgc                                                                                    <  2:2603934/78‑1 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg            >  2:604513/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg            >  2:1859942/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAgg                                                                                         <  2:2980565/72‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAgg                                                                                         >  1:2980565/1‑72 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                                 <  1:2934373/112‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                                 >  2:2934373/1‑112 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt           >  1:2701213/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt           >  1:611545/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt           >  1:686249/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt           >  1:2455385/1‑150 (MQ=255)
                                                                                                                                           ccGCCCCTACAGCCGACCGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt           >  2:209883/1‑150 (MQ=38)
                                                                                                                                             gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgggg                                                >  2:1808302/1‑111 (MQ=255)
                                                                                                                                             gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgggg                                                <  1:1808302/111‑1 (MQ=255)
                                                                                                                                             gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg         >  1:1551278/1‑150 (MQ=255)
                                                                                                                                               tCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGtt       <  1:1185791/150‑1 (MQ=255)
                                                                                                                                                ccTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTg      >  1:785718/1‑150 (MQ=255)
                                                                                                                                                 cTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGGCGACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCCGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc     >  1:2091616/1‑150 (MQ=255)
                                                                                                                                                  tACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  <  1:3068295/104‑1 (MQ=255)
                                                                                                                                                  tACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  >  2:3068295/1‑104 (MQ=255)
                                                                                                                                                  tACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgcg    <  1:604513/150‑1 (MQ=255)
                                                                                                                                                    cGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                       >  1:117520/1‑113 (MQ=255)
                                                                                                                                                    cGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                       <  2:117520/113‑1 (MQ=255)
                                                                                                                                                    cGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGcc  >  2:1720730/1‑150 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  <  1:2040465/101‑1 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                  >  2:2040465/1‑101 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCC  >  NC_002947/5984349‑5984646

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