Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,311,048 G→A 14.3% V267I (GTA→ATA) ‡ PP_2033 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,311,0480GA14.3% 101.8 / 8.6 42V267I (GTA→ATA) ‡PP_2033hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base A (4/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

ACGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAG  >  NC_002947/2310905‑2311194
                                                                                                                                               |                                                                                                                                                   
aCGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTAGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGc                                                                                                                                                <  7:227585/149‑1 (MQ=255)
           cATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGa                                                                                                                                    <  2:200257/149‑1 (MQ=255)
             tAACTTCCCGCTGGTGCGCAGTGGCCACGCCATCGTCCTGCCGGGCGTGCACTCGAATATTCGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGACTGTAGGC‑ACTCGTGTCgcg                                                                                                                                   >  8:222915/1‑146 (MQ=255)
                     cGCTGGTGCGCAGTGGCCACGACAGCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCcgc                                                                                                                          <  7:240231/149‑1 (MQ=255)
                         ggTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTAc                                                                                                                      <  6:256155/149‑1 (MQ=255)
                            gcgcAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATg                                                                                                                   >  4:307689/1‑149 (MQ=255)
                            gcgcAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATg                                                                                                                   >  4:302810/1‑149 (MQ=255)
                               cAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCa                                                                                                                <  8:89392/149‑1 (MQ=255)
                                       aCGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGt                                                                                                                               >  7:104037/1‑126 (MQ=255)
                                       aCGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGt                                                                                                                               <  8:104037/126‑1 (MQ=255)
                                                tgctgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGt                                                                                                                               >  4:125201/1‑117 (MQ=255)
                                                tgctgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGt                                                                                                                               <  3:125201/117‑1 (MQ=255)
                                                  ctgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGc                                                                                             <  3:302810/149‑1 (MQ=255)
                                                  ctgcCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGc                                                                                             <  3:307689/149‑1 (MQ=255)
                                                     ccGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGAAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTc                                                                                          <  2:94755/149‑1 (MQ=255)
                                                          cGTGCACTAGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAg                                                                                                                                                  >  5:299535/1‑89 (MQ=255)
                                                               aCTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCgg                                                                                        >  3:239054/1‑141 (MQ=255)
                                                               aCTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCgg                                                                                        <  4:239054/141‑1 (MQ=255)
                                                                 tCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGa                                                                              >  6:76124/1‑149 (MQ=255)
                                                                                    aCCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCa                                                                <  7:19758/144‑1 (MQ=255)
                                                                                    aCCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCa                                                                >  8:19758/1‑144 (MQ=255)
                                                                                              ccTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGCAAGGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtgat                                                   >  1:250171/1‑148 (MQ=255)
                                                                                              ccTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtgatga                                                 >  7:192042/1‑149 (MQ=255)
                                                                                               cTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCcgcgc                                                        <  4:285844/141‑1 (MQ=255)
                                                                                               cTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCcgcgc                                                        >  3:285844/1‑141 (MQ=255)
                                                                                                     gAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTgg                                            >  4:62236/1‑147 (MQ=255)
                                                                                                            tgtTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGcc                                    >  1:255658/1‑148 (MQ=255)
                                                                                                                      cAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTg                         >  6:217436/1‑149 (MQ=255)
                                                                                                                        aGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTggg                       >  3:217208/1‑149 (MQ=255)
                                                                                                                          ccGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGggcg                     >  3:95050/1‑149 (MQ=255)
                                                                                                                                  ccAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGAc              >  3:229362/1‑148 (MQ=255)
                                                                                                                                    cacgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTg                         >  7:189282/4‑135 (MQ=255)
                                                                                                                                    cacgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTg                         <  8:189282/132‑1 (MQ=255)
                                                                                                                                    cacgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAg           >  5:132082/4‑149 (MQ=255)
                                                                                                                                    aGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAg           >  1:313973/1‑149 (MQ=255)
                                                                                                                                    aGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAg           <  4:217208/149‑1 (MQ=255)
                                                                                                                                     acgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtgat                                                   >  7:298805/3‑108 (MQ=255)
                                                                                                                                     acgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGtgat                                                   <  8:298805/106‑1 (MQ=255)
                                                                                                                                      tgcgtgcCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGgc         <  5:76124/149‑1 (MQ=255)
                                                                                                                                      cgcgtgcCTACAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCTCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGgc         >  5:76020/2‑149 (MQ=255)
                                                                                                                                        cgtgcCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGgcgc       <  2:255658/149‑1 (MQ=255)
                                                                                                                                          tgcCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCa      >  3:12337/1‑148 (MQ=255)
                                                                                                                                               gTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAg  >  3:90644/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
ACGAGCGGCGGCATAACTTCCCGCTGGTGCGCAGTGGCCACGACATCGTGCTGCCGGGCGTGCACTCGAATATTGGGGGTGGTTACCCGGAGACCCTTCTGGAGCAGGTGTTGCTGTGCAAGCCGCAGTCCCAGTGCGTGCCTGTAGGC‑ACGCGTGTCGAGCAGTGCCGCGTACATGCCAAGGTGAGTGCCTTGCTGGCTTCGGCGTTTGCCGATATGGGCGAACCCAGGCCGCGCGTGATGACCTGGGAGGAGCCGATTGCCGGTGGGCGGCGGGACGAGGCGCAAAAG  >  NC_002947/2310905‑2311194

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