Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,331,722 G→A 9.6% L126L (CTG→CTA)  PP_3803 → cation ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,331,7220GA 9.6% 131.5 / 3.7 49L126L (CTG→CTA) PP_3803cation ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (3/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCTGTACAAC  >  NC_002947/4331576‑4331866
                                                                                                                                                  |                                                                                                                                                
cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                                                                                                               >  1:414108/1‑150 (MQ=255)
cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                                                                                                               >  2:1766181/1‑150 (MQ=255)
  cTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGCCGCTTGGCGCCAGCctgctgct                                                                                                                                             >  1:1001336/1‑150 (MQ=255)
     taCCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGAt                                                                                                                                          <  1:1426120/150‑1 (MQ=255)
                  ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAg                                                                                                                             >  1:1029451/1‑150 (MQ=255)
                   gCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                                                                                                                            <  1:1590777/150‑1 (MQ=255)
                           ttCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                                                                    >  1:1632788/1‑150 (MQ=255)
                                      cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCga                                                                                                         >  1:432273/1‑150 (MQ=255)
                                      cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAAACGTACCGGGCTGTAGCCGGATACACTGATTGGCCTGTTCCTGTCGGGGTGGCCGGCGCCAGCCACCAGCCTGCTGATGCTGCTGAGATGCAAGGTAAACGGACACAtacaccg                                                                                                         >  2:620046/1‑144 (MQ=255)
                                             ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                                                                  <  2:190901/150‑1 (MQ=255)
                                             ttCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCCGTCGGTGTCGCTGGCGCTTGGCCCCTGCGTACTACGGAGACGGTCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                                                                  >  1:1400394/1‑150 (MQ=255)
                                             ttCGGCATCCTGCTCAACTGCCTGCGCCACCGTCCCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg                                                                                                  <  2:167236/150‑1 (MQ=255)
                                                ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTg                                                                                               >  1:1131447/1‑150 (MQ=255)
                                                  cATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgtt                                                                                             <  1:1766181/150‑1 (MQ=255)
                                                      cTGCTCAACTTCCTGCGCAACCGTACCCGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcggt                                                                                         >  2:710172/1‑150 (MQ=255)
                                                        gCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACAACCTCGAGAACGTGCTgttcggttc                                                                                       >  2:64396/1‑150 (MQ=255)
                                                            aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGt                                                                                                                     >  1:1583223/1‑116 (MQ=255)
                                                            aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGt                                                                                                                     <  2:1583223/116‑1 (MQ=255)
                                                            aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTg                                                                                   >  2:863294/1‑150 (MQ=255)
                                                             aCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGc                                                                                  <  1:710172/150‑1 (MQ=255)
                                                              cTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCt                                                                                 <  2:1632788/150‑1 (MQ=255)
                                                                tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGa                                                                               >  1:1053558/1‑150 (MQ=255)
                                                                tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCATGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCt                                                                                                            <  2:1655998/121‑1 (MQ=255)
                                                                tCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCATGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCt                                                                                                            >  1:1655998/1‑121 (MQ=255)
                                                                      gcAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTa                                                                         >  1:1448977/1‑150 (MQ=255)
                                                                          ccGTACCGGGCTGTCGCCGGATACACTGATTAGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                                                                    <  1:1257438/103‑1 (MQ=255)
                                                                          ccGTACCGGGCTGTCGCCGGATACACTGATTAGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTg                                                                                                                    >  2:1257438/1‑103 (MQ=255)
                                                                                cGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGa                                                               >  2:333678/1‑150 (MQ=255)
                                                                                 gggCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGAc                                                              >  2:25380/1‑150 (MQ=255)
                                                                                   gCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGTAcca                                                            >  2:77622/1‑148 (MQ=255)
                                                                                     tGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                                                                                                                >  2:1730629/1‑80 (MQ=255)
                                                                                     tGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                                                                                                                <  1:1730629/80‑1 (MQ=255)
                                                                                     gatctccGGATACGCTGATTGGCGTGTTCCTGTCGGTGCTGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGc                                                          <  1:18181/145‑1 (MQ=255)
                                                                                     gatctccGGATACGCTGATTGGCGTGTTCCTGTCGGTATGCTTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGc                                                          <  1:18312/145‑1 (MQ=255)
                                                                                          ccGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGAACTTGCTGGTa                                                     >  1:636714/1‑150 (MQ=255)
                                                                                           cGGATACGCTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTAGAGAACGTGATGTTTGGTTAGGTGGTGGTGGTTAGTGTCGAGTTCTTGCtagatc                                                    >  2:18181/1‑145 (MQ=255)
                                                                                           cGGATACGCTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTACTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTAGAGAACGTACTATTAGGTTCTGTTCATACCTTTAGGGCCCAGtgcttgctagatc                                                    >  2:18312/1‑145 (MQ=255)
                                                                                                    tGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTc                                                                                                           <  2:1799161/86‑1 (MQ=255)
                                                                                                    tGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTc                                                                                                           >  1:1799161/1‑86 (MQ=255)
                                                                                                       ttGGCGTGTTCCTGTCGGTGTCGCTGGAGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGCTAGTGCCCAGGACTTGCTGGTACTGGGCAtcgtcg                                        <  1:620046/150‑1 (MQ=255)
                                                                                                            gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGACCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtg                                   <  2:432273/150‑1 (MQ=255)
                                                                                                             tgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgc                                  >  2:1150423/1‑150 (MQ=255)
                                                                                                             tgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgc                                  <  2:1131447/150‑1 (MQ=255)
                                                                                                                tCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTtcggttcggt                                                                                    <  1:397968/97‑1 (MQ=255)
                                                                                                                tCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTtcggttcggt                                                                                    >  2:397968/1‑97 (MQ=255)
                                                                                                                  cTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtc                                            <  2:1692846/135‑1 (MQ=255)
                                                                                                                  cTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtc                                            >  1:1692846/1‑135 (MQ=255)
                                                                                                                          gtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCtt                     <  1:863294/150‑1 (MQ=255)
                                                                                                                           tCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCGGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTg                    <  1:9332/150‑1 (MQ=255)
                                                                                                                             gCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCGTTGgc                  >  1:414521/1‑150 (MQ=255)
                                                                                                                                    cctGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCg           <  1:64396/148‑1 (MQ=255)
                                                                                                                                             aGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCTGTACAAc  >  1:767650/1‑150 (MQ=255)
                                                                                                                                                 tgctgctgATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGc          <  1:1433377/138‑1 (MQ=255)
                                                                                                                                                 tgctgctgATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGc          >  2:1433377/1‑138 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCTGTACAAC  >  NC_002947/4331576‑4331866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑