Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,731 C→A 8.7% M329I (ATG→ATT)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294785,7310CA 8.7% 136.9 / 3.7 46M329I (ATG→ATT) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

GAACAATGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGAA  >  NC_002947/85584‑85875
                                                                                                                                                    |                                                                                                                                                
gAACAATGCACGGATCCCTCGGTGGTCCACCACCGCTCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCAt                                                                                                                                                >  1:640953/1‑150 (MQ=255)
    aaTGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGt                                                                                                                                            >  2:136282/1‑150 (MQ=255)
      tGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGa                                                                                                                                          <  1:269415/150‑1 (MQ=255)
           ggATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACggg                                                                                                                                     >  1:360940/1‑150 (MQ=255)
                 cTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccg                                                                                                                               >  1:871523/1‑150 (MQ=255)
                            accaACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGG‑CACAGGCGAACCCCGCGTCCTTCATCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTc                                                                                                                    <  1:136282/150‑1 (MQ=255)
                               aaCCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAgg                                                                                                                 <  1:113959/150‑1 (MQ=255)
                               aaCCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAgg                                                                                                                 >  2:113959/1‑150 (MQ=255)
                                    ccAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGAGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGa                                                                                                            >  1:906683/1‑150 (MQ=255)
                                    ccAGCAGATCAGCAGAATGCTGCCGGCCATGGAATGTTTCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGa                                                                                                            <  2:640506/150‑1 (MQ=255)
                                       gcggATCAGCAGACTGCTTGCAGCCATGGCCGGTAGCGGCGGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCg                                                                                                         <  2:360940/147‑1 (MQ=255)
                                        cagATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCgg                                                                                                        <  2:871523/150‑1 (MQ=255)
                                          gATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGt                                                                                                      <  2:35886/150‑1 (MQ=255)
                                          gATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGt                                                                                                      >  2:683893/1‑150 (MQ=255)
                                           aTCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTg                                                                                                     <  2:951689/150‑1 (MQ=255)
                                                 agAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCAGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGaa                                                                                               >  2:17382/1‑150 (MQ=255)
                                                  gAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAAc                                                                                              <  2:766044/150‑1 (MQ=255)
                                                          ccGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc                                                                                      <  1:859487/150‑1 (MQ=255)
                                                               cATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTa                                                                                 >  2:968962/1‑150 (MQ=255)
                                                                 tGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGCTAGGTAAc                                                                               >  2:748505/1‑150 (MQ=255)
                                                                    ccGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCag                                                                            <  1:328125/150‑1 (MQ=255)
                                                                         aTCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGaaa                                                                       >  2:712280/1‑150 (MQ=255)
                                                                            ggCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAAt                                                                    <  2:241302/150‑1 (MQ=255)
                                                                             gCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGGGCGGGTGAAGGAACCGGCAGACAGGTAACCAGCGAAAAATg                                                                   >  1:1083895/1‑150 (MQ=255)
                                                                                 gTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAGCGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCAcc                                                               >  1:105139/1‑150 (MQ=255)
                                                                                   gcgcGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAgg                                                                                                 >  1:799/1‑114 (MQ=255)
                                                                                   gcgcGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAgg                                                                                                 <  2:799/114‑1 (MQ=255)
                                                                                          cAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCAGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATAtt                                                      <  1:17382/150‑1 (MQ=255)
                                                                                           aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCATATAAAAATGCACCGGCAATATTc                                                     <  2:210625/150‑1 (MQ=255)
                                                                                           aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTc                                                     >  2:116047/1‑150 (MQ=255)
                                                                                           aGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTc                                                     >  2:666816/1‑150 (MQ=255)
                                                                                                   gcCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccgcc                                                                                                                             >  1:657884/1‑70 (MQ=255)
                                                                                                   gcCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGcgcccgcc                                                                                                                             <  2:657884/70‑1 (MQ=255)
                                                                                                                gAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAACGCACCGGCAATATTGGACAGGCCCTGTGCGCGTaa                                 >  2:9724/1‑149 (MQ=255)
                                                                                                                      ccccaccctcAGCGCTGAAACCGCCTCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTg                           <  2:906683/148‑1 (MQ=255)
                                                                                                                           cccacAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCg                      <  2:713258/150‑1 (MQ=255)
                                                                                                                                     tGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGaa                                                                                                   <  1:63868/63‑1 (MQ=255)
                                                                                                                                     tGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGaa                                                                                                   >  2:63868/1‑63 (MQ=255)
                                                                                                                                      gAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgt           <  1:116047/150‑1 (MQ=255)
                                                                                                                                      gAAAACACCCCGGCAATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgt           <  1:666816/150‑1 (MQ=255)
                                                                                                                                        aaaCACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAgttg         >  2:965455/1‑150 (MQ=255)
                                                                                                                                             ccccGGCCATCGGTGAACGGGCGCCCGGTGCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTg    >  1:308950/1‑150 (MQ=255)
                                                                                                                                              cccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc                                                                                      <  1:185516/67‑1 (MQ=255)
                                                                                                                                              cccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGAc                                                                                      >  2:185516/1‑67 (MQ=255)
                                                                                                                                               ccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGaa  >  1:1087719/1‑150 (MQ=255)
                                                                                                                                               ccGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAAGGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGaa  <  2:646206/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAACAATGCACGGATCCCTCGGTGGTCCACCAACCCCCAGCAGATCAGCAGAATGCTGCCGGCCATGGCCGGTATCGGCAGGTGCGCGATCAGGCCGGCGCCGGTCAC‑GGCGAACAACGCCACCCACAGCGCTGAAAACACCCCGGCCATCGGTGAACGGGCGCCCGCCTCGTAGCTCAGGCCCGAGCGGGTGAAGGAACCGGCAGACAGGTAACCAGAGAAAAATGCACCGGCAATATTCGACAGGCCCTGTGCGCGTATTTCCTGGTCCGGGTCGATCAGTTGTTCTGAA  >  NC_002947/85584‑85875

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