Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,716,352 G→T 5.5% S528I (AGT→ATT)  PP_4173 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,716,3520GT 5.5% 234.0 / 3.4 72S528I (AGT→ATT) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (30/38);  new base T (2/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

GAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACG  >  NC_002947/4716204‑4716486
                                                                                                                                                     |                                                                                                                                      
gAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAg                                                                                                                                        >  7:133594/1‑149 (MQ=255)
gAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAg                                                                                                                                        >  5:18368/1‑149 (MQ=255)
 aaGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAgt                                                                                                                                       <  3:108455/149‑1 (MQ=255)
    ggcTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCATCCCTGCACCCGTCCCGTTGACGCCGGCAAAGAATGTACGAACGACGGTGCCGGACAATGATGAGTTACTGATGTTGGAGGACAATGCGGTCAACCGgtgtgt                                                                                                                                     <  3:158949/146‑1 (MQ=255)
         ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTgg                                                                                                                                >  5:28041/1‑148 (MQ=255)
         ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTgg                                                                                                                                >  5:22518/1‑148 (MQ=255)
           ctCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAg                                                                                                                             <  8:133594/148‑1 (MQ=255)
               aCCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCGCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                          >  1:263826/1‑149 (MQ=255)
               aCCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGATTACGCCGCCAACGACTGTACCAACGACGGTGCCGGACAATGATCAGCTACTGCTGGTGGAGGACAATCCGGTCAACAGGAGTGTCTTGGAAGCCAt                                                                                                                         >  7:311911/1‑149 (MQ=255)
                      cccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTgcg                                                                                                                   >  7:10473/1‑148 (MQ=255)
                       ccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTc                                                                                                                                    <  1:228976/130‑1 (MQ=255)
                       ccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTc                                                                                                                                    >  2:228976/1‑130 (MQ=255)
                              aTGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcc                                                                                                                                   >  8:342328/1‑123 (MQ=255)
                              aTGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcc                                                                                                                                   <  7:342328/124‑2 (MQ=255)
                              aTGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTggg                                                                                                          >  1:65204/1‑149 (MQ=255)
                              aTGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTggg                                                                                                          <  8:350741/149‑1 (MQ=255)
                               tGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                         <  7:2804/149‑1 (MQ=255)
                               tGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                         <  2:263826/149‑1 (MQ=255)
                                ggcGTTGGCCCTTGGCAACCGTGCGCCCGTGCCGTTGACGGCGCCAACGACTGTACGATGGACGGTGCCGGACAATGTTCAGGTACTGCTGGTGGATGTCAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                         <  8:311911/146‑1 (MQ=255)
                                  cGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTg                                                                                                      <  8:138918/149‑1 (MQ=255)
                                   gTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTgg                                                                                                     <  8:160006/149‑1 (MQ=255)
                                      ggCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCa                                                                                                   <  8:343107/148‑1 (MQ=255)
                                       gCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAg                                                                                                  >  2:180383/1‑148 (MQ=255)
                                         ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGc                                                                                                                            >  4:229708/1‑120 (MQ=255)
                                         ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                               >  8:131706/1‑149 (MQ=255)
                                         ccTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                               >  4:276779/1‑149 (MQ=255)
                                               cAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACgg                                                                                         <  4:58349/149‑1 (MQ=255)
                                                   ccTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                      >  3:185436/1‑148 (MQ=255)
                                                   ccTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                      >  4:312414/1‑148 (MQ=255)
                                                       cccccGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGAt                                                                                 <  7:137866/149‑1 (MQ=255)
                                                        ccccGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATg                                                                                <  7:314260/149‑1 (MQ=255)
                                                           ccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                          <  7:58622/105‑1 (MQ=255)
                                                           ccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                          >  8:58622/1‑105 (MQ=255)
                                                           ccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGACGTGGCAGCCATGCTGCGCAGCCTGCGGCTGGCAGTATGCACGGCAGAAGATgcc                                                                              >  8:267891/1‑147 (MQ=255)
                                                             gtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCt                                                                            <  1:180383/149‑1 (MQ=255)
                                                              tcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCtt                                                                           <  2:301178/149‑1 (MQ=255)
                                                                ccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTg                                                                          >  1:292513/1‑148 (MQ=255)
                                                                      aCGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGt                                                                   <  3:189006/149‑1 (MQ=255)
                                                                      aCGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGt                                                                   <  6:18368/149‑1 (MQ=255)
                                                                      aCGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGt                                                                   >  1:370293/1‑149 (MQ=255)
                                                                         ccgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGa                                                                <  2:292513/149‑1 (MQ=255)
                                                                         ccgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCATCGa                                                                <  3:276779/149‑1 (MQ=255)
                                                                          cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGc                                                                                        <  1:128118/124‑1 (MQ=255)
                                                                          cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGc                                                                                        >  2:128118/1‑124 (MQ=255)
                                                                          cgccAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAg                                                               >  2:246082/1‑149 (MQ=255)
                                                                              aaCGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTgg                                                           <  5:301754/149‑1 (MQ=255)
                                                                                      tACGAACGACGGTGCCGGACAATGAGCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTcggccgg                                                   <  7:267891/149‑1 (MQ=255)
                                                                                         gaacgaCTGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTc                                                                                                                                    >  2:105495/1‑65 (MQ=255)
                                                                                         gaacgaCGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTc                                                                                                                                    <  1:105495/65‑1 (MQ=255)
                                                                                                ggTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTGAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTtg                                         <  6:28041/149‑1 (MQ=255)
                                                                                                ggTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGGCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGATGCCGTCGAGCTGGTCGGCCGGCAGCGCTTtg                                         <  6:22518/149‑1 (MQ=255)
                                                                                                     cGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGcc                                    <  1:246082/149‑1 (MQ=255)
                                                                                                         cAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGc                                >  5:108255/1‑149 (MQ=255)
                                                                                                          aaTGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTcggcc                                                     >  8:197541/1‑127 (MQ=255)
                                                                                                          aaTGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTcggcc                                                     <  7:197541/127‑1 (MQ=255)
                                                                                                            tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCa                                                 <  4:209497/129‑1 (MQ=255)
                                                                                                            tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCa                                                 >  3:209497/1‑129 (MQ=255)
                                                                                                            tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                             <  7:131706/149‑1 (MQ=255)
                                                                                                            tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                             <  2:370293/149‑1 (MQ=255)
                                                                                                            tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                             <  8:10473/149‑1 (MQ=255)
                                                                                                               tCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGAAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGGGCGGATg                           >  5:214936/1‑148 (MQ=255)
                                                                                                                      ctgctgGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGTCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTgcc                    <  2:65204/148‑1 (MQ=255)
                                                                                                                      ctgctgGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTgcc                    <  3:56219/148‑1 (MQ=255)
                                                                                                                       tgctgGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTgccgc                  >  2:80579/1‑149 (MQ=255)
                                                                                                                           ggtggAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCAGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCAGGCAGCGCTTTGATGCCGTGCTGATGGATTGCCgcctgc              >  2:239571/1‑149 (MQ=255)
                                                                                                                            gtggAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCgcctgcc             >  3:190483/1‑149 (MQ=255)
                                                                                                                            gtggAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCgcctgcc             >  5:74901/1‑149 (MQ=255)
                                                                                                                                  gaCAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCt       <  1:80579/149‑1 (MQ=255)
                                                                                                                                  gaCAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCt       <  3:312414/149‑1 (MQ=255)
                                                                                                                                     aaTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGc                                                              <  2:136468/91‑1 (MQ=255)
                                                                                                                                      aTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                             <  3:182229/123‑1 (MQ=255)
                                                                                                                                      aTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                             >  4:182229/1‑123 (MQ=255)
                                                                                                                                      aTCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGAc   <  4:185436/149‑1 (MQ=255)
                                                                                                                                       tCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACg  >  3:189805/1‑149 (MQ=255)
                                                                                                                                                   gAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTtgctg                                      >  3:182226/1‑101 (MQ=255)
                                                                                                                                                   gAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTtgctg                                      <  4:182226/101‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
GAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGACG  >  NC_002947/4716204‑4716486

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