Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 203,677 T→C 9.3% G3061G (GGT→GGC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947203,6770TC 9.3% 107.5 / 3.1 43G3061G (GGT→GGCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base C (2/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAAT  >  NC_002947/203537‑203825
                                                                                                                                            |                                                                                                                                                    
gagaACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGACTCTTTGGCGCGCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGacc                                                                                                                                             >  1:149290/1‑148 (MQ=14)
     ccTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTGGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCGGATcaacaa                                                                                                                                        >  1:183553/1‑150 (MQ=255)
      cTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAATTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCaacaac                                                                                                                                       >  2:269285/1‑150 (MQ=255)
                    ggCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                                                                                         <  2:407186/149‑1 (MQ=255)
                        gcgAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTc                                                                                                                     >  2:692356/1‑150 (MQ=40)
                          gAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCAc                                                                                                                   >  2:27225/1‑150 (MQ=40)
                              gTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACcaca                                                                                                                                   >  1:143584/1‑130 (MQ=255)
                              gTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACcaca                                                                                                                                   <  2:143584/130‑1 (MQ=255)
                                 gATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGc                                                                                                            >  2:63625/1‑150 (MQ=40)
                                               caccGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                                                              <  1:355002/150‑1 (MQ=40)
                                                    acgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtcatca                                                                                         <  1:27225/150‑1 (MQ=255)
                                                           cGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCC                                                                                  <  1:508198‑M1/150‑1 (MQ=255)
                                                            gTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCg                                                                                 >  2:160865/1‑150 (MQ=255)
                                                                 cAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACTCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCaa                                                                            >  2:226986/1‑150 (MQ=255)
                                                                       gACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGa                                                                                                                          <  1:751161/98‑1 (MQ=255)
                                                                       gACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGa                                                                                                                          >  2:751161/1‑98 (MQ=255)
                                                                         ccGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCgttggccggcacgg                                                                    <  1:706604‑M1/150‑15 (MQ=255)
                                                                           gCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCAAGTCCCCGAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCgttggccggcacggtg                                                                  <  1:692356‑M1/150‑17 (MQ=255)
                                                                            ccACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTgg                                                                 <  1:519694/150‑1 (MQ=37)
                                                                             gacCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCgttggccggcacggtgat                                                                <  1:695566‑M1/149‑19 (MQ=255)
                                                                               ctccGTCGGTCACCGAAGGCGGCGAGATCACTCACACCATCACCCTAACCAACAAAGATGGTCTGCCCATGAACAACCACAGGGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCa                                                              <  2:183553/148‑1 (MQ=11)
                                                                               ccccGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACTCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCa                                                              <  1:226986/150‑1 (MQ=17)
                                                                                  cGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTg                                                           <  1:474611/150‑1 (MQ=35)
                                                                                        tCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                     >  2:725699/1‑150 (MQ=39)
                                                                                         cACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGt                                                    <  1:63625/150‑1 (MQ=39)
                                                                                          aCCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTg                                                   >  2:850329/1‑150 (MQ=38)
                                                                                          aCCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTg                                                   >  1:150622/1‑150 (MQ=39)
                                                                                           ccGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGa                                                  <  2:572660/150‑1 (MQ=40)
                                                                                             gAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCgttggccggcacggtgatgacggtcttgcc                                                    <  2:182555‑M1/146‑31 (MQ=255)
                                                                                             gAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCgttggccggcacggtgatgacggtcttgcc                                                    >  1:182555‑M1/1‑116 (MQ=255)
                                                                                                        gaTCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACaa                                     <  2:150622/150‑1 (MQ=255)
                                                                                                                   ccATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCg                          <  1:160865/150‑1 (MQ=39)
                                                                                                                    cATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACGGTGCCgttggccggcacggtgatgacgg                                                           >  2:118977‑M1/1‑93 (MQ=255)
                                                                                                                    cATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGACGGCAAAACCGTCATCACGGTGCCgttggccggcacggtgatgacgg                                                           <  1:118977‑M1/116‑24 (MQ=255)
                                                                                                                    cATCACCCTGAACAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGGCGGCAAAACCGTCATAACCGTGACGACCAACGGCACCGTGGGCACTGCAACCGTGACTACCCCGGACAACGTCTACGTCgg                         >  2:145840/1‑150 (MQ=11)
                                                                                                                        aCCCTGCCCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTAcc                     <  2:647220/150‑1 (MQ=37)
                                                                                                                               ccAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACAGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCACGGACAACGTCTACGTCGGTACCAACGAcc              >  2:306681/1‑150 (MQ=255)
                                                                                                                                caacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACGGTGCCgttggccggcacggtga                                                                 >  2:287849‑M1/1‑81 (MQ=255)
                                                                                                                                caacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACGGTGCCgttggccggcacggtga                                                                 <  1:287849‑M1/98‑18 (MQ=255)
                                                                                                                                 aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCACCCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCt            <  1:410853/150‑1 (MQ=37)
                                                                                                                                 aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCattagccggcacggtga                                                                 <  2:748831‑M1/97‑18 (MQ=255)
                                                                                                                                 aacaaAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCattagccggcacggtga                                                                 >  1:748831‑M1/1‑80 (MQ=255)
                                                                                                                                           gTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAat  <  1:725699/150‑1 (MQ=34)
                                                                                                                                            |                                                                                                                                                    
GAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAAT  >  NC_002947/203537‑203825

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: