Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,985,967 A→G 10.6% I82V (ATC→GTC)  oplB → 5‑oxoprolinase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,985,9670AG10.6% 123.4 / 4.4 47I82V (ATC→GTC) oplB5‑oxoprolinase B
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (2/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGCCC  >  NC_002947/3985820‑3986115
                                                                                                                                                   |                                                                                                                                                    
tCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACAt                                                                                                                                                     >  4:273540/1‑149 (MQ=255)
  cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATc                                                                                                                                                    >  1:156001/1‑148 (MQ=255)
  cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATc                                                                                                                                                    <  2:156001/148‑1 (MQ=255)
       tgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                              <  4:248717/149‑1 (MQ=255)
            aaCCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGa                                                                                                                                         <  3:273540/149‑1 (MQ=255)
                   accACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATc                                                                                                                                  <  6:180726/148‑1 (MQ=255)
                    ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATc                                                                                                                                  >  1:189310/1‑148 (MQ=255)
                      aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                              >  1:294387/1‑134 (MQ=255)
                      aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                              <  2:294387/134‑1 (MQ=255)
                      aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAc                                                                                                                               >  7:246229/1‑149 (MQ=255)
                             ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                                                              >  5:55246/1‑127 (MQ=255)
                             ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                                                              <  6:55246/127‑1 (MQ=255)
                             ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGCACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                         <  1:50947/148‑1 (MQ=255)
                                cAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAc                                                                                                                      <  8:246229/148‑1 (MQ=255)
                                 aaCTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAccc                                                                                                                    >  4:112599/1‑149 (MQ=255)
                                     cGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGcc                                                                                                                >  7:255133/1‑149 (MQ=255)
                                         gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTag                                                                                                                               >  3:8323/1‑129 (MQ=255)
                                         gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAc                                                                                                                               <  4:8323/130‑1 (MQ=255)
                                            cGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAg                                                                                                          <  2:54490/148‑1 (MQ=255)
                                                 aCTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGc                                                                                                    <  2:189310/149‑1 (MQ=255)
                                                             tGATCCGGCAGGGCCAGGACATGCCCATCACCCTGGGCGCGCTTATCCCGGCCATTAAGTTCGTGGCCGAGTTATTCGGCGACGACATCGCCGAAGGCGACGTTATATACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAc                                                                                        <  6:270457/149‑1 (MQ=255)
                                                                        ggCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGCCCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTa                                                                             >  4:235103/1‑149 (MQ=255)
                                                                               aCATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCg                                                                      <  2:300395/149‑1 (MQ=255)
                                                                                      cATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTAtt                                                               <  3:112599/149‑1 (MQ=255)
                                                                                           aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                                                             <  3:237332/82‑1 (MQ=255)
                                                                                           aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                                                             >  4:237332/1‑82 (MQ=255)
                                                                                           aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAg                                                          >  5:214124/1‑149 (MQ=255)
                                                                                            ccTGGGCTCGCTGATCCCGGCCATGAAGTGGGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGTCGACGTGATCTACCACAACGACCCAGCCTACAAGGGCAGACATATTCTCGACTGATGCATGTACAAACCGGTGTATTACCAg                                                          >  1:26565/1‑148 (MQ=255)
                                                                                            ccTGGCCTGCCTGATCCCGGCCATGACCGGCGGGGCCGAGTTCTTCGGCCACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAgg                                                         <  8:103747/149‑1 (MQ=255)
                                                                                                ggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGAc                                                       >  1:44078/1‑147 (MQ=255)
                                                                                                   tCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtgt                                                                  >  4:202805/1‑133 (MQ=255)
                                                                                                   tCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtgt                                                                  <  3:202805/133‑1 (MQ=255)
                                                                                                     gCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTgg                                                <  8:255133/149‑1 (MQ=255)
                                                                                                                 cATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACg                                                                                                                                       >  5:331153/1‑50 (MQ=255)
                                                                                                                 cATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACg                                                                                                                                       <  6:331153/50‑1 (MQ=255)
                                                                                                                      aGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGgt                                               <  8:178004/133‑1 (MQ=255)
                                                                                                                      aGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGgt                                               >  7:178004/1‑133 (MQ=255)
                                                                                                                          cGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAg                            <  2:44078/148‑1 (MQ=255)
                                                                                                                                cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCAcc                     >  5:177283/1‑149 (MQ=255)
                                                                                                                                    ttcttcGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGAc                 >  7:44884/1‑149 (MQ=255)
                                                                                                                                       ttcGGCGACGACATCGCCGAAGGCGACGTGATCTACAACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTATGGACCGTTTGAAAGGCCCACCTGAcct               >  8:264183/1‑147 (MQ=255)
                                                                                                                                         cGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAg                                                    <  1:119165/109‑1 (MQ=255)
                                                                                                                                         cGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAg                                                    >  2:119165/1‑109 (MQ=255)
                                                                                                                                          ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                                 <  5:101546/111‑1 (MQ=255)
                                                                                                                                          ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                                 >  6:101546/1‑111 (MQ=255)
                                                                                                                                                gacATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATcggcgg     >  7:300704/1‑149 (MQ=255)
                                                                                                                                                  cATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGAcc                                     >  2:257349/1‑115 (MQ=255)
                                                                                                                                                   aTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGccc  >  8:38362/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGGCCC  >  NC_002947/3985820‑3986115

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