Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,989,447 G→A 26.6% 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,989,4470GA26.6% 320.8 / 142.0 211intergenic (+38/+51)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (66/85);  new base A (22/32);  total (94/117)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

Only 100 of 220 total aligned reads displayed.

CTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCA  >  NC_002947/5989298‑5989596
                                                                                                                                                     |                                                                                                                                                     
cTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACg                                                                                                                                                       <  1:5495162/150‑1 (MQ=255)
 ttGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgg                                                                                                                                                      >  2:414238/1‑150 (MQ=255)
      gATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGc                                                                                                                                                 <  1:5134954/150‑1 (MQ=255)
       aTCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTGAAACCCCTCCAACGggcggcg                                                                                                                                                >  1:5022145/1‑149 (MQ=255)
       aTCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                                <  2:6881511/150‑1 (MQ=255)
           gTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAg                                                                                                                                            <  2:6945566/150‑1 (MQ=39)
            tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAgg                                                                                                                                           <  2:4982876/150‑1 (MQ=255)
               ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                        <  1:935778/150‑1 (MQ=38)
               ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCTCCGACGGGCGGCGGGGGGGGc                                                                                                                                        >  2:1933622/1‑150 (MQ=39)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                      <  2:5739155/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                      <  2:1544546/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                      <  1:1271858/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                      <  1:80129/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                      <  2:6154199/150‑1 (MQ=255)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  1:4969776/150‑1 (MQ=38)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  2:565953/150‑1 (MQ=38)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  2:165725/150‑1 (MQ=38)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  1:4928948/150‑1 (MQ=38)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  2:2939618/150‑1 (MQ=38)
                 cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  1:888370/150‑1 (MQ=38)
                  ggTGGCTCAGGCTCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACGGTCGGCGGTAGGGGCggg                                                                                                                                     >  2:7050461/1‑150 (MQ=255)
                  ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                                     <  2:2434826/150‑1 (MQ=255)
                  ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                                     <  1:3489852/150‑1 (MQ=255)
                  ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                                     <  2:3003158/150‑1 (MQ=38)
                   gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                    <  1:4151947/150‑1 (MQ=255)
                   gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                    <  1:7259745/150‑1 (MQ=255)
                   gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                                    <  1:6410843/150‑1 (MQ=255)
                   gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGt                                                                                                                                    <  1:57958/150‑1 (MQ=38)
                    tGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCGCCGACGGGCGGCGGGTgggggggggt                                                                                                                                   >  1:3262680/1‑148 (MQ=255)
                      gCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTa                                                                                                                                 >  2:4844582/1‑150 (MQ=255)
                       cTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                                <  2:3581707/150‑1 (MQ=255)
                        tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCTCGCATAGCCTGTGCCGGCCTCTTCGCGGGTAATCCCGCTCCTACGGTCGGGTGTAGGGGGGGGttttac                                                                                                                               >  2:5442906/1‑147 (MQ=39)
                        tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCGACGGTGGGCGGTAGGGGCGGGTTTAcc                                                                                                                               >  1:5891952/1‑150 (MQ=255)
                        tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCCCCCCCACAGCCGACCGGAGGAGCGGGTTTAcc                                                                                                                               >  1:2177865/1‑150 (MQ=25)
                         cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                              <  2:6025299/150‑1 (MQ=255)
                         cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTAccc                                                                                                                              <  1:5653096/150‑1 (MQ=37)
                         cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTAccc                                                                                                                              <  1:3121626/150‑1 (MQ=25)
                             cccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGa                                                                                                                          >  2:795423/1‑150 (MQ=255)
                                aCGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCAACGGTCGGCTGTAGGGGCGGGTTTTCCCGCGAaga                                                                                                                       >  1:1804248/1‑150 (MQ=255)
                                 cggcggGCAGCTGCAGGCGCTGGTGTCGCAGTGGGGTGGTTGCCGGATTGGTGTCAGTCTGCCTTTTAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGTCGGGGGGGGTTTACCCGCGAagag                                                                                                                      <  2:5020709/150‑1 (MQ=38)
                                  ggcggGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGTTTTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAgg                                                                                                                     <  2:4179050/150‑1 (MQ=255)
                                   gcggGCAGCTGCCGGGGCGGGTGTCGCTGTTGGGTGGATGCCGGATTCGTGTCATTCGGCCTTTGAATCACGCATAGCGTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTGGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGc                                                                                                                    <  2:3866410/150‑1 (MQ=37)
                                     gggCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCAACGGTCGGCTGTAGGGGCGGGTTTAACCCCGCAGAGGCCa                                                                                                                  >  1:3589730/1‑150 (MQ=255)
                                           ctgctgGCGGTGGTGTCGCTGTGGGGGGGATGCGGGATTCGTGTCAGTCGGCCTTTGAATCACGCATCGCCTGTGCCGGCCTCTTGGCGGGTAACCCCGCTCCTCCGGTGGGCTGTGGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAg                                                                                                            <  1:6732880/150‑1 (MQ=25)
                                           cggccggCGCTGGTGTCGCTGTTGGGTGGCTGCGGGTTTCGTTTCAGTCTGCCTTTTAATCCCGCATAGCCTGTGCCGGCCACTGCGGGGGTAAACCCGCTCCTCCGGTCGGCTGTCGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAg                                                                                                            <  2:2076422/145‑1 (MQ=255)
                                           cTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACGGTCGGCTGTAGGGGCGGGGTTACCCGCGACGAGGCCAGTACAg                                                                                                            >  2:6463619/1‑150 (MQ=255)
                                            tGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      >  1:7356314/1‑123 (MQ=37)
                                            tGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                      <  2:7356314/123‑1 (MQ=37)
                                                     tGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCa                                                                                                  >  2:1078155/1‑150 (MQ=255)
                                                        tgtCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAgcg                                                                                               <  2:7298160/150‑1 (MQ=255)
                                                                     tGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                                     <  2:5872363/99‑1 (MQ=21)
                                                                             gttttCGTGTCCGTCGGCATTTGACTCCCGCAGCGCCTGGGCCGGCCTCTGCGGGGTTAACCCCGCCCCACCGGTCGGCTGTGGGGGGGGTTTCACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCt                                                                          <  1:1020083/147‑1 (MQ=255)
                                                                             ggATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACGGTCGGCGGGTGGGGCGGGTTTACCCGCGGCGAGGCCACTACAGGTTTACATCAGCCCGAGGCCAGCGCTGCGGAGcg                                                                          >  2:2712766/1‑149 (MQ=38)
                                                                              gATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCCCGAACAGGCCAGTACAGGTTTAAAACAGCGCGAGGCCAGCGCAGCGGAGcgt                                                                         >  2:1435544/1‑148 (MQ=38)
                                                                                  cGTGTCAGTCTGCCTTTGAATCACGCCTCGCCTGTGCCGGCCTCTCCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtca                                                                     <  2:1964824/150‑1 (MQ=38)
                                                                                   gtgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcat                                                                    >  2:1749960/1‑150 (MQ=255)
                                                                                     gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                                  <  2:1365214/150‑1 (MQ=255)
                                                                                     gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGGAAACCCCCGCCGACGGTCGGCGGGAGGGGCGGGTTTACCCCCGAAGAGGCCACGACAGGTTTACAACAGCCCCCGGCCCGCCCGGCGTAGCGGttaaaca                                                                  >  1:4832254/1‑145 (MQ=255)
                                                                                      tcagtcagCCTTTGAATCACGCATAGCCTGTGCCGGCCTCGTCGCGGGTAAACCCGCTCCTACGCTCGGCGGGAGGGGCGGGGTTACCCCCCACCAGGCCACTACCGGGTTACATCAGCGCGCGGCCACCGCTGCGGAGCGGTTCATCAg                                                                 >  2:2826648/1‑150 (MQ=25)
                                                                                          tCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACGGTCGGGTGTAGGGGCGGGTTTAACCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTGATAAGGTTg                                                             >  2:1629928/1‑150 (MQ=38)
                                                                                           cTGCCTTTTAATCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTTGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGc                                                            <  2:3262680/150‑1 (MQ=255)
                                                                                             gCCTTTTAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCgg                                                          <  1:6059999/150‑1 (MQ=255)
                                                                                             gCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAGGTCGGGTGGAGGGGGGGGTTTTGCCGCGGAGTGGGCAGTTCAGGTTTTGATGAGCGGGAGGGCAGCGGTGGGTAGGTGTTGATGGGGTTgggg                                                          >  2:1407420/1‑150 (MQ=16)
                                                                                               cTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTa                                                        >  2:6781059/1‑150 (MQ=255)
                                                                                                tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGGGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAg                                                       >  2:2364748/1‑150 (MQ=255)
                                                                                                    aaacACGCACCGCCTGGGCCGGCCTCTTCGGGTTTAACCCGGCTCATCGGGCGGGCTGTCGGGGCGGTTTAACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                   <  1:1407420/147‑1 (MQ=255)
                                                                                                    aaTCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                                   <  2:2897538/150‑1 (MQ=255)
                                                                                                        aCGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCGGGAGGGGCGGGTTTTCCCGCGAAGAGGCCCGTAAAGGTTTACATCAGCGCGAGGCCCGCGCTGCGGTGCTGTTGCTCCGGTTGCGGGAGTTGGGGAt                                               >  2:1597869/1‑150 (MQ=37)
                                                                                                                   gtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTTGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                                    <  1:1729426/150‑1 (MQ=255)
                                                                                                                   gtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                                    <  2:3316216/150‑1 (MQ=255)
                                                                                                                       cGGCCTCTTCGCGGGTAAACCCGCTACTAAGGTAGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTATTTGGGGATACGCTGCGAcagcag                                >  2:538666/1‑150 (MQ=255)
                                                                                                                         cccactTCGCGGGAAAACCCGAACCTACGGACGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGt                              <  1:5442906/146‑1 (MQ=255)
                                                                                                                           ctctTCGCGGCTAACCCCGCCCCTCCGGTCGGCTTTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTa                            <  1:414238/150‑1 (MQ=255)
                                                                                                                            tctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAc                           >  2:5565455/1‑150 (MQ=255)
                                                                                                                             cgtCGCGGGTCAACCCGCTCCGCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAcc                          <  2:6980756/148‑1 (MQ=255)
                                                                                                                                 gcgGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                      <  1:2986847/150‑1 (MQ=255)
                                                                                                                                   gtttAACCCCGCGCCTTGGGGCGGGTGTTGGGGGGGGTTCACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc                    <  1:1362935/147‑1 (MQ=255)
                                                                                                                                   gggTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc                    >  1:525338/1‑150 (MQ=37)
                                                                                                                                    ggTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCt                   >  1:4125119/1‑150 (MQ=37)
                                                                                                                                       aaaCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCg                >  2:7422553/1‑150 (MQ=255)
                                                                                                                                        accccGCTCCTCCGGGGGGCTGTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGa               <  1:3995623/148‑1 (MQ=255)
                                                                                                                                         aCCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt              <  2:7209413/150‑1 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtc                                                                      >  2:3122161/1‑93 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg             >  1:7235557/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg             >  2:4061847/1‑150 (MQ=255)
                                                                                                                                          cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                                  <  1:2036451/113‑1 (MQ=25)
                                                                                                                                          cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                                  >  2:2036189/1‑113 (MQ=25)
                                                                                                                                           ccGCTCCTCCGGTCGGCTGTCGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt            <  1:804302/150‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt            >  1:1125305/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt            >  1:6942390/1‑150 (MQ=255)
                                                                                                                                           ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt              >  1:70016/1‑148 (MQ=38)
                                                                                                                                             gccccTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg          >  2:3812247/4‑150 (MQ=255)
                                                                                                                                              ctcctcCGGTGGGCTGTGGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGt         <  1:1953950/150‑1 (MQ=39)
                                                                                                                                               tCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGtt        >  2:4757363/1‑150 (MQ=255)
                                                                                                                                                 cTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc      >  1:6290620/1‑150 (MQ=255)
                                                                                                                                                 cTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc      >  2:2648699/1‑150 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTa                                                        >  2:6435977/1‑96 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTa                                                        <  1:6435977/96‑1 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCa  >  2:5173929/1‑150 (MQ=255)
                                                                                                                                                     ggTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCa  >  2:1343261/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGCGCCA  >  NC_002947/5989298‑5989596

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