Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,682,324 A→C 6.6% R1523R (CGT→CGG PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,682,3240AC 6.6% 264.4 / 7.9 88R1523R (CGT→CGGPP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base C (4/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CGGTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGCCCCCGATAACCCGG  >  NC_002947/5682190‑5682472
                                                                                                                                      |                                                                                                                                                     
cGGTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCgg                                                                                                                                         >  4:122034/1‑149 (MQ=255)
cGGTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCgg                                                                                                                                         <  2:15364/149‑1 (MQ=255)
cGGTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCgg                                                                                                                                         <  6:195247/149‑1 (MQ=255)
 ggTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        <  3:272654/149‑1 (MQ=255)
      tctTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        <  4:184979/144‑1 (MQ=255)
      tctTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        >  3:184979/1‑144 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCCCGGcgccgcgccgtg                                                                                                                                        <  2:52248/142‑3 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCCCGGcgccgcgccgtg                                                                                                                                        >  1:52248/1‑140 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACTGCGCGGCGCCGGGCCAAGc                                                                                                                                  >  2:382864/1‑148 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAg                                                                                                                                   <  1:279525/147‑1 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGc                                                                                                                                  >  2:192505/1‑148 (MQ=255)
        ttCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGc                                                                                                                                  >  3:213835/1‑148 (MQ=255)
          ccAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAAc                                                                                                                               >  5:252136/1‑149 (MQ=255)
           cagcagGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        <  7:183017/139‑1 (MQ=255)
           cagcagGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        >  8:183017/1‑139 (MQ=255)
            agcagGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTg                                                                                                                             >  7:303015/1‑149 (MQ=255)
             gcagGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTgg                                                                                                                            >  8:310322/1‑149 (MQ=255)
             gcagGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTgg                                                                                                                            >  2:208230/1‑149 (MQ=255)
              cagGGCCATGGCATGTACAGCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTAACTTGCGTACGGTTACCGAATAGTCCACCACCATCACCAGCAAGGCAcgg                                                                                                                                                    <  8:329169/124‑1 (MQ=255)
              cagGGCCATGGCATCTACTCCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGc                                                                                                                           >  4:345959/1‑149 (MQ=255)
              cagGGCCATGGCATCTACACCAACCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCAcgg                                                                                                                                                    >  7:329169/1‑124 (MQ=255)
                 ggCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGtt                                                                                                                        <  7:183386/149‑1 (MQ=255)
                  gCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTg                                                                                                                       >  5:152557/1‑149 (MQ=255)
                  gCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGCCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATTGTCCACCACCATCACCAACAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTg                                                                                                                       <  1:315771/149‑1 (MQ=255)
                    cATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGc                                                                                                                     >  6:255812/1‑149 (MQ=255)
                     aTGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCACCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGc                                                                                                                     >  2:196706/1‑148 (MQ=255)
                       ggCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTAc                                                                                                                  >  3:351874/1‑149 (MQ=255)
                        gCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGTGTTGGCTAc                                                                                                                  <  1:208230/148‑1 (MQ=255)
                        gCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTAc                                                                                                                  <  2:130386/148‑1 (MQ=255)
                            cTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGcca                                                                                                            >  6:262627/1‑149 (MQ=255)
                              acacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccacc                                                                                                          >  3:337473/1‑149 (MQ=255)
                              acacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccacc                                                                                                          >  6:70550/1‑149 (MQ=255)
                               cacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccacca                                                                                                         >  8:111314/1‑149 (MQ=255)
                               cacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccacca                                                                                                         <  1:192505/149‑1 (MQ=255)
                               cacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccac                                                                                                           >  5:16487/1‑147 (MQ=255)
                               cacCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGccac                                                                                                           >  2:263819/1‑147 (MQ=255)
                                 ccATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAgg                                                                                                       <  4:351874/149‑1 (MQ=255)
                                  cATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAgg                                                                                                       >  8:254794/1‑148 (MQ=255)
                                  cATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGc                                                                                                      <  3:254686/149‑1 (MQ=255)
                                          gCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGc                                                                                              <  5:70550/149‑1 (MQ=255)
                                          gCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGc                                                                                              >  6:189208/1‑149 (MQ=255)
                                            cGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCa                                                                                            <  6:16487/149‑1 (MQ=255)
                                                  cACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTc                                                                                      >  8:366697/1‑149 (MQ=255)
                                                      cTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCgctg                                                                                  <  3:174245/149‑1 (MQ=255)
                                                         aTGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCgctgct                                                                                <  1:196706/148‑1 (MQ=255)
                                                          tGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTgg                                                                              >  4:91290/1‑149 (MQ=255)
                                                           gCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGc                                                                             >  8:222129/1‑149 (MQ=255)
                                                                 gcgcTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCcgcg                                                                         >  3:151728/1‑147 (MQ=255)
                                                                   gcTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGc                                                                     >  2:320553/1‑149 (MQ=255)
                                                                    cTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCt                                                                    <  7:310322/149‑1 (MQ=255)
                                                                       cagcagGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACTATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTg                                                                   >  5:130977/1‑147 (MQ=255)
                                                                       cagcagGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACTATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTg                                                                   <  6:130977/147‑1 (MQ=255)
                                                                              cggcTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCcagcag                                                          <  7:111314/149‑1 (MQ=255)
                                                                                    gTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCaa                                                     <  3:122034/148‑1 (MQ=255)
                                                                                     ttACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAg                                                    >  8:31187/1‑148 (MQ=255)
                                                                                     ttACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAg                                                    >  1:167929/1‑148 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAg                                                                                           <  5:398608/106‑1 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAg                                                                                           >  6:398608/1‑106 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAg                                                    <  1:55001/145‑1 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAg                                                    >  2:55001/1‑145 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATc                                                <  6:152557/149‑1 (MQ=255)
                                                                                        ccTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATc                                                <  5:189208/149‑1 (MQ=255)
                                                                                         cTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAAt                                                 >  2:255146/1‑147 (MQ=255)
                                                                                            gCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGc                                             <  5:36886/148‑1 (MQ=255)
                                                                                                  ggTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCAccc                                      >  5:418751/1‑149 (MQ=255)
                                                                                                    ttACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCACGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCCGCAGGTCAAGACTCAGCACCACCCgg                                    >  8:370327/1‑149 (MQ=255)
                                                                                                      aCCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGt                                                        <  7:130262/127‑1 (MQ=255)
                                                                                                      aCCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGt                                                        >  8:130262/1‑127 (MQ=255)
                                                                                                          aaTCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCg                              <  3:345959/149‑1 (MQ=255)
                                                                                                           aTCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGc                             <  8:303015/149‑1 (MQ=255)
                                                                                                            tCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        >  6:338115/1‑42 (MQ=255)
                                                                                                            tCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCggg                                                                                                                                        <  5:338115/42‑1 (MQ=255)
                                                                                                              gTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGcc                            >  7:238918/1‑147 (MQ=255)
                                                                                                              gTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCca                          <  4:213835/149‑1 (MQ=255)
                                                                                                                 caccaccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCcagca                       <  7:366697/149‑1 (MQ=255)
                                                                                                                   ccaccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAgt                     <  4:151728/149‑1 (MQ=255)
                                                                                                                      ccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTTGCCCAGCAGGTCAACAATCAGCGCCACCCGGCCAACGccc                           >  8:141565/1‑140 (MQ=255)
                                                                                                                      ccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGccc                           <  7:141565/140‑1 (MQ=255)
                                                                                                                      ccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTc                  <  7:31187/149‑1 (MQ=255)
                                                                                                                      ccaTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTc                  <  8:238918/149‑1 (MQ=255)
                                                                                                                       caTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGTCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTc                  >  1:110373/1‑149 (MQ=255)
                                                                                                                       caTCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCg                 <  5:255812/149‑1 (MQ=255)
                                                                                                                           aCCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGcccc             <  6:252136/149‑1 (MQ=255)
                                                                                                                            ccAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGccccc            >  6:274827/1‑149 (MQ=255)
                                                                                                                                cttggCCCGGCGCCGCGCCGTGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTc                  <  3:348322/136‑1 (MQ=255)
                                                                                                                                cttggCCCGGCGCCGCGCCGTGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTc                  >  4:348322/4‑139 (MQ=255)
                                                                                                                                 ttggCCCGGCGCCGCGCCGTGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGCCCCCGATaa       >  4:99875/3‑149 (MQ=255)
                                                                                                                                 ttggCCCGGCGCCGCGCCGTGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGCCCCCGATaa       >  1:3961/3‑149 (MQ=255)
                                                                                                                                      aCGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGCCCCCGATAACCCgg  >  2:46750/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                     
CGGTGCTCTTCCAGCAGGGCCATGGCATCTACACCATCCTTGGCCGTAAGCACGCTCATGCCGTGGCGCTCCAGCAGGCGGCTGGTTACCTTGCGTACGGTTACCGAATCGTCCACCACCATCACCAGCAAGGCACGGCGCGGCGCCGGGCCAAGCAACTGGCGTTGGCTACCGCCACCAGGCAAGCGGGCCAGCCGTCGCTGCTGGCCGCGCAGCTGGCCCAGCAGGTCAAGAATCAGCACCACCCGGCCATCGCCCAGCAGTGTCGCCCCCGATAACCCGG  >  NC_002947/5682190‑5682472

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: