Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 198,777 T→A 11.4% L1428Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947198,7770TA11.4% 119.3 / 7.5 44L1428Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base A (2/3);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGA  >  NC_002947/198631‑198915
                                                                                                                                                  |                                                                                                                                          
cAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTgg                                                                                                                                          <  1:754630/149‑1 (MQ=255)
cAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGt                                                                                                                                         <  1:89140/150‑1 (MQ=255)
cAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATAGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTgg                                                                                                                                          >  2:754630/1‑149 (MQ=255)
 aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                        <  2:1246066/150‑1 (MQ=255)
   cACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGt                                                                                                                                      <  1:263252/150‑1 (MQ=255)
     ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGa                                                                                                                                    <  1:5867/150‑1 (MQ=255)
     ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGa                                                                                                                                    >  2:60029/1‑150 (MQ=255)
     ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTCGTGa                                                                                                                                    <  1:1189538/150‑1 (MQ=255)
     ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTCGTGa                                                                                                                                    <  2:440606/150‑1 (MQ=255)
           tCACCGACGAGCCAGGTACGCCAGGAAACGACGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                              <  1:398180/150‑1 (MQ=255)
               cgacgaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                          >  2:918167/1‑150 (MQ=255)
               cgacgaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                          >  1:1124738/1‑150 (MQ=255)
                         ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                >  1:488586/1‑150 (MQ=255)
                            aCGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTc                                                                                                             <  1:110263/150‑1 (MQ=255)
                              gCCAGGCAACGAAGGCACAAGCGGCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                           <  1:1029930/150‑1 (MQ=17)
                                cAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTAAGGTCACCa                                                                                                         >  2:1186790/1‑150 (MQ=255)
                                  ggCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATc                                                                                                       <  2:591533/150‑1 (MQ=255)
                                  ggCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATc                                                                                                       <  1:918167/150‑1 (MQ=255)
                                     aaCGAAGGCGATATCGTCAAGGTCACCATCACGGGCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                    <  1:447725/150‑1 (MQ=255)
                                     aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                    <  2:1124738/150‑1 (MQ=255)
                                     aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                    <  1:279889/150‑1 (MQ=255)
                                       cGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGc                                                                                                  <  1:60029/150‑1 (MQ=255)
                                               atatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGaaa                                                                                          <  1:240095/150‑1 (MQ=255)
                                                  tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                       <  2:893431/150‑1 (MQ=255)
                                                     tCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCCCGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                       <  2:1041530/147‑1 (MQ=255)
                                                     tCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTAACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                       >  1:1041530/1‑147 (MQ=255)
                                                                  cACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCa                                                                       <  2:1081373/150‑1 (MQ=255)
                                                                                 ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGa                                                        >  1:316081/1‑150 (MQ=255)
                                                                                    ggTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGa                                                     >  1:759637/1‑150 (MQ=255)
                                                                                            agaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATc                                             >  2:231532/1‑150 (MQ=255)
                                                                                             gaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCa                                            <  2:789909/150‑1 (MQ=255)
                                                                                                 gTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGAt                                        >  2:927098/1‑150 (MQ=255)
                                                                                                     aaCCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACTACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                    <  1:68111/150‑1 (MQ=255)
                                                                                                            cTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGAt                             <  1:231532/150‑1 (MQ=255)
                                                                                                                  tgtgCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGACGGCGTAACGAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCt                       >  1:477216/1‑150 (MQ=255)
                                                                                                                       acaTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCa                  >  1:828307/1‑150 (MQ=255)
                                                                                                                         atgaaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATc                >  1:1282491/4‑150 (MQ=255)
                                                                                                                         aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATc                >  1:653461/1‑150 (MQ=255)
                                                                                                                               acCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcg          >  1:736322/1‑150 (MQ=255)
                                                                                                                                      agggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                    >  2:262787/5‑117 (MQ=14)
                                                                                                                                      agggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                    <  1:262787/113‑1 (MQ=14)
                                                                                                                                       gggtCACGACCAGGTCGTAACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGa  >  2:1198698/4‑150 (MQ=255)
                                                                                                                                             cGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTg    >  1:929732/1‑142 (MQ=255)
                                                                                                                                             cGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTg    <  2:929732/142‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
CAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGA  >  NC_002947/198631‑198915

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