Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 657,272 T→C 10.2% V234A (GTT→GCT)  PP_0564 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947657,2720TC10.2% 132.5 / 5.0 49V234A (GTT→GCT) PP_0564hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/29);  new base C (3/2);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GCATCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTGGGCGCTGT  >  NC_002947/657125‑657418
                                                                                                                                                   |                                                                                                                                                  
gCATCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTg                                                                                                                                                  <  2:255905/150‑1 (MQ=255)
    cAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGAGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGTTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCg                                                                                                                                              <  1:699015/150‑1 (MQ=255)
       gTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGt                                                                                                                                           <  1:970472/150‑1 (MQ=255)
                ccGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCaa                                                                                                                                  <  1:485408/150‑1 (MQ=255)
                         gACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCAAGCGTTGGCCGGGTCACGCTCAACACCACGca                                                                                                                         >  2:63405/1‑149 (MQ=255)
                             cGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAg                                                                                                                     <  2:54066/150‑1 (MQ=255)
                               tCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGaa                                                                                                                   <  2:903233/150‑2 (MQ=255)
                                              cAAGGGCGCGATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGc                                                                                                    <  1:356841/150‑1 (MQ=255)
                                               aaGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGcc                                                                                                   <  1:63405/150‑1 (MQ=255)
                                                aGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAATACCGTTTCGTTCGTTGCCg                                                                                                  <  2:1013291/150‑1 (MQ=255)
                                                  ggCCCTATCTTCTATACCCGCTCTGACAACTTGAACAACTGGCTGGGTGACGTCGCCACCGGTCTGAGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGc                                                                                                <  2:604348/150‑1 (MQ=255)
                                                           ttctATACCCGCGCCGACAAACTGGACAACTGGATGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCATGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCAGGGCAGGcgccg                                                                                       >  2:1035666/1‑150 (MQ=255)
                                                             ctATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCa                                                                                     >  1:1092557/1‑150 (MQ=255)
                                                                taCCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCAGCCGGTTTGCTCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGt                                                                                  <  2:634989/150‑1 (MQ=255)
                                                                 aCCCGCGCCGACAACCTGAACAAATGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTg                                                                                 >  1:452181/1‑150 (MQ=255)
                                                                   ccGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACGCTGGCCgaca                                                                                                                                           <  1:344391/90‑4 (MQ=255)
                                                                   ccGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACGCTGGCCgaca                                                                                                                                           >  2:344391/1‑87 (MQ=255)
                                                                     gcgcCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCg                                                                                                                                              <  2:971424/85‑1 (MQ=255)
                                                                     gcgcCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCg                                                                                                                                              >  1:971424/1‑85 (MQ=255)
                                                                      cgcCGACAACCTGAACAAATGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATga                                                                            <  2:452181/150‑1 (MQ=255)
                                                                       gcCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATgag                                                                           >  1:1016248/1‑150 (MQ=255)
                                                                             aaCCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtg                                                                     >  2:426986/1‑150 (MQ=255)
                                                                               ccTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCTCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGATCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtggt                                                                   <  1:1152266/150‑1 (MQ=255)
                                                                               ccTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtggt                                                                   <  2:1092557/150‑1 (MQ=255)
                                                                               ccTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtggt                                                                   >  2:297514/1‑150 (MQ=255)
                                                                                cTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCAACACGCCGCAGGTGGATGAGGAGCtggtg                                                                  <  1:32300/150‑1 (MQ=255)
                                                                                 tGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGAt                                                                              >  2:549032/1‑137 (MQ=255)
                                                                                 tGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGAt                                                                              <  1:549032/137‑1 (MQ=255)
                                                                                         tggctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACGCTGGCCgac                                                                                                                                            >  1:338748/1‑65 (MQ=255)
                                                                                         tggctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACGCTGGCCgac                                                                                                                                            <  2:338748/67‑3 (MQ=255)
                                                                                           gctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATg                                                       <  2:1193776/150‑1 (MQ=255)
                                                                                                gTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGc                                                  <  1:297514/150‑1 (MQ=255)
                                                                                                 tGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCa                                                 >  1:748663/1‑150 (MQ=255)
                                                                                                     gTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATc                                             <  2:748663/150‑1 (MQ=255)
                                                                                                             ccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATgaggag                                                                        >  2:1023194/1‑115 (MQ=255)
                                                                                                             ccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATgaggag                                                                        <  1:1023194/115‑1 (MQ=255)
                                                                                                                   gggCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATgg                                                      <  2:341306/127‑1 (MQ=255)
                                                                                                                   gggCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATgg                                                      >  1:341306/1‑127 (MQ=255)
                                                                                                                      cTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCgg                                                                                                                                             >  2:483422/1‑37 (MQ=255)
                                                                                                                      cTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCgg                                                                                                                                             <  1:483422/37‑1 (MQ=255)
                                                                                                                              cGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTg                    <  1:128839/150‑1 (MQ=255)
                                                                                                                                cAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATc                                             <  2:389859/123‑1 (MQ=255)
                                                                                                                                cAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATc                                             >  1:389859/1‑123 (MQ=255)
                                                                                                                                  gcgcCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCCGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGt                     <  1:833056/145‑1 (MQ=255)
                                                                                                                                  gcgcCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCCGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGt                     >  2:833056/1‑145 (MQ=255)
                                                                                                                                       acccGGCCAACGCTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCt           >  2:764518/4‑150 (MQ=255)
                                                                                                                                        gTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGATCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTg          <  1:942030/150‑1 (MQ=255)
                                                                                                                                          cGGCCAGCGTTGGCCGGGGCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGCGTCAACGGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTggg        <  1:74844/150‑1 (MQ=255)
                                                                                                                                                gCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGATTCCGTCTTTGCCGGCCAGGCGCCCCAGGTGGATGAAGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACCTGTTCTATGAAGCCCGTGGCCAAGCCTGGGCGCTGt  <  1:1035666/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCATCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTGGGCGCTGT  >  NC_002947/657125‑657418

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