Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,716,337 A→C 5.7% N523T (AAT→ACT)  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,3370AC 5.7% 223.3 / 3.9 70N523T (AAT→ACT) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base C (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

GAAAGCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGAC  >  NC_002947/4716195‑4716485
                                                                                                                                               |                                                                                                                                                    
gAAAGCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGt                                                                                                                                                >  7:343107/1‑149 (MQ=255)
         gAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAg                                                                                                                                       >  5:18368/1‑149 (MQ=255)
         gAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAg                                                                                                                                       >  7:133594/1‑149 (MQ=255)
          aaGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAgt                                                                                                                                      <  3:108455/149‑1 (MQ=255)
             ggcTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCATCCCTGCAC‑CCGTCCCGTTGACGCCGGCAAAGAATGTACGAACGACGGTGCCGGACAATGATGAGTTACTGATGTTGGAGGACAATGCGGTCAACCGgtgtgt                                                                                                                                    <  3:158949/146‑1 (MQ=255)
                  ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTgg                                                                                                                               >  5:22518/1‑148 (MQ=255)
                  ggTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTgg                                                                                                                               >  5:28041/1‑148 (MQ=255)
                    ctCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAg                                                                                                                            <  8:133594/148‑1 (MQ=255)
                        aCCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCCGCCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                         >  1:263826/1‑149 (MQ=255)
                        aCCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGATTACGCCGCCAACGACTGTACCAACGACGGTGCCGGACAATGATCAGCTACTGCTGGTGGAGGACAATCCGGTCAACAGGAGTGTCTTGGAAGCCAt                                                                                                                        >  7:311911/1‑149 (MQ=255)
                               cccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTgcg                                                                                                                  >  7:10473/1‑148 (MQ=255)
                                ccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTc                                                                                                                                   >  2:228976/1‑130 (MQ=255)
                                ccTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTc                                                                                                                                   <  1:228976/130‑1 (MQ=255)
                                       aTGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcc                                                                                                                                  >  8:342328/1‑123 (MQ=255)
                                       aTGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACACTCCGGTTGACCGGATTGTcc                                                                                                                                  <  7:342328/124‑2 (MQ=255)
                                       aTGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTggg                                                                                                         >  1:65204/1‑149 (MQ=255)
                                       aTGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTggg                                                                                                         <  8:350741/149‑1 (MQ=255)
                                        tGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                        <  2:263826/149‑1 (MQ=255)
                                        tGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                        <  7:2804/149‑1 (MQ=255)
                                         ggcGTTGGCCCTTGGCAACCGTGCGC‑CCGTGCCGTTGACGGCGCCAACGACTGTACGATGGACGGTGCCGGACAATGTTCAGGTACTGCTGGTGGATGTCAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTgggg                                                                                                        <  8:311911/146‑1 (MQ=255)
                                           cGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTg                                                                                                     <  8:138918/149‑1 (MQ=255)
                                            gTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTgg                                                                                                    <  8:160006/149‑1 (MQ=255)
                                               ggCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCa                                                                                                  <  8:343107/148‑1 (MQ=255)
                                                gCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAg                                                                                                 >  2:180383/1‑148 (MQ=255)
                                                  ccTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGc                                                                                                                           >  4:229708/1‑120 (MQ=255)
                                                  ccTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                              >  8:131706/1‑149 (MQ=255)
                                                  ccTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTAt                                                                                              >  4:276779/1‑149 (MQ=255)
                                                        cAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACgg                                                                                        <  4:58349/149‑1 (MQ=255)
                                                            ccTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                     >  3:185436/1‑148 (MQ=255)
                                                            ccTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCag                                                                                     >  4:312414/1‑148 (MQ=255)
                                                                ccc‑ccGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGAt                                                                                <  7:137866/149‑1 (MQ=255)
                                                                 cc‑ccGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATg                                                                               <  7:314260/149‑1 (MQ=255)
                                                                    ccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                         >  8:58622/1‑105 (MQ=255)
                                                                    ccgtcccgtTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCa                                                                                                                         <  7:58622/105‑1 (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                                                                  >  1:370293/1‑149 (MQ=255)
                                                                               aCGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGt                                                                  <  6:18368/149‑1 (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                            <  8:10473/149‑1 (MQ=255)
                                                                                                                     tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                            <  7:131706/149‑1 (MQ=255)
                                                                                                                     tGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGa                            <  2:370293/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            >  5:74901/1‑149 (MQ=255)
                                                                                                                                     gtggAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCgcctgcc            >  3:190483/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)
                                                                                                                                               |                                                                                                                                                    
GAAAGCCGTGAAGGCCTGGGTTCGACCTTCACCCTGGAGATGCCGTTGGCCCTTGCCAACCCTGCCC‑CCGTCCCGTTGACGCCGCCAACGACTGTACGAACGACGGTGCCGGACAATGATCAGGTACTGCTGGTGGAGGACAATCCGGTCAACCGGAGTGTCGTGGAAGCCATGCTGCGCAGCCTGGGGCTGGCAGTATGCACGGCAGAAGATGGCCTTGAAGCCGTCGAGCTGGTCGGCCGGCAGCGCTTTGCTGCCGTGCTGATGGATTGCCGCCTGCCTCACCTGGAC  >  NC_002947/4716195‑4716485

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