Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,079,678 C→T 8.7% intergenic (‑370/+64) PP_2688 ← / ← PP_2689 hypothetical protein/endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,079,6780CT 8.7% 119.1 / 2.7 46intergenic (‑370/+64)PP_2688/PP_2689hypothetical protein/endoribonuclease
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base T (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGG  >  NC_002947/3079541‑3079827
                                                                                                                                         |                                                                                                                                                     
gggCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGTATCGAGTGCCGCccgcgcg                                                                                                                                           >  1:1163189/1‑150 (MQ=255)
             aaaaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGTATCGAGTGCCGCCCGCGCGGCGCTCGAtgcat                                                                                                                              <  2:1163189/150‑5 (MQ=255)
              aaaaaGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTgcgc                                                                                                                             <  1:1019253/150‑1 (MQ=255)
                     gcgcTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGAGCGAACGcc                                                                                                                      <  2:1245257/150‑2 (MQ=255)
                         ttttGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAat                                                                                                                  >  2:738112/1‑150 (MQ=255)
                         ttttGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAat                                                                                                                  >  1:1117956/1‑150 (MQ=255)
                              aCCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                             >  2:1248115/1‑150 (MQ=255)
                                   ggCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCAc                                                                                                                                       <  1:1095793/119‑1 (MQ=255)
                                   ggCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCAc                                                                                                                                       >  2:1095793/1‑119 (MQ=255)
                                      gcgGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGCAGACCCGTGCAGGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACggg                                                                                                     >  2:495627/1‑150 (MQ=255)
                                       cgGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGgc                                                                                                    <  1:475294/150‑1 (MQ=255)
                                          gTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGcc                                                                                                 >  1:643968/1‑150 (MQ=255)
                                           ttttGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCt                                                                                                >  1:1268152/1‑150 (MQ=255)
                                               gAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAg                                                                                            <  1:1248115/150‑1 (MQ=255)
                                               gAACCGAAATGGGGCGCAGGCCGGTCCCTGGCGGAAGAAAACGGGAGTATTCCGGGAGAATGAAGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCGGAATATCGTGAGACGGGCGCGTTCac                                                                                            >  2:513004/1‑149 (MQ=255)
                                                  ccGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACgg                                                                                                      >  2:1057750/1‑137 (MQ=255)
                                                  ccGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACgg                                                                                                      <  1:1057750/137‑1 (MQ=255)
                                                  ccGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGc                                                                                         >  2:1094388/1‑150 (MQ=255)
                                                   cGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCg                                                                                        >  2:541365/1‑150 (MQ=255)
                                                    gAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGt                                                                                       <  2:499384/150‑1 (MQ=255)
                                                        gggggCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCaga                                                                                   <  2:350022/149‑1 (MQ=255)
                                                                aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGAATCGAGTGCCGCCCGCGCGGCGCTCGAtgcat                                                                                                                              >  2:403027/1‑95 (MQ=37)
                                                                aGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGAATCGAGTGCCGCCCGCGCGGCGCTCGAtgcat                                                                                                                              <  1:403027/99‑5 (MQ=37)
                                                                 ggCCGGTGCCGGGCGGGAGAAAAAGGGAGGATTCCGGGTGCAGGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAAGGCGTCAGAGGCAGGCGt                                                                          <  1:107306/150‑1 (MQ=255)
                                                                  gCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTa                                                                         <  2:643968/150‑1 (MQ=255)
                                                                   ccGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAg                                                                        <  2:1117956/150‑1 (MQ=255)
                                                                   ccGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACCTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAg                                                                        <  2:927526/150‑1 (MQ=255)
                                                                    cGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGa                                                                       >  2:1158122/1‑150 (MQ=255)
                                                                         ccTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGaccacc                                                                  <  1:1094388/150‑1 (MQ=255)
                                                                           tGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                             <  1:543904/105‑1 (MQ=255)
                                                                           tGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                             >  1:637298/1‑105 (MQ=255)
                                                                           tGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                             <  2:637298/105‑1 (MQ=255)
                                                                           tGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGt                                                                                                             >  2:543904/1‑105 (MQ=255)
                                                                                 gagaAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGGGCGGGGGGCACTCGATTGCGCGAACGCTGAATAATGTGAGACGGGCGCCTTCCGGGCGTCAGAGGCAGGAGTAGACCACCAGTTccac                                                          >  1:442798/1‑150 (MQ=255)
                                                                                 gagaAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAgg                                                                                           <  2:438854/117‑1 (MQ=255)
                                                                                 gagaAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAgg                                                                                           >  1:438854/1‑117 (MQ=255)
                                                                                      aaaGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTga                                                                                                           >  2:329656/1‑96 (MQ=255)
                                                                                      aaaGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTga                                                                                                           <  1:329656/96‑1 (MQ=255)
                                                                                         gggAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCAttt                                                  >  2:47405/1‑150 (MQ=255)
                                                                                         gggAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCAttt                                                  <  2:1268152/150‑1 (MQ=255)
                                                                                                   cGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGc                                        <  2:726337/150‑1 (MQ=255)
                                                                                                           tGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCAtt                                                   >  2:538364/1‑131 (MQ=255)
                                                                                                           tGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCAtt                                                   <  1:538364/131‑1 (MQ=255)
                                                                                                                 ggggCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGccacca                          >  1:202885/1‑150 (MQ=255)
                                                                                                                                 gcatCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTAAGTGGCAGGCGTAGACCCTCACGACCGCTTCCATTTGTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGtt          <  1:348719/150‑1 (MQ=255)
                                                                                                                                    tCGAGCGCCGCGCGGGCGGCGCTCGATTGCGCGAACGCTGAATATCGTGAGACGGGTGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTggg       <  1:738112/150‑1 (MQ=255)
                                                                                                                                         cgccgcGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGGAgg  <  2:442798/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
GGGCCGGCGCATCAAAAAAGGGCGCTTTTGACCCGGGCGCGGGTTTTGAACCGAAATGGGGCGCAGGCCGGTGCCTGGCGGGAGAAAAAGGGAGTATTCCGGGTGCATGACGTGGGGCGGAGGGGGCATGCATCGAGCGCCGCGCGGGCGGCACTCGATTGCGCGAACGCTGAATATCGTGAGACGGGCGCCTTCAGGGCGTCAGAGGCAGGCGTAGACCACCAGTTCCACCAGCATTTCTTCCCGGGCGAGGCCAGCCACCACTACGGCCGCGCGGTTGGGGTAGG  >  NC_002947/3079541‑3079827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: