Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,842 T→A 29.5% intergenic (+79/+674) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8420TA29.5% 48.2 / 27.2 44intergenic (+79/+674)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (20/11);  new base A (5/8);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

CGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTG  >  NC_002947/4740702‑4740983
                                                                                                                                            |                                                                                                                                              
cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                       <  6:34476/148‑1 (MQ=38)
  ggcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                     <  1:58777/148‑1 (MQ=38)
  cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                     <  3:236583/149‑1 (MQ=255)
   gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                    <  7:176402/149‑1 (MQ=38)
     cAGCTGTACACCGGCGAGCAGAACCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTCGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGc                                                                                                                                  <  1:142056/149‑1 (MQ=37)
              acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAg                                                                                                                          >  3:292348/1‑148 (MQ=255)
                    gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                    <  5:274112/132‑1 (MQ=255)
                    gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                    >  6:274112/1‑132 (MQ=255)
                      gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TGTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGGCACccca                                                                                                                  >  8:146251/1‑148 (MQ=255)
                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                      >  6:290736/1‑125 (MQ=37)
                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                      <  5:290736/125‑1 (MQ=37)
                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGTAGAGGCCGGCACAGACACCACaaaa                                                                                                               >  2:243362/1‑149 (MQ=255)
                          aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                      <  6:147623/124‑1 (MQ=37)
                          aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                                                                      >  5:147623/1‑124 (MQ=37)
                                gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCCGACACCACAAAACCTCA‑Gtt                                                                                                      >  3:174876/1‑149 (MQ=255)
                                                     cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAg                                                                                 >  6:95965/1‑149 (MQ=255)
                                                       cAATTACGCTGCAGTAGCCGAACGCAAGAAGGCGGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTGCCGCCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGcc                                                                               <  7:146251/149‑1 (MQ=255)
                                                           tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCCCCACAAAACCGCA‑GGTGTTCTCCGCCCGCTCACGCAGcccgtg                                                                           >  2:120239/1‑146 (MQ=255)
                                                                tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGCCCGCGAAGAGGGCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGACCCCTCACGCAGCCCTTTCAggg                                                                      >  5:153905/1‑149 (MQ=255)
                                                                     tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                                                  >  6:3279/1‑148 (MQ=255)
                                                                        ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTGATCGGGGGCTTGGCCGGGGAGGGGGCGGGCCAGCCCCCCCAAAACCGGC‑GGTTTTGGGCGGCCGCGGACGCCGCCCGTTGCGGGTATTGAAc                                                              >  2:93509/1‑149 (MQ=255)
                                                                          gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                                              <  7:231083/147‑1 (MQ=255)
                                                                          gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAc                                                              >  8:231083/1‑147 (MQ=255)
                                                                           aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaa                                                                                                                >  4:273699/1‑97 (MQ=37)
                                                                           aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaa                                                                                                                <  3:273699/97‑1 (MQ=37)
                                                                              ccacgatggCTGCCCGCGTGCTG‑CTTTTTTTTTTTGGCCGGTGCGGGCCTCTGCGCGGCATTGCCCGCTATGTGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCa                                                        <  8:198430/142‑1 (MQ=255)
                                                                                    atgCTGCCCTCCGGCATCCTTTTTTTTTTGTCCTTG‑GCCGGCCGCATCGGGGCCTTGCCCGGGAATAGGCCGGGAAATTAAACAAAAAACCTCATCTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                                                   <  1:93509/147‑1 (MQ=255)
                                                                                           ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTATCCGCCCGCTCACGCAGCCCGTTCAGGGTATTGAACGGCGCATCCACCACGAAcg                                           >  2:7092/1‑148 (MQ=255)
                                                                                                     gccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca                                  >  6:343702/4‑149 (MQ=255)
                                                                                                                     ctgTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGc                   <  5:3279/146‑1 (MQ=255)
                                                                                                                          ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccac             >  3:261087/1‑149 (MQ=255)
                                                                                                                           cGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCGTCGCCGGCACGCTaccacc            >  8:107287/1‑149 (MQ=255)
                                                                                                                            ggCCTCTTGGCGGCGTTGCCCGCGATGTTGCCGGCAGTGACACCTCAAAACCTCA‑GTTTTTCTCCGCACGCTGACGCAGCCCTTTCAGGGTAATGAACGGCGCATCCACCACTAACTTGTTGGCGACCTTGGGCGGTACGCTaccagc            <  6:153905/148‑3 (MQ=255)
                                                                                                                             gCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg          >  7:339226/1‑149 (MQ=255)
                                                                                                                               ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                         >  1:188846/1‑132 (MQ=255)
                                                                                                                               ctctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCg                         <  2:188846/132‑1 (MQ=255)
                                                                                                                                tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                       >  2:199988/1‑117 (MQ=39)
                                                                                                                                tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                       <  1:199988/117‑1 (MQ=39)
                                                                                                                                  ttCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg     >  8:20477/1‑149 (MQ=255)
                                                                                                                                  ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCT‑CTCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg    >  1:28546/1‑149 (MQ=255)
                                                                                                                                   tCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg    >  7:3969/1‑149 (MQ=255)
                                                                                                                                    cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc      <  5:343702/146‑1 (MQ=255)
                                                                                                                                    cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt   <  4:292348/149‑1 (MQ=255)
                                                                                                                                     gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCGGGGTATTGAACGGCGCATCCACCACGCACTTGTTGGCCACCATCGCCGGCCCGCTACCACCTGGCTCGgtg  >  4:229571/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCA‑GTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTG  >  NC_002947/4740702‑4740983

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