Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,120,471 T→A 10.6% intergenic (‑213/‑89) pepA ← / → PP_0982 cytosol aminopeptidase/YjgP/YjgQ family permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,120,4710TA10.6% 129.6 / 6.2 47intergenic (‑213/‑89)pepA/PP_0982cytosol aminopeptidase/YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base A (2/3);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CGCGTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCC  >  NC_002947/1120324‑1120610
                                                                                                                                                   |                                                                                                                                           
cgcgTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTg                                                                                                                                            >  6:124630/1‑149 (MQ=255)
 gcgTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTgg                                                                                                                                           <  1:165387/149‑1 (MQ=255)
   gtcgtcGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCTCAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGc                                                                                                                                          >  8:226538/1‑148 (MQ=255)
           gcACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCa                                                                                                                                 <  3:62080/149‑1 (MQ=255)
            cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctccat                                                                                                                                <  2:93062/149‑6 (MQ=255)
            cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctccat                                                                                                                                <  8:203081/149‑6 (MQ=255)
            cACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCAAAGGCGAAcctccat                                                                                                                                <  6:7417/149‑6 (MQ=255)
             aCCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCAtt                                                                                                                               >  4:77068/1‑149 (MQ=255)
                        aGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATg                                                                                                                    <  5:50072/149‑1 (MQ=255)
                        aGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATg                                                                                                                    <  2:11129/149‑1 (MQ=255)
                             ccAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACt                                                                                                                >  7:48401/1‑148 (MQ=255)
                             ccAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACt                                                                                                                >  4:143393/1‑148 (MQ=255)
                              cAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTaa                                                                                                              <  8:1299/149‑1 (MQ=255)
                              cAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTaa                                                                                                              <  3:143393/149‑1 (MQ=255)
                               aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                             <  3:77068/149‑1 (MQ=255)
                               aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                             <  8:120924/149‑1 (MQ=255)
                               aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                             >  5:116729/1‑149 (MQ=255)
                               aGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAg                                                                                                             <  5:124630/149‑1 (MQ=255)
                                gTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGt                                                                                                            >  1:65956/1‑149 (MQ=255)
                                     aGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCAtt                                                                                                        <  5:121649/148‑1 (MQ=255)
                                      gCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGt                                                                                                      <  2:65956/149‑1 (MQ=255)
                                       ccTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                         >  3:269003/1‑129 (MQ=255)
                                       ccTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCtt                                                                                                                         <  4:269003/129‑1 (MQ=255)
                                       ccTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGtt                                                                                                     >  6:142753/1‑149 (MQ=255)
                                           tgGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTgg                                                                                                  <  3:288774/148‑1 (MQ=255)
                                                    tGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCa                                                                                        <  1:173423/149‑1 (MQ=255)
                                                      cTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCAt                                                                                       <  6:116729/148‑1 (MQ=255)
                                                         aCAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATgg                                                                                     <  5:142753/147‑1 (MQ=255)
                                                             ccccGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGcc                                                                                              <  4:128408/134‑1 (MQ=255)
                                                             ccccGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACGTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGcc                                                                                              >  3:128408/1‑134 (MQ=255)
                                                              cccGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCg                                                                              >  5:134257/1‑149 (MQ=255)
                                                                  aTGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCt                                                                          <  2:172887/149‑1 (MQ=255)
                                                                    gtgCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGa                                                                        >  4:248368/1‑149 (MQ=255)
                                                                      gcgcCTTCGTACGCGTCAAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTga                                                                                         <  4:82251/130‑2 (MQ=255)
                                                                      gcgaCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTAGCCTTTGGCGAACCGGCATTGGACTTGGATGTACTAAGTCATTGTTTGGCGCCATTGc                                                                                         >  3:82251/1‑130 (MQ=255)
                                                                                     gCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTg                                                       <  6:134257/149‑1 (MQ=255)
                                                                                        aaGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGtt                                                    <  3:248368/149‑1 (MQ=255)
                                                                                          gTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGa                                                              <  8:29402/137‑1 (MQ=255)
                                                                                          gTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGa                                                              >  7:29402/1‑137 (MQ=255)
                                                                                                         aaTCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTa                                    <  7:226538/148‑1 (MQ=255)
                                                                                                         aaTCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTAt                                   >  3:190927/1‑149 (MQ=255)
                                                                                                             acaGGATAATGCGCGCACTTTACATGGAGGTTCGCGTTTGTCGATCCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCAATTGAATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATGTTCTTTAGTTATCTGt                               >  7:207930/1‑149 (MQ=255)
                                                                                                                         gcgcACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAccc                                                                  >  6:249480/1‑102 (MQ=255)
                                                                                                                         gcgcACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAAccc                                                                  <  5:249480/102‑1 (MQ=255)
                                                                                                                                        cggTTCGCCAAAGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAgcg    >  5:85762/2‑149 (MQ=255)
                                                                                                                                          gTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGACAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGcc  >  7:111623/1‑149 (MQ=255)
                                                                                                                                          gTTCGCCAAAGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGcc  >  5:13639/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CGCGTCGTCGCGCACCAGCGGAAAAGCTGCCAGTTTGAGCCTGTGGCAAGCGTGCTGACAACCCCGATGTGCGCCTTCGTACGCGGCCAAGTGACAGGGGCGGTCAATCACAGGATAATGCGCGCACTTTACATGGAGGTTCGCCTTTGGCGAACCGGCATTGGTCTTGGATGTACTAAGTCATTGTTTGGCGCCATTGCATGGCAGTCCGCCCTGACAACCCAGGAGTGTCTGGTTTGATCGTCTTTCGTTATCTGTCCCGCGAGGTCCTGGTCACCTTGAGCGCC  >  NC_002947/1120324‑1120610

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: