Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,866,154 T→A 13.3% M206K (ATG→AAG)  PP_1669 → NLP/P60 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,866,1540TA13.3% 120.7 / 10.3 45M206K (ATG→AAG) PP_1669NLP/P60 family protein
Reads supporting (aligned to +/- strand):  ref base T (25/14);  new base A (4/2);  total (29/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.82e-01

CGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGGCTGCAG  >  NC_002947/1866026‑1866294
                                                                                                                                |                                                                                                                                            
cGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAATACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGgcgc                                                                                                                                         <  1:1154972/134‑1 (MQ=255)
cGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGgcgc                                                                                                                                         >  2:1154972/1‑134 (MQ=255)
               tGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCaaaa                                                                                                          <  1:1353617/150‑1 (MQ=255)
                      cAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGata                                                                                                                <  1:771132/137‑1 (MQ=255)
                      cAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGata                                                                                                                >  2:771132/1‑137 (MQ=255)
                            aTCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGc                                                                                             >  1:168546/1‑150 (MQ=255)
                                aCTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCcgcg                                                                                         >  1:244746/1‑150 (MQ=255)
                                          ccgccgCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGc                                                                               >  1:861529/1‑150 (MQ=255)
                                             ccgcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                            >  2:767576/1‑150 (MQ=255)
                                             ccgcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                            <  2:168546/150‑1 (MQ=255)
                                             ccgcAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                            <  1:721488/150‑1 (MQ=255)
                                                     gtggtgTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCCGTGATGCTTTGaaaaaaa                                                                    >  2:817348/1‑150 (MQ=255)
                                                              gcATCGACAGCCTCGGCGAACACTGCTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGAAACTCATATCCAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTccacc                                                           >  2:487141/1‑146 (MQ=255)
                                                                    aCAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTgc                                                     <  1:1394649/150‑1 (MQ=255)
                                                                         cTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                                >  2:1265272/1‑150 (MQ=255)
                                                                            ggtggtCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTAc                                             >  2:428427/1‑150 (MQ=255)
                                                                             gTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACa                                            <  2:861529/150‑1 (MQ=255)
                                                                                     gaCTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGt                                    >  2:977000/2‑150 (MQ=255)
                                                                                      aCTGGAGCAAGACCTTCATTGAAGCCAAGGGCGCCCTGGCCATGGCGCCGACCACCATCGCGCGCAACTAATATCAAAATGATGTACCATGCCGCGGCGATGATGCTTTGAAAAAAAAGCTCAACGCCGTGCGCATTTACaaaagacgcc                                   >  2:394029/1‑144 (MQ=255)
                                                                                        tGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtcta                                 >  2:465801/1‑150 (MQ=255)
                                                                                            gCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGtctaatct                             >  1:895226/1‑150 (MQ=255)
                                                                                               aGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTaaa                          >  1:729455/1‑150 (MQ=255)
                                                                                               aGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTCATCTTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTaaa                          <  1:533081/150‑1 (MQ=255)
                                                                                                gACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCAACGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAAt                         >  1:4519/1‑150 (MQ=255)
                                                                                                   cTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                            >  1:1196486/1‑96 (MQ=255)
                                                                                                   cTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCttt                                                                            <  2:1196486/96‑1 (MQ=255)
                                                                                                    ttCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGc                                      <  2:950354/133‑1 (MQ=255)
                                                                                                    ttCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGc                                      >  1:950354/1‑133 (MQ=255)
                                                                                                    ttCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtctc                     >  1:110894/1‑150 (MQ=255)
                                                                                                    ttCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATtctc                     >  2:476976/1‑150 (MQ=255)
                                                                                                       aTTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                                >  2:1162061/1‑120 (MQ=255)
                                                                                                       aTTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCAtt                                                <  1:1162061/120‑1 (MQ=255)
                                                                                                        ttGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTTCGCATTTACAGAAAGCGTCTAATCTAAATTCGCAAct                 >  1:517990/1‑150 (MQ=255)
                                                                                                        ttGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAAct                 >  2:336352/1‑150 (MQ=255)
                                                                                                            aGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttt             <  1:841621/150‑1 (MQ=255)
                                                                                                              ccAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCttttt           >  2:80158/1‑150 (MQ=255)
                                                                                                                aaGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCg         >  2:1318579/1‑150 (MQ=255)
                                                                                                                 tcggcgcCATGGCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCg                                     >  1:1417063/4‑121 (MQ=255)
                                                                                                                 tcggcgcCATGGCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCgg        >  1:906440/4‑150 (MQ=255)
                                                                                                                 tcggcgcCATGGCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTGACATAAAGCg                                     <  2:1417063/118‑1 (MQ=255)
                                                                                                                  gcgcgcCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGGc       <  1:476976/150‑1 (MQ=255)
                                                                                                                      gcCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGGCTgca   <  1:350332/150‑1 (MQ=255)
                                                                                                                       ccatGGCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGGCTgcag  >  2:1203256/4‑150 (MQ=255)
                                                                                                                           ggCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGcc                                                                                            <  2:449706/56‑1 (MQ=255)
                                                                                                                           ggCCAAGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGcc                                                                                            >  1:449706/1‑56 (MQ=255)
                                                                                                                                |                                                                                                                                            
CGGCATCTACCTGGGTGACAACCAGTTCATCCACTCCAGCAGCCGCCGCAGTGGTGGTGTGCGCATCGACAGCCTCGGTGATCGCTACTGGAGCAAGACCTTCATTGAAGCCAAGCGCGCCTTGGCCATGGCGCCGACCAACATCGCGCGCAACTGATATCAAAATGATGTACCCTGCCGCGGCGGTGATGCTTTGAAAAAAAAGCTCACCGCCGTGCGCATTTACAGAAAGCGTCTAATCTAAATTCTCAACTCTTTTTCGGCTGCAG  >  NC_002947/1866026‑1866294

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