Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,634 C→T 49.3% T380T (ACC→ACT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,6340CT49.3% ‑6.4 / 47.4 37T380T (ACC→ACTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (9/9);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

AGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGG  >  NC_002947/195495‑195780
                                                                                                                                           |                                                                                                                                                  
aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCC                                                                                                                                          >  2:361097‑M1/1‑150 (MQ=255)
         ccAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGANCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGa                                                                                                                                 <  1:53990/150‑1 (MQ=40)
         ccAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATAGCTACACAGTCAGGTGCTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGa                                                                                                                                 <  2:324618/150‑1 (MQ=11)
                   cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCtgacgggctcgccgtcagg                                                                                                                       <  2:360543‑M1/150‑20 (MQ=255)
                         caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGa                                                                                                                 >  1:329516/1‑150 (MQ=39)
                                   aCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg                                                                                                       >  2:552215/1‑150 (MQ=38)
                                      gTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGATCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                    >  1:61115/1‑150 (MQ=255)
                                      gTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                    <  1:429190/150‑1 (MQ=17)
                                       tCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGCGCAACTGACCCCGAATCCGCCCCCGGCCCACACCACCACCAACGACTCGGTCGACACCACCACCCCGACCCTGACGGCGACCCCGTCCACCACCGAAGgcggcgccctcc                                                                                                   >  2:112937/1‑144 (MQ=11)
                                       tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCtgacgggctcgccgtcaggg                                                                                                                      <  2:611262‑M1/131‑21 (MQ=255)
                                       tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCtgacgggctcgccgtcaggg                                                                                                                      >  1:611262‑M1/1‑111 (MQ=255)
                                                 cATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         <  2:352625/150‑1 (MQ=255)
                                                  aTCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTg                                                                                        >  1:292987/1‑150 (MQ=40)
                                                  aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCTACTCATACCACGATGATCCAGTCGGTCCACACCACTATGGCGACCCTGACGGCGCGCCCGTCGGTCACCGAAGGCGGCGTGATTACTTACACCGTg                                                                                        >  1:141336/1‑150 (MQ=11)
                                                  aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCcgacgggctcgccgtcagggtcgcggtggt                                                                                                            >  1:92029‑M1/1‑100 (MQ=255)
                                                  aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCcgacgggctcgccgtcagggtcgcggtggt                                                                                                            <  2:92029‑M1/130‑31 (MQ=255)
                                                      agaACGCCACAGGCGGTAATTTCGAGCAACTGATCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAccc                                                                                    <  2:61115/150‑1 (MQ=255)
                                                      agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCAAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         >  2:586399/1‑145 (MQ=255)
                                                      agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCAAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt                                                                                         <  1:586399/145‑1 (MQ=255)
                                                         aCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTAGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa                                                                                 <  1:343287/150‑1 (MQ=9)
                                                                              aGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCgg                                                            <  1:361097/150‑1 (MQ=25)
                                                                               gCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGt                                                           >  1:245527/1‑150 (MQ=17)
                                                                                   cTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                       >  2:592016/1‑150 (MQ=255)
                                                                                   cTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAcc                                                       <  1:552215/150‑1 (MQ=255)
                                                                                       ccccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGa                                                   >  1:433438/1‑150 (MQ=255)
                                                                                         ccGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGAGACCCTGAGCAACCACGCCCGGACTCCGGTCACCGGGaca                                                 >  1:59629/1‑148 (MQ=11)
                                                                                                 gACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCaa                                         >  1:133480/1‑150 (MQ=255)
                                                                                                    cccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCcgacgggctcgccgtcagggtcgcggtagtggtgt                                                                                                       >  1:421734‑M1/1‑50 (MQ=255)
                                                                                                    cccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCcgacgggctcgccgtcagggtcgcggtagtggtgt                                                                                                       <  2:421734‑M1/85‑36 (MQ=255)
                                                                                                      cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCcgacgggctcgccgtcagggtcgcggtagtggtgt                                                                                                       <  1:535274‑M1/83‑36 (MQ=255)
                                                                                                      cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCcgacgggctcgccgtcagggtcgcggtagtggtgt                                                                                                       >  2:535274‑M1/1‑48 (MQ=255)
                                                                                                             gACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAt                             <  2:329516/150‑1 (MQ=255)
                                                                                                              aCCACGATCAACGACTCGGTCGACACCACAACAACGCCCTTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGAtt                            <  2:357228/150‑1 (MQ=14)
                                                                                                                cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTAc                          >  1:285937/1‑150 (MQ=255)
                                                                                                                cACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGCTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTAc                          >  1:408393/1‑150 (MQ=255)
                                                                                                                       aaCGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGaa                   <  1:399814/150‑1 (MQ=255)
                                                                                                                           aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg               <  1:112937/150‑1 (MQ=255)
                                                                                                                                        cactaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg               >  2:15018/1‑137 (MQ=255)
                                                                                                                                        cactaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg               <  1:15018/137‑1 (MQ=255)
                                                                                                                                        caccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAAACGTCAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAggg  >  1:346004/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
AGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGG  >  NC_002947/195495‑195780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: