Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,661,483 T→G 8.2% V52G (GTC→GGC)  PP_1456 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,661,4830TG 8.2% 133.3 / 2.7 46V52G (GTC→GGC) PP_1456membrane protein
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base G (2/2);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGC  >  NC_002947/1661337‑1661622
                                                                                                                                                  |                                                                                                                                           
tCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTc                                                                                                                                            <  7:137333/148‑1 (MQ=255)
   gATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGACCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGcc                                                                                                                                         <  3:72632/148‑1 (MQ=255)
   gATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGACCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCATGTACGcc                                                                                                                                         >  4:72632/1‑148 (MQ=255)
     ttGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCg                                                                                                                                           <  6:275055/144‑1 (MQ=255)
     ttGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTAGGTCATGGCGCCCTGGTCg                                                                                                                                           >  5:275055/1‑144 (MQ=255)
        ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATgac                                                                                                                                    <  5:173814/148‑3 (MQ=255)
            cATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggt                                                                                                                                 <  2:18803/147‑1 (MQ=255)
                agcagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtgg                                                                                                                               >  6:27804/1‑145 (MQ=255)
                agcagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtgg                                                                                                                               <  5:27804/145‑1 (MQ=255)
                       tAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTg                                                                                                                    >  8:109038/1‑149 (MQ=255)
                         gTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                  <  7:54422/149‑1 (MQ=255)
                               aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                            <  2:220798/149‑1 (MQ=255)
                               aTCAACTAAGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGCCCATTTTTGTTCATGTCGCCCTGCTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                            >  4:191980/1‑149 (MQ=255)
                                                cAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCa                                                                                           >  1:124659/1‑149 (MQ=255)
                                                cAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATga                                                                                                                                     >  2:193229/1‑106 (MQ=255)
                                                       ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                               >  1:157096/1‑138 (MQ=255)
                                                       ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                               <  2:157096/138‑1 (MQ=255)
                                                                   atcatcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCCCCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCAc                                                                        >  4:13741/1‑149 (MQ=255)
                                                                      atcTGGTCTGCCTGGGGAGGCCTGGTGATCATCCTGCGCGACATGGAGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGAGGCACTGGAGCTGTCCAATGTCACCTCGGAGCCGGCCATCGCCATGGCCAGCGACACTGc                                                                     >  4:230102/1‑149 (MQ=255)
                                                                         tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGtt                                                                  >  7:296973/1‑149 (MQ=255)
                                                                                   ggggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTcc                                                        >  8:19168/1‑149 (MQ=255)
                                                                                          ccTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                 >  4:96311/1‑149 (MQ=255)
                                                                                               tgtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGAc                                            <  3:13741/149‑1 (MQ=255)
                                                                                               tgtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGAc                                            <  8:36530/149‑1 (MQ=255)
                                                                                                     tCCTGCGCGACATGGTGGATACCCGTTTTGGTCATGGCGCCGTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTTTGTAACCTCGGAGCCGGCCATCGCCCTGGGCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                      <  3:191980/149‑1 (MQ=255)
                                                                                                     tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                      <  8:296973/149‑1 (MQ=255)
                                                                                                     tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                      >  1:217795/1‑149 (MQ=255)
                                                                                                       cTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                    >  5:309177/1‑149 (MQ=255)
                                                                                                       cTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                    >  4:93633/1‑149 (MQ=255)
                                                                                                        tGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATc                                                                                         <  6:162933/95‑1 (MQ=255)
                                                                                                        tGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATc                                                                                         >  5:162933/1‑95 (MQ=255)
                                                                                                          cgcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGtt                                 <  7:109038/149‑1 (MQ=255)
                                                                                                            cgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGtt                                 >  5:151162/1‑147 (MQ=255)
                                                                                                              aCATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGc                             <  3:93633/149‑1 (MQ=255)
                                                                                                                    tGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAcc                       >  8:297192/1‑149 (MQ=255)
                                                                                                                         aCCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCg                  >  6:91762/1‑149 (MQ=255)
                                                                                                                           ccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggct                                       >  6:250284/1‑126 (MQ=255)
                                                                                                                           ccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggct                                       <  5:250284/126‑1 (MQ=255)
                                                                                                                           ccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGc                >  2:23090/1‑149 (MQ=255)
                                                                                                                              ttttGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGACCGc             >  6:24771/1‑149 (MQ=255)
                                                                                                                                    tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTgcgc        >  3:99575/1‑148 (MQ=255)
                                                                                                                                      aTGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGa      <  5:91762/148‑1 (MQ=255)
                                                                                                                                      aTGGCGACCAGGGCGCCATGGAGGAGGCGCTGGTGCTGTCCTATCTCCGCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAc     >  4:318504/1‑149 (MQ=255)
                                                                                                                                       tGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAcc    >  7:124768/1‑149 (MQ=255)
                                                                                                                                       tGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAcc    >  6:325513/1‑149 (MQ=255)
                                                                                                                                         gcgcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGc  <  3:230102/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGC  >  NC_002947/1661337‑1661622

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