Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,799,638 T→C 12.3% V65A (GTA→GCA)  PP_5664 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,799,6380TC12.3% 77.8 / 3.4 33V65A (GTA→GCA) PP_5664hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base C (2/2);  total (17/16)
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.25e-01

GATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAA  >  NC_002947/4799489‑4799777
                                                                                                                                                     |                                                                                                                                           
gATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGt                                                                                                                                             >  2:1282419/1‑150 (MQ=255)
               ggggcgtcctTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAAGACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCtcattcat                                                                                                                              <  1:81968/144‑1 (MQ=255)
                         tCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACt                                                                                                                    <  2:86970/150‑1 (MQ=255)
                             gCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTgcgc                                                                                                                >  2:1029001/1‑150 (MQ=255)
                                tAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCg                                                                                                             <  1:119390/150‑1 (MQ=255)
                                      gcTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGc                                                                                                       <  2:401717/150‑1 (MQ=255)
                                                  aaTTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATc                                                                                           <  1:422692/150‑1 (MQ=255)
                                                       ccGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAg                                                                                                                            >  2:686420/1‑112 (MQ=255)
                                                       ccGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAg                                                                                                                            <  1:686420/112‑1 (MQ=255)
                                                         gaagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAgg                                                                                                                        <  2:408102/114‑1 (MQ=255)
                                                         gaagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAgg                                                                                                                        >  1:408102/1‑114 (MQ=255)
                                                          aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTa                                                                                                                                            <  1:1239262/93‑1 (MQ=255)
                                                          aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTa                                                                                                                                            >  2:1239262/1‑93 (MQ=255)
                                                          aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATc                                                                                   >  1:681227/1‑150 (MQ=255)
                                                            gaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTAt                                                                                            >  1:370878/1‑139 (MQ=255)
                                                            gaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTAt                                                                                            <  2:370878/139‑1 (MQ=255)
                                                                      gAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGcc                                                                       >  2:424676/1‑150 (MQ=255)
                                                                              gACCGAGCACGGATGGTGCGGACATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGa                                                               >  2:227730/1‑150 (MQ=255)
                                                                              gACCGAGCACGGATGGTGCGGACATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGa                                                               >  2:228030/1‑150 (MQ=255)
                                                                               aCCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGAt                                                              >  2:175217/1‑150 (MQ=255)
                                                                               aCCGAGCACGGATGGTGCGGATATGGTGATGGAATATCCGGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGAt                                                              <  1:1231755/150‑1 (MQ=255)
                                                                                        ggATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTg                                                     >  2:459875/1‑150 (MQ=255)
                                                                                           tGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTAc                                                  <  1:175217/150‑1 (MQ=255)
                                                                                             gTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACgg                                                >  1:984792/1‑150 (MQ=255)
                                                                                                   atatGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGacaa                                          <  1:1282419/150‑1 (MQ=255)
                                                                                                                      tGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTgg                                                                  >  2:345780/1‑107 (MQ=255)
                                                                                                                      tGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTgg                                                                  <  1:345780/107‑1 (MQ=255)
                                                                                                                            gTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTa                 >  2:487690/1‑150 (MQ=255)
                                                                                                                             tCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAg                >  1:1041461/1‑150 (MQ=255)
                                                                                                                              cAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGc               >  2:992768/1‑150 (MQ=255)
                                                                                                                                     cAAGACTACGAGGGCGGAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACt        >  2:690564/1‑150 (MQ=255)
                                                                                                                                        gACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGCAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTca     <  1:228030/150‑1 (MQ=255)
                                                                                                                                        gACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTca     <  1:227730/150‑1 (MQ=255)
                                                                                                                                           tACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACaa  <  2:1041461/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAA  >  NC_002947/4799489‑4799777

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: