Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,177,331 C→T 9.8% V52M (GTG→ATG)  atpG ← ATP synthase subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,177,3310CT 9.8% 140.7 / 4.7 51V52M (GTG→ATG) atpGATP synthase subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/3);  total (24/27)
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.15e-01

CCAGGGCCTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGCGAATCTCTTTT  >  NC_002947/6177187‑6177473
                                                                                                                                                 |                                                                                                                                              
ccAGGGCCTTGAACAGGTTGGTGTTCAAGCCACCGCACGGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTgg                                                                                                                                            <  8:164561/149‑1 (MQ=255)
ccAGGGCCTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTgg                                                                                                                                            >  2:19066/1‑149 (MQ=255)
     ggcTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGGTGATCGGTCATCTGGccaac                                                                                                                                       <  2:111464/147‑5 (MQ=255)
      ccTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCt                                                                                                                                              >  8:155656/1‑141 (MQ=255)
      ccTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCt                                                                                                                                              <  7:155656/141‑1 (MQ=255)
            aCAGGTTGGAGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTcc                                                                                                                                >  3:33502/1‑149 (MQ=255)
               ggTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCcgcg                                                                                                                             <  6:106442/149‑1 (MQ=255)
                     tgtTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGa                                                                                                                       >  5:180329/1‑149 (MQ=255)
                         cAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACgg                                                                                                                    <  2:23016/148‑1 (MQ=255)
                          aaGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGc                                                                                                                   >  6:50900/1‑148 (MQ=255)
                          aaGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGc                                                                                                                   <  3:51089/148‑1 (MQ=255)
                            gCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTag                                                                                                                <  6:180329/149‑1 (MQ=255)
                                 cGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTagcagc                                                                                                            >  6:204010/1‑148 (MQ=255)
                                   cacaGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCAt                                                                                                         <  5:204010/149‑1 (MQ=255)
                                          cACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCAtt                                                                                                  >  6:187466/1‑149 (MQ=255)
                                             ggTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTtgt                                                                                               <  2:299404/149‑1 (MQ=255)
                                                cACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGcc                                                                                            <  4:33502/149‑1 (MQ=255)
                                                  cTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCtt                                                                                          <  3:317955/149‑1 (MQ=255)
                                                   tGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGGACGGGTGGCGGTATTCCGGGTTGGCGTTTGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTtgt                                                                                               <  5:238427/143‑1 (MQ=255)
                                                   tGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGGACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTtgt                                                                                               >  6:238427/1‑143 (MQ=255)
                                                   tGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTg                                                                                         >  5:87687/1‑149 (MQ=255)
                                                     cTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTgc                                                                                        >  7:41102/1‑148 (MQ=255)
                                                     cTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTgc                                                                                        <  8:41102/148‑1 (MQ=255)
                                                        accacGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCAt                                                                                    <  7:183483/149‑1 (MQ=255)
                                                          cacGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCt                                                                                  >  8:4019/1‑149 (MQ=255)
                                                             gATATAACCAGCGCGCTTTACAGGGCGGCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTg                                                                                >  1:4407/1‑149 (MQ=255)
                                                                      aGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCa                                                                      >  5:270802/1‑149 (MQ=255)
                                                                                aCAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCtt                                                                                 >  3:988/1‑128 (MQ=255)
                                                                                aCAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCtt                                                                                 <  4:988/128‑1 (MQ=255)
                                                                                     gc‑gcTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGgcgc                                                       <  4:169506/149‑1 (MQ=255)
                                                                                        gcTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTg                                                     >  2:164956/1‑149 (MQ=255)
                                                                                        gcTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTg                                                     <  7:4019/149‑1 (MQ=255)
                                                                                            gATCATGAACGGGTGGCGGTAATCCGGCATGCCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATATGTGCCTTGCGCATCTTGCTCACCGCCACTTTCTCCATGGCCCTGGTaa                                                 >  2:319372/1‑149 (MQ=255)
                                                                                              tCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAAtt                                               <  1:19066/149‑1 (MQ=255)
                                                                                                    aCGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGc                                         <  6:87687/149‑1 (MQ=255)
                                                                                                                 ttCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGc                            >  6:125056/1‑149 (MQ=255)
                                                                                                                    cGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTc                          <  2:4407/148‑1 (MQ=255)
                                                                                                                       gttggcgttggcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCaa                      >  7:120782/1‑149 (MQ=255)
                                                                                                                            cgttggcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTa                 <  2:172099/149‑1 (MQ=255)
                                                                                                                            cgttggcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTa                 <  5:50900/149‑1 (MQ=255)
                                                                                                                             gttggcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTAc                >  3:186069/1‑149 (MQ=255)
                                                                                                                              ttggcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACt               <  6:270802/149‑1 (MQ=255)
                                                                                                                              tccgcCAGGTGATCGGTCATCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACt               >  7:267121/4‑149 (MQ=255)
                                                                                                                                 gcCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGc             >  5:185309/1‑148 (MQ=255)
                                                                                                                                   cAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAAtttttt                                           >  3:263248/1‑116 (MQ=255)
                                                                                                                                   cAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAAtttttt                                           <  4:263248/116‑1 (MQ=255)
                                                                                                                                   cAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGCGaa          >  8:282486/1‑149 (MQ=255)
                                                                                                                                     ggtGATCGGTCATCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCttt                                  <  3:116150/121‑1 (MQ=255)
                                                                                                                                     ggtGATCGGTCATCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCttt                                  >  4:116150/3‑123 (MQ=255)
                                                                                                                                          ccGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGCGAATCTCttt   <  4:186069/149‑1 (MQ=255)
                                                                                                                                            ggtCATCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGCGAATCTCtttt  <  8:267121/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCAGGGCCTTGAACAGGTTGGTGTTCAAGCCACCGCACAGACCACGGTCACTGCTCACCACGATATAACCAGCGCGCTTTACAGGGC‑GCTCGATCATGAACGGGTGGCGGTATTCCGGGTTGGCGTTGGCCAGATGACCGATCACCTGGCGGATACGTTCCGCGTAAGGACGGCTAGCAGCCATGCGCATTTGTGCCTTGCGCATCTTGCTGACCGCCACTTTCTCCATGGCGCTGGTAATTTTTTGCGTGCTTTTGATGCTCGCAATCTTACTGCGAATCTCTTTT  >  NC_002947/6177187‑6177473

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