Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,188 G→A 9.2% E368K (GAG→AAG)  PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,1880GA 9.2% 210.4 / 5.8 75E368K (GAG→AAG) PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (28/40);  new base A (5/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGA  >  NC_002947/1622042‑1622333
                                                                                                                                                  |                                                                                                                                                 
aaGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAg                                                                                                                                                 >  1:196052/1‑149 (MQ=255)
      gTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttgtt                                                                                                                                            >  5:184848/1‑148 (MQ=255)
       tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAATGAGTTGTTGa                                                                                                                                          <  6:182955/149‑1 (MQ=255)
       tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                          <  7:160660/149‑1 (MQ=255)
        gcgcAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAg                                                                                                                                         >  8:241716/1‑149 (MQ=255)
           cAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCa                                                                                                                                       >  1:186605/1‑148 (MQ=255)
           cAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCa                                                                                                                                       <  2:186605/148‑1 (MQ=255)
            aGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCCCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAAc                                                                                                                                     <  7:267728/149‑1 (MQ=255)
             gCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAAcc                                                                                                                                    >  8:252561/1‑149 (MQ=255)
              ggTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCaatgt                                                                                                                                   <  7:39587/148‑4 (MQ=255)
                 cTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggt                                                                                                                                >  2:85681/1‑149 (MQ=255)
                  tGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCACCCtggtg                                                                                                                               >  8:139687/1‑149 (MQ=255)
                    gACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggtg                                                                                                                               <  8:295847/147‑1 (MQ=255)
                     aCTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGa                                                                                                                             >  2:315076/1‑148 (MQ=255)
                       ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtg                                                                                                                                  <  6:31659/141‑1 (MQ=255)
                       ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtg                                                                                                                                  >  5:31659/1‑141 (MQ=255)
                       ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACaa                                                                                                                          >  1:49525/1‑149 (MQ=255)
                       ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAGGTGGACCCCAAATAAATAAAAAAAAAATTGAGCAACCTGGTGGACaa                                                                                                                          >  6:278690/1‑149 (MQ=255)
                           gTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGcc                                                                                                                      >  3:310874/1‑149 (MQ=255)
                              aGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAGGTGGACCCCAAATAAATAAAAAAAAAATTGAGCAACCTGGTGGACAATGCCttt                                                                                                                   >  4:293158/1‑148 (MQ=255)
                               gTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTgg                                                                                                                  >  2:42692/1‑149 (MQ=255)
                               gTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTgg                                                                                                                  >  3:162204/1‑149 (MQ=255)
                               gTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTgg                                                                                                                  <  2:49525/149‑1 (MQ=255)
                                ttGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgc                                                                                                                 <  7:247206/149‑1 (MQ=255)
                                  ggCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGgcg                                                                                                                >  6:212930/1‑148 (MQ=255)
                                     ccGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATa                                                                                                            <  1:315076/149‑1 (MQ=255)
                                            cTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGc                                                                                                     <  7:99888/149‑1 (MQ=255)
                                            cTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGGGGCGCTGGCGCTGGAGGCCGTGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGc                                                                                                     <  3:204831/149‑1 (MQ=255)
                                              gggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCgg                                                                                                   >  2:262997/1‑149 (MQ=255)
                                                gCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGt                                                                                                 <  2:196052/149‑1 (MQ=255)
                                                 catggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                >  2:136501/1‑149 (MQ=255)
                                                  atggccatggcTCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTgg                                                                                               <  4:162204/149‑1 (MQ=255)
                                                    ggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGc                                                                                              >  5:216273/1‑148 (MQ=255)
                                                     gccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCaa                                                                                            <  5:80853/149‑1 (MQ=255)
                                                      ccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAt                                                                                           <  1:85681/149‑1 (MQ=255)
                                                        atggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAtgt                                                                                         <  4:269991/149‑1 (MQ=255)
                                                        atggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGGGGCAAtgt                                                                                         >  7:984/1‑149 (MQ=255)
                                                           gcgcCTTTGGCGCACAAGCGGGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGa                                                                                       <  3:293158/148‑1 (MQ=255)
                                                            cgcCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGTCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGAtt                                                                                     <  1:262997/149‑1 (MQ=255)
                                                            cgcCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGAtt                                                                                     <  6:216273/149‑1 (MQ=255)
                                                             gcCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTAGATCAAGGTGGACCCGACATTAATAAAAAAATTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCg                                                                                                    >  6:284256/1‑133 (MQ=255)
                                                             gcCTTTCGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCg                                                                                                    <  5:284256/133‑1 (MQ=255)
                                                                tttGGCGCACAAGGGGGGGGTGGCGCTGGCGGTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGc                                                                                 <  4:310874/149‑1 (MQ=255)
                                                                  tGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGc                                                                                 <  8:9655/147‑1 (MQ=255)
                                                                  tGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCg                                                                                <  1:42692/148‑1 (MQ=255)
                                                                  tGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCg                                                                                <  6:184848/148‑1 (MQ=255)
                                                                   ggCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCgg                                                                               <  8:111174/148‑1 (MQ=255)
                                                                   ggCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCgg                                                                               <  4:87763/148‑1 (MQ=255)
                                                                   ggCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCgg                                                                               <  2:176374/148‑1 (MQ=255)
                                                                      gcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTaa                                                                           >  7:258469/1‑149 (MQ=255)
                                                                      gcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTaa                                                                           >  5:321280/1‑149 (MQ=255)
                                                                            gcgcGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggt                                                                                                                                <  3:112272/90‑1 (MQ=255)
                                                                            gcgcGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTCGATGAAGGAGGAGCCGACATTACTCAACAACTAGTTTAGCAACAtggt                                                                                                                                >  4:112272/1‑90 (MQ=255)
                                                                                        gctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtg                                                          >  5:275441/1‑148 (MQ=255)
                                                                                            gcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                     <  4:259257/149‑1 (MQ=255)
                                                                                             cgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGGGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                    <  8:984/149‑1 (MQ=255)
                                                                                                tggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtggaa                                                 <  6:275441/149‑1 (MQ=255)
                                                                                                tggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtggaa                                                 <  8:258469/149‑1 (MQ=255)
                                                                                                  gcggCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAga                                               <  5:278690/147‑1 (MQ=255)
                                                                                                     gCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                    >  6:206288/1‑141 (MQ=255)
                                                                                                     gCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                    <  5:206288/141‑1 (MQ=255)
                                                                                                      ccGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAgacgac                                           <  5:212930/149‑1 (MQ=255)
                                                                                                             gcCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAg                                    <  7:241716/149‑1 (MQ=255)
                                                                                                                   gtgGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGACTTGGCGCATACGCCAGCGGGTGGCAATTTGATTTTGCGGGTAATGGCGCAGGCAGAGCTGGAAGTGGAAGACGACGGGCCAGGTATTc                              >  5:205188/1‑149 (MQ=255)
                                                                                                                     ggTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTccc                            >  2:189558/1‑149 (MQ=255)
                                                                                                                      gTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAtgtg                                                                                        >  6:114617/1‑88 (MQ=255)
                                                                                                                      gTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAtgtg                                                                                        <  5:114617/88‑1 (MQ=255)
                                                                                                                      gTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAgacg                                             >  3:96423/1‑131 (MQ=255)
                                                                                                                        tGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTg                         >  5:145608/1‑149 (MQ=255)
                                                                                                                        tGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAAGTGATTTTGCGGGTAATGGCTCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCCGGTATTCACGTg                         >  6:188821/1‑149 (MQ=255)
                                                                                                                          acgggggAGCCGACATTGCTCAACGAGTTGTTCAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAAGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTgg                        <  7:139687/146‑1 (MQ=255)
                                                                                                                                     gACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGTCCGAGCGCGAAc            <  6:321280/149‑1 (MQ=255)
                                                                                                                                         ttGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgt        <  1:189558/149‑1 (MQ=255)
                                                                                                                                             tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc    >  5:287937/1‑149 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGa  <  7:252561/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AAGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGA  >  NC_002947/1622042‑1622333

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: