Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,266,125 G→A 11.1% intergenic (+42/+17) PP_2863 → / ← PP_2864 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,266,1250GA11.1% 100.0 / 2.9 39intergenic (+42/+17)PP_2863/PP_2864membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base A (2/2);  total (18/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 = 2.48e-01

AAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCA  >  NC_002947/3265982‑3266264
                                                                                                                                               |                                                                                                                                             
aaGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCg                                                                                                                                          <  8:31528/149‑1 (MQ=38)
           ccAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGa                                                                                                                               >  6:171283/1‑149 (MQ=255)
            cAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGAt                                                                                                                              >  4:61379/1‑149 (MQ=255)
                  gggtgCTGGGGGATTGGGCTGCTGCACCAGGTGTTCTGCGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACAt                                                                                                                         <  8:121263/146‑1 (MQ=255)
                   tgtgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGa                                                                                                                       <  7:177656/149‑1 (MQ=255)
                   tgtgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGa                                                                                                                       <  4:251782/149‑1 (MQ=255)
                     tgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      <  5:171283/148‑1 (MQ=255)
                      gCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAcc                                                                                                                     >  8:283047/1‑148 (MQ=255)
                         gggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAcc                                                                                                                 >  7:299468/1‑149 (MQ=255)
                          ggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      >  8:106302/1‑143 (MQ=255)
                          ggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      <  7:106302/143‑1 (MQ=255)
                                    ctgctgcTCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCaa                                                                                                       <  2:143550/148‑1 (MQ=255)
                                          cACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCtt                                                                                                 <  7:283047/148‑1 (MQ=255)
                                           aCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTgg                                                                                               <  8:299468/149‑1 (MQ=255)
                                                gtgtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAAc                                                                                          >  3:22479/1‑149 (MQ=255)
                                                      tGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATc                                                                                                                             >  8:33364/1‑108 (MQ=255)
                                                        ggCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCAc                                                                                <  3:61379/149‑1 (MQ=255)
                                                               tGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAg                                                                          >  6:195911/1‑148 (MQ=255)
                                                                 gCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAg                                                                          <  5:195911/146‑1 (MQ=255)
                                                                 gCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCt                                                                        <  4:22479/148‑1 (MQ=255)
                                                                   tGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACa                                                                                                                                     <  1:283086/87‑1 (MQ=255)
                                                                   tGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACa                                                                                                                                     >  2:283086/1‑87 (MQ=255)
                                                                     tttCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGc                                                                   >  2:259509/1‑149 (MQ=255)
                                                                          tGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATc                                                                                                                             >  4:247680/1‑88 (MQ=255)
                                                                            aCGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTCCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTCGGCAACCG‑‑ACGCACTACCAGCTGGTTCGACCGCt                                                            <  8:141053/149‑1 (MQ=255)
                                                                                gCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCaa                                                        <  2:44993/149‑1 (MQ=255)
                                                                                      tCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAcccc                                                                                                               >  2:273345/1‑90 (MQ=255)
                                                                                             cTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                              >  2:262924/1‑100 (MQ=255)
                                                                                             cTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                              <  1:262924/100‑1 (MQ=255)
                                                                                                             ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           >  2:269677/1‑149 (MQ=255)
                                                                                                             ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           <  1:269677/149‑1 (MQ=255)
                                                                                                               gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                                    >  4:64802/4‑106 (MQ=21)
                                                                                                               gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                                    <  3:64802/103‑1 (MQ=21)
                                                                                                                tACCCGCCTCATCGGCGGTAAACCCGCGTAATGGCCGGTAAATGTCGTATACATGTACAAACCCACCTTCAACCGCTTTGGCAAACGTCACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           <  1:259509/148‑1 (MQ=255)
                                                                                                                   cGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAgg                     <  2:252781/149‑1 (MQ=255)
                                                                                                                     gccctttCGCGGCTTTACCCGCGAAGAGGCCGGTACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggc                   >  6:16345/6‑149 (MQ=38)
                                                                                                                                ggTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCt        <  6:112072/149‑1 (MQ=255)
                                                                                                                                     aGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAcc    >  2:163688/1‑148 (MQ=255)
                                                                                                                                      gCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCa  >  7:167605/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                             
AAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCA  >  NC_002947/3265982‑3266264

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: ‑