Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,559 A→C 9.7% V174V (GTA→GTC)  PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,5590AC 9.7% 105.5 / 2.2 39V174V (GTA→GTC) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base C (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCG  >  NC_002947/4141414‑4141688
                                                                                                                                                 |                                                                                                                                 
ccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCGTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTCCTg                                                                                                                                >  8:238043/1‑149 (MQ=255)
ccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCACTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCACCGCCATTCGCCAGGTACTg                                                                                                                                >  1:105390/1‑149 (MQ=255)
 cGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTgg                                                                                                                               >  3:19970/1‑149 (MQ=255)
      aaaGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAg                                                                                                                           >  7:221056/1‑148 (MQ=255)
       aaGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGTACCTGGCGaatg                                                                                                                         <  5:53400/149‑4 (MQ=255)
        aGTGCCGGGGTGATCGCCCGCGCTGCCCAACCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCaa                                                                                                                        >  4:157059/1‑149 (MQ=255)
                ggTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGc                                                                                                                <  4:257089/149‑1 (MQ=255)
                  tGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAg                                                                                                              <  8:221056/149‑1 (MQ=255)
                     tCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGc                                                                                                           >  3:247959/1‑149 (MQ=255)
                           gcgcTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTacca                                                                                                      >  4:36420/1‑148 (MQ=255)
                             gcTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTaccacac                                                                                                   <  6:301078/149‑1 (MQ=255)
                                   cAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGTACCTGGCGaat                                                                                                                          >  2:305764/1‑118 (MQ=255)
                                   cAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGTACCTGGCGaat                                                                                                                          <  1:305764/120‑3 (MQ=255)
                                    aaGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGAcc                                                                                                                   >  4:273515/1‑126 (MQ=255)
                                    aaGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGAcc                                                                                                                   <  3:273515/126‑1 (MQ=255)
                                      gcggcgGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTc                                                                                          <  4:19970/149‑1 (MQ=255)
                                      gcggcgGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTc                                                                                          <  1:111353/149‑1 (MQ=255)
                                         gcggAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCt                                                                                       <  1:307751/149‑1 (MQ=255)
                                                 gggCCCAGACGCCAACGTCGATGCTGCCCTGGACCACTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCACCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCa                                                                               <  2:105390/149‑1 (MQ=255)
                                                    cccAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTc                                                                            <  3:157059/149‑1 (MQ=255)
                                                                    gATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGc                                                                                 <  8:187200/128‑1 (MQ=255)
                                                                    gATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCgcgc                                                            >  7:165019/1‑149 (MQ=255)
                                                                    gATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCgcgc                                                            >  4:182057/1‑149 (MQ=255)
                                                                         tgcCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGcc                                                                                <  8:230893/124‑1 (MQ=255)
                                                                         tgcCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGcc                                                                                >  7:230893/1‑124 (MQ=255)
                                                                         tgcCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGt                                                                   <  6:308547/137‑1 (MQ=255)
                                                                         tgcCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGt                                                                   >  5:308547/1‑137 (MQ=255)
                                                                           cccTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCgacgac                                                     <  8:165019/149‑1 (MQ=255)
                                                                             cTGGACCGCTGGTTCAGCCGTGAATACAAGGCCGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACAACTAATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCt                                                   >  5:257550/1‑149 (MQ=255)
                                                                                   cGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAg                                             >  1:76285/1‑149 (MQ=255)
                                                                                   cGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAg                                             >  2:69047/1‑149 (MQ=255)
                                                                                             gCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTg                                   <  3:182057/149‑1 (MQ=255)
                                                                                               cGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGcc                                 <  1:69047/149‑1 (MQ=255)
                                                                                                gTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACAACACCACTTATTCTCTGTTCGCCACTCAGGCCATGTACCGAGCCGACGTCCTGGGCAGTATCCAGGTTCCg                                >  3:106555/1‑149 (MQ=255)
                                                                                                                 aaCCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCgg               >  2:91616/1‑149 (MQ=255)
                                                                                                                  aCCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGc              <  3:36420/149‑1 (MQ=255)
                                                                                                                      cgcgCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCATCGACCCGCAGGGGTCCCCCCGCACTTATTCGGTTTTCGCCACTCCGGTCATGTCCGGCGCCGACGACCTGGGCAGTATCCTGGTGCCGACGCTGATCGCCACCGGCGAac          <  2:238316/149‑1 (MQ=255)
                                                                                                                             gTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAactcgactc   <  7:238043/149‑1 (MQ=255)
                                                                                                                              tCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAactcgactcg  >  4:60636/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                 
CCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAGGGCTACCACACCACTTATTCGCTGTTCGCCACTCAGGACATGTACCGCGCCGACGACCTGGGCAGTATCCAGGTGCCGACGCTGATCGCCACCGGCGAACTCGACTCG  >  NC_002947/4141414‑4141688

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: ‑