Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,328,758 C→A 8.3% A89D (GCC→GAC)  PP_2927 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,328,7580CA 8.3% 139.1 / 2.9 48A89D (GCC→GAC) PP_2927hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base A (2/2);  total (25/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 = 4.76e-01

GGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATCACCACCCGTGCATGTTCAGC  >  NC_002947/3328621‑3328903
                                                                                                                                         |                                                                                                                                                  
ggTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTgga                                                                                                                                        <  1:522103/150‑1 (MQ=255)
      cccTGAGCAACTGGCGGAAAAGACCGGCTGCCCACCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGACCTGGCGCTGGAGGACTc                                                                                                                                  >  1:1001387/1‑150 (MQ=255)
                 tGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACa                                                                                                                       >  2:21208/1‑150 (MQ=255)
                   gCGGAAAAGACCGGCTGCCACCCGCGGGTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATg                                                                                                                     >  1:585723/1‑150 (MQ=255)
                    cGGAAAAGACCGGCTGCCACCCGCGGCTGGGGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATgg                                                                                                                    <  1:292833/150‑1 (MQ=255)
                          aGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCcggcg                                                                                                              >  2:88543/1‑150 (MQ=255)
                           gACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCcggcgg                                                                                                             >  2:99219/1‑150 (MQ=255)
                                gCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCGTTTGTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGTCCGgcggcgcgg                                                                                                        >  1:79500/1‑150 (MQ=255)
                                 cTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGgcggcgcggc                                                                                                       >  1:656457/1‑150 (MQ=255)
                                     cACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAgtg                                                                                                   >  1:39642/1‑150 (MQ=255)
                                      aCCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAgtgg                                                                                                  >  2:800232/1‑150 (MQ=255)
                                          gcgGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTggc                                                                                              >  1:189169/1‑150 (MQ=255)
                                                   gCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGc                                                                                                                                                    <  2:961865/87‑1 (MQ=255)
                                                   gCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGc                                                                                                                                                    >  1:961865/1‑87 (MQ=255)
                                                   gCGAGAATGGCTGAATGCCGAGGCGGCATCGGGCTACATGGTCCATGACACGGGCCGTGTTCTGCGCGCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGtt                                                                                     <  1:85124/150‑1 (MQ=255)
                                                    cGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTc                                                                                    <  2:414053/150‑1 (MQ=255)
                                                      agaATGGTTGAATGCCCAGGCGGCATCGGGCAATATGGTTCATTACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGACCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCa                                                                                  <  2:1001387/150‑1 (MQ=255)
                                                       gaATGTTTGAATGCCCAGGCGGCATCGCGCTACATGGTCCATGACAGGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCACTGGCGCGGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCAt                                                                                 <  2:474838/150‑1 (MQ=255)
                                                       gaATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTTGTGGCGGCGCTGTTCCAt                                                                                 <  1:891877/150‑1 (MQ=255)
                                                       gaATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCAt                                                                                 <  2:1035536/150‑1 (MQ=255)
                                                       gaATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGGGGCGGCGCTGTTCCAt                                                                                 >  1:474838/1‑150 (MQ=255)
                                                            gTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaa                                                                            >  1:1021738/1‑150 (MQ=255)
                                                             ttGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTggcggc                                                                                           <  2:331826/134‑1 (MQ=255)
                                                             ttGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTggcggc                                                                                           >  1:331826/1‑134 (MQ=255)
                                                             ttGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaag                                                                           >  2:107438/1‑150 (MQ=255)
                                                             ttGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaag                                                                           >  1:491008/1‑150 (MQ=255)
                                                             ttGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaag                                                                           <  1:88543/150‑1 (MQ=255)
                                                               gAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaagga                                                                         >  1:137128/1‑150 (MQ=255)
                                                                                cGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCAt                                                                                 <  2:622545/125‑1 (MQ=255)
                                                                                cGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCAt                                                                                 >  1:622545/1‑125 (MQ=255)
                                                                                         tGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTg                                                                  <  1:301939/131‑1 (MQ=255)
                                                                                         tGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTg                                                                  >  2:301939/1‑131 (MQ=255)
                                                                                           gTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑At                                            >  1:1086487/1‑150 (MQ=255)
                                                                                            tCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaaggacaa                                                                      <  1:600898/124‑1 (MQ=255)
                                                                                            tCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaaggacaa                                                                      >  2:600898/1‑124 (MQ=255)
                                                                                             ccATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATgg                                          <  1:320673/150‑1 (MQ=255)
                                                                                               aTGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaaggac                                                                        <  1:794203/119‑1 (MQ=255)
                                                                                               aTGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATgacaaggac                                                                        >  2:794203/1‑119 (MQ=255)
                                                                                               aTGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATTCGCGGCG‑ATGGCg                                        <  2:1079818/150‑1 (MQ=255)
                                                                                                                  gctgcCGGAAGAACAGGCGATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGGTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATca                     <  1:590573/150‑1 (MQ=255)
                                                                                                                    tgcCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTggg                            >  2:20402/1‑141 (MQ=255)
                                                                                                                    tgcCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTggg                            <  1:20402/141‑1 (MQ=255)
                                                                                                                        tGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGCGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATg                                           >  2:760010/1‑122 (MQ=255)
                                                                                                                        tGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGCGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATg                                           <  1:760010/122‑1 (MQ=255)
                                                                                                                              aCAGGCCATGGACCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTg                                                                  >  1:393875/1‑94 (MQ=255)
                                                                                                                              aCAGGCCATGGACCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTg                                                                  <  2:393875/94‑1 (MQ=255)
                                                                                                                               cAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGCCGCCGCACCATTTTTGGCGGTGCTGTTCTATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATCACCACCCGTGCATg        <  2:585723/150‑1 (MQ=255)
                                                                                                                                    cATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCAT‑CGCCGCGCATGGCgg                                       >  2:1100860/1‑114 (MQ=255)
                                                                                                                                    cATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCAT‑CGCCGCGCATGGCgg                                       <  1:1100860/114‑1 (MQ=255)
                                                                                                                                     aTGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATCACCACCCGTGCATGTTCAGc  <  1:800232/150‑1 (MQ=255)
                                                                                                                                         cccTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATcacca                  <  2:1013458/130‑1 (MQ=255)
                                                                                                                                         cccTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATcacca                  >  1:1013458/1‑130 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCGGGCTACATGGTCCATGACAAGGGCCGTTTTGTGCTGCCTGAAGAACAGGCCATGGCCCTGGCGCTGGAGGACTCCCCGGCCTACATGGCCGGCGGCGCGGCAGTGGTGGCGGCGCTGTTCCATGACAAGGACAAGCTGGTAGCCGCCATGCGCGGCG‑ATGGCGGGTTGGCCTGGGGCGATCACCACCCGTGCATGTTCAGC  >  NC_002947/3328621‑3328903

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