Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,336,179 G→A 6.7% R870H (CGC→CAC)  PP_4695 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,336,1790GA 6.7% 222.0 / 4.6 75R870H (CGC→CAC) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base A (3/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

AGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACA  >  NC_002947/5336034‑5336322
                                                                                                                                                 |                                                                                                                                               
aGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCt                                                                                                                                              <  1:255755/149‑1 (MQ=255)
aGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCt                                                                                                                                              <  1:260836/149‑1 (MQ=255)
 gCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTg                                                                                                                                             <  5:199998/149‑1 (MQ=255)
 gCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTg                                                                                                                                             <  5:205109/149‑1 (MQ=255)
    accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTgg                                                                                                                                            <  5:29096/147‑1 (MQ=255)
    accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTgg                                                                                                                                            >  6:29096/1‑147 (MQ=255)
    accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGcc                                                                                                                                          <  3:12619/149‑1 (MQ=255)
       aGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGAAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAg                                                                                                                                       >  1:296968/1‑149 (MQ=255)
           cAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGc                                                                                                                                   >  6:95439/1‑149 (MQ=255)
           cAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGc                                                                                                                                   <  6:79328/149‑1 (MQ=255)
              cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCATCCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                <  2:251499/149‑1 (MQ=255)
               gaggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGAt                                                                                                                               >  6:134777/1‑149 (MQ=255)
                   agCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAc                                                                                                                           <  2:16531/149‑1 (MQ=255)
                    gCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAcc                                                                                                                          <  8:59504/149‑1 (MQ=255)
                             ccTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTcc                                                                                                                  >  4:188688/1‑148 (MQ=255)
                                 gCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGc                                                                                                             >  8:67884/1‑149 (MQ=255)
                                     aaGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCcg                                                                                                          <  2:51723/148‑1 (MQ=255)
                                            acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCt                                                                                                  >  1:343185/1‑149 (MQ=255)
                                            acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCt                                                                                                  >  7:253749/1‑149 (MQ=255)
                                             caAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCCTCGCCCGCGATGCCTc                                                                                                 >  4:300811/1‑149 (MQ=255)
                                                  tGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                >  6:23199/1‑112 (MQ=255)
                                                  tGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                <  5:23199/113‑2 (MQ=255)
                                                      aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCgg                                                                                        >  1:162537/1‑149 (MQ=255)
                                                         ggTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAg                                                                                     >  3:17906/1‑149 (MQ=255)
                                                             tgctgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGcc                                                                                 >  4:183882/1‑149 (MQ=255)
                                                             tgctgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGcc                                                                                 >  4:178054/1‑149 (MQ=255)
                                                               ctgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgt                                                                               <  3:188688/149‑1 (MQ=255)
                                                               ctgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgt                                                                               >  4:344052/1‑149 (MQ=255)
                                                                tgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtg                                                                              >  6:82211/1‑149 (MQ=255)
                                                                tgGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtg                                                                              >  5:13669/1‑149 (MQ=255)
                                                                  gCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgc                                                                             <  8:78537/148‑1 (MQ=255)
                                                                   cGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgcgc                                                                           <  5:134777/149‑1 (MQ=255)
                                                                          ccGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgc                                                                             <  1:19329/140‑1 (MQ=255)
                                                                          ccGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgc                                                                             >  2:19329/1‑140 (MQ=255)
                                                                           cGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCa                                                                                                                 <  6:307019/103‑1 (MQ=255)
                                                                           cGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCa                                                                                                                 >  5:307019/1‑103 (MQ=255)
                                                                             ccgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCCGACCATTCGGACGAAGGTGACCAGCAACGCCTTGCCCAGGTTCTCATCAACCTGCTATCCAACGACCGCGATGCCTCTCCGCCCCCCAGCCCCCTGCGCGTGCGCAGc                                                                  >  2:116559/2‑148 (MQ=255)
                                                                              cgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGa                                                                <  7:67884/149‑1 (MQ=255)
                                                                               gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                               >  8:64427/1‑149 (MQ=255)
                                                                               gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                               >  7:285685/1‑149 (MQ=255)
                                                                                    ttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTaa                                                          >  7:63941/1‑149 (MQ=255)
                                                                                    ttCGAAGTACAGTTCTTCAACCGTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTaa                                                          >  7:55313/1‑149 (MQ=255)
                                                                                             cAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCacac                                                 >  6:141603/1‑149 (MQ=255)
                                                                                              aGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                <  5:141603/149‑1 (MQ=255)
                                                                                               gTTCTTCAACCTTTGCGACCCGGAACATTGGGCCGAAGGTGACCCGCAGCTCCTGGCCCAGGTGCTGATCAACCAGCTCTCCAACGCCCGCGATGCCTGGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCg                                               <  1:116559/149‑1 (MQ=255)
                                                                                                ttcttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGt                                              >  6:137391/1‑149 (MQ=255)
                                                                                                    tcAACCTTTGCGACCCGGACCATTGGGACGAAGGTGACCCGCAGCGCCAGGCCGAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGGGCAGGGAGGTTTACGAGCACACCGTCGAc                                          >  6:148582/1‑149 (MQ=255)
                                                                                                     cAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGAcc                                         >  5:26529/1‑149 (MQ=255)
                                                                                                     cAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGAc                                          <  8:253749/148‑1 (MQ=255)
                                                                                                     cAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGAc                                          <  5:221973/148‑1 (MQ=255)
                                                                                                       aCCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCt                                        >  5:318883/1‑148 (MQ=255)
                                                                                                            tGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtc                                    <  5:95439/147‑1 (MQ=255)
                                                                                                                  ccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                            >  5:50356/1‑149 (MQ=255)
                                                                                                                  ccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                            <  3:344052/149‑1 (MQ=255)
                                                                                                                    ggACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAg                                                                   <  8:119982/108‑1 (MQ=255)
                                                                                                                    ggACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAg                                                                   >  7:119982/1‑108 (MQ=255)
                                                                                                                    ggACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACg                          <  8:55313/149‑1 (MQ=255)
                                                                                                                    ggACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACg                          <  4:296095/149‑1 (MQ=255)
                                                                                                                     gACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGa                         <  8:285685/149‑1 (MQ=255)
                                                                                                                      aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtc                                 >  4:298834/1‑140 (MQ=255)
                                                                                                                      aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtc                                 <  3:298834/140‑1 (MQ=255)
                                                                                                                      aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                        <  3:178054/149‑1 (MQ=255)
                                                                                                                      aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                        <  4:17906/149‑1 (MQ=255)
                                                                                                                      aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                        <  3:183882/149‑1 (MQ=255)
                                                                                                                        cATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAgg                      <  5:82211/149‑1 (MQ=255)
                                                                                                                          ttGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAgg                      >  3:80680/1‑147 (MQ=255)
                                                                                                                              gCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                        >  6:117768/1‑141 (MQ=255)
                                                                                                                              gCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                        <  5:117768/141‑1 (MQ=255)
                                                                                                                                cGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGa              >  2:28789/1‑149 (MQ=255)
                                                                                                                                     gcgACCCGCAGCGCCTGGCCCAGCTGCTGACCACCCTCCTCTCCAACCCCCGCGACGCCTCGCCGCCCGCCCGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGCTCGTCGAAGACGAAGGCACCGGGATTCCg         <  6:50356/147‑1 (MQ=255)
                                                                                                                                      tGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCga        <  1:28789/149‑1 (MQ=255)
                                                                                                                                       gACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCg         >  8:86230/1‑147 (MQ=255)
                                                                                                                                          gatcAGCACCTGGGCCAGGCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCgaagaa    >  8:297449/4‑149 (MQ=255)
                                                                                                                                            tcAGCACCTGGGCCAGGCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                            <  7:280965/122‑1 (MQ=255)
                                                                                                                                             cAGCACCTGGGCCAGGCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCCAGGTTAACGAGCCCCCCGTCGACCTGAACGACGAAGACGAAGGCACCGGGATGCCCAAAAAca  >  7:231459/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACA  >  NC_002947/5336034‑5336322

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