Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,231,527 C→T 8.7% V86V (GTG→GTA)  PP_3709 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,231,5270CT 8.7% 123.2 / 3.3 46V86V (GTG→GTA) PP_3709hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CCCTTTTCATGCCAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGAACGCACTACCAAGGTA  >  NC_002947/4231383‑4231676
                                                                                                                                                |                                                                                                                                                     
cccTTTTCATGCCAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCa                                                                                                                                                  >  1:1313087/1‑150 (MQ=255)
  cTTTTCATGCCAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGa                                                                                                                                                >  1:424767/1‑150 (MQ=255)
          gCCAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACct                                                                                                                                        <  1:1173424/150‑1 (MQ=255)
           ccAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGAc                                                                                                                                          <  1:392490/147‑1 (MQ=255)
           ccAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGAc                                                                                                                                          >  2:392490/1‑147 (MQ=255)
             aGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACcttct                                                                                                                                     >  1:1396722/1‑150 (MQ=255)
                gTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCttctttc                                                                                                                                  <  2:1396722/150‑1 (MQ=255)
                 tGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGTGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTTAGTACTGGCCACCCAGAATAGACCttctttct                                                                                                                                 <  1:683257/150‑1 (MQ=255)
                    ccaTCACACCACATAAGGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCAGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGa                                                                                                                              <  1:587468/148‑1 (MQ=255)
                             cacaTAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACga                                                                                                                     >  1:1179921/1‑150 (MQ=255)
                              acaTAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACgag                                                                                                                    >  1:1135670/1‑150 (MQ=255)
                                  aaTGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGAccc                                                                                                                <  2:1135670/150‑1 (MQ=255)
                                         ccTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCt                                                                                                         <  2:1179921/150‑1 (MQ=255)
                                          cTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTg                                                                                                        >  1:1289120/1‑150 (MQ=255)
                                           tAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGa                                                                                                       <  2:424767/150‑1 (MQ=255)
                                             aaTCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGaaa                                                                                                     <  2:1128131/150‑1 (MQ=255)
                                             aaTCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGaaa                                                                                                     <  2:1313087/150‑1 (MQ=255)
                                                  aCTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCTACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGAc                                                                                                >  2:817911/1‑150 (MQ=255)
                                                  aCTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCTACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGAc                                                                                                >  2:817741/1‑150 (MQ=255)
                                                                 tAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATc                                                                                         >  1:485485/1‑142 (MQ=255)
                                                                 tAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATc                                                                                         <  2:485485/142‑1 (MQ=255)
                                                                       aaTCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGa                                                                           >  2:1020739/1‑150 (MQ=255)
                                                                                    tCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTcccc                                                              >  2:489022/1‑150 (MQ=255)
                                                                                      tACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTc                                                                 >  1:721389/1‑145 (MQ=255)
                                                                                      tACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTc                                                                 <  2:721389/145‑1 (MQ=255)
                                                                                        cgtccgtcACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCtga                                                          <  1:1020739/150‑1 (MQ=255)
                                                                                         gtccgtcACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGATACCCGAAGCCTGAAAGAGACAGCCATCCTCaa                                                                                    <  2:873491/123‑1 (MQ=255)
                                                                                         gtccgtcACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCaa                                                                                    >  1:873491/1‑123 (MQ=255)
                                                                                         gtccgtcACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCtgat                                                         >  2:1344328/1‑150 (MQ=255)
                                                                                            cgtcACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGAc                                                      >  1:841100/1‑150 (MQ=255)
                                                                                                  cccGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCTACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCaaaaa                                                <  1:817911/150‑1 (MQ=255)
                                                                                                  cccGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCTACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCaaaaa                                                <  1:817741/150‑1 (MQ=255)
                                                                                                            cGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAgc                                      <  2:1289120/150‑1 (MQ=255)
                                                                                                             gAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAATGGTCTTCCCCCTGATGACCAAAAATGATCCCAgcc                                     >  2:421043/1‑150 (MQ=255)
                                                                                                             gAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAgcc                                     >  1:333406/1‑150 (MQ=255)
                                                                                                                    tAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATc                              >  2:789807/1‑150 (MQ=255)
                                                                                                                        gTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAgtcgtc                                                                              >  1:785241/1‑98 (MQ=255)
                                                                                                                        gTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAgtcgtc                                                                              <  2:785241/98‑1 (MQ=255)
                                                                                                                           gCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCAtt                       <  1:582806/150‑1 (MQ=255)
                                                                                                                             ggCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCa                     >  1:201133/1‑150 (MQ=255)
                                                                                                                               cAGTCCTCAGTACTGGCCACCCAGAATAGACCATCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCCTCATACACTCGTCCCAAAGTGCCACCTTCCGTAAGACCAAAAACGCTCCCAACCGCCGATCGCCCATTCATg                   >  2:311870/1‑150 (MQ=255)
                                                                                                                                 gTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGaa                 >  1:461033/1‑150 (MQ=255)
                                                                                                                                    cTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAg                                                                                   <  1:775413/81‑1 (MQ=255)
                                                                                                                                    cTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAg                                                                                   >  2:775413/1‑81 (MQ=255)
                                                                                                                                       aGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGAACGCACt           >  1:98406/1‑150 (MQ=255)
                                                                                                                                            tGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGAACGCACTACCaa      <  1:489022/150‑1 (MQ=255)
                                                                                                                                               ccACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCtg                                                           <  2:648992/94‑1 (MQ=255)
                                                                                                                                               ccACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCtg                                                           >  1:648992/1‑94 (MQ=255)
                                                                                                                                                cACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTc                      <  1:1320161/130‑1 (MQ=255)
                                                                                                                                                cACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTc                      >  2:1320161/1‑130 (MQ=255)
                                                                                                                                                cACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGAACGCACTACCAAGGTa  >  1:962335/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCCTTTTCATGCCAGCGTGGCGATCACACCACATAATGTGACCTAAATCGACTCACCATAGCGGCTAAGAAAATCTCGAAGACTTCTACGTCCGTCACCCCGTTTAAACGAAGCGATAACGTAGCGGCAGTCCTCAGTACTGGCCACCCAGAATAGACCTTCTTTCTTGACGCTCACGAGACCCGAAGCCTGAAAGAGACAGCCATCCTCAAGTCGTCCGAAAGGGTCTTCCCCCTGATGACCAAAAATGATCCCAGCCGCCGATCGCCCATTCATGAACGCACTACCAAGGTA  >  NC_002947/4231383‑4231676

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