Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,110 A→G 15.7% intergenic (+26/‑472) PP_4308 → / → PP_4309 AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,900,1100AG15.7% 114.5 / 14.1 51intergenic (+26/‑472)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base G (4/4);  total (26/25)
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.91e-01

GCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCA  >  NC_002947/4899964‑4900251
                                                                                                                                                  |                                                                                                                                             
gCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTGGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                            <  1:679484/150‑1 (MQ=255)
            gATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGcac                                                                                                                                          <  1:175185/140‑3 (MQ=255)
            gATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGcac                                                                                                                                          >  2:175185/1‑138 (MQ=255)
                 agcagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGcc                                                                                                                           >  1:692825/1‑150 (MQ=255)
                     gcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTa                                                                                                                       >  1:16708/1‑150 (MQ=255)
                      cTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTAc                                                                                                                      <  2:185546/150‑1 (MQ=255)
                            acCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaaa                                                                                                                >  2:213398/1‑150 (MQ=255)
                              ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                              >  1:462617/1‑150 (MQ=255)
                              ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                              >  1:348582/1‑150 (MQ=255)
                               cTGCTCCACCTCATGAATTTGCTGGCCGAGACCAATACCTCGATTATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGc                                                                                                             <  1:354776/150‑1 (MQ=255)
                                           aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTAtttttt                                                                                                 <  2:462617/150‑1 (MQ=255)
                                              aaTGGGCTGGCCGAGACCAATACCGCGAGCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTCATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTTACCGCTAAGAGGCCACCACAGGAAAAGCTGAGTATTTTTTGTa                                                                                              <  2:1051341/150‑1 (MQ=255)
                                               aTGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                              >  1:413666/1‑133 (MQ=255)
                                               aTGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAg                                                                                                              <  2:413666/133‑1 (MQ=255)
                                                   gCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACaa                                                                                         <  2:16708/150‑1 (MQ=255)
                                                        ccGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACg                                                                                    >  2:1118293/1‑150 (MQ=255)
                                                         cGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGt                                                                                   >  2:716244/1‑150 (MQ=255)
                                                          gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAgcac                                                                                                                      >  2:627443/1‑111 (MQ=25)
                                                          gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTACTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGCAc                                                                                                                      <  1:627443/114‑1 (MQ=25)
                                                                aaTACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATa                                                                            >  1:738832/1‑150 (MQ=255)
                                                                   aCCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                         <  1:213398/150‑1 (MQ=255)
                                                                    ccGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTg                                                                        >  2:681288/1‑150 (MQ=255)
                                                                     cgcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCa                                                                                                                          >  1:1054336/1‑99 (MQ=37)
                                                                     cgcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCa                                                                                                                          <  2:1054336/99‑1 (MQ=37)
                                                                     cgcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaaaa                                                                                                               <  1:952091/110‑1 (MQ=255)
                                                                     cgcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaaaa                                                                                                               >  2:952091/1‑110 (MQ=255)
                                                                             tcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtg                                                               >  2:814069/1‑150 (MQ=255)
                                                                               aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtg                                                             >  1:310080/1‑150 (MQ=255)
                                                                                agaagaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                            >  1:640043/1‑150 (MQ=255)
                                                                                       cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTa                                                                                  >  2:675959/1‑121 (MQ=255)
                                                                                       cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTa                                                                                  <  1:675959/121‑1 (MQ=255)
                                                                                       cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACttt                                                     <  2:640043/150‑1 (MQ=255)
                                                                                          ggTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttctt                                                  <  2:348582/150‑1 (MQ=255)
                                                                                                 cggcggTTGCCGCCAATGGATTAATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgc                                           <  2:777853/150‑1 (MQ=255)
                                                                                                           cgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTc                                 <  1:681288/150‑1 (MQ=255)
                                                                                                            tccAATGTATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                              <  2:864968/87‑1 (MQ=255)
                                                                                                            gccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                              >  1:864968/1‑88 (MQ=255)
                                                                                                            gccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGAACGAGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGGACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCg                                <  2:738832/150‑1 (MQ=255)
                                                                                                             ccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                               <  2:692825/150‑1 (MQ=255)
                                                                                                                    attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTacacac                        >  1:1109740/1‑150 (MQ=255)
                                                                                                                        attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                               <  2:1196960/139‑1 (MQ=255)
                                                                                                                        attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                               >  1:1196960/1‑139 (MQ=255)
                                                                                                                        attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                             <  2:908004/141‑1 (MQ=255)
                                                                                                                        attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTc                    >  1:1148142/1‑150 (MQ=255)
                                                                                                                        attgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTAAGTGGCGCCCCAATTCGTTa                             >  1:908004/1‑141 (MQ=255)
                                                                                                                                   tgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtt         <  2:1148142/150‑1 (MQ=255)
                                                                                                                                      tGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAcc                <  2:704216/140‑1 (MQ=39)
                                                                                                                                      tGGCCTCTTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCAcc                >  1:704216/1‑140 (MQ=39)
                                                                                                                                        gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTATACAAAAACGTAACCATAGGTGCCATTTTtgt                                                              >  2:443001/1‑92 (MQ=255)
                                                                                                                                        gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTATACAAAAACGTAACCATAGGTGCCATTTTtgt                                                              <  1:443001/92‑1 (MQ=255)
                                                                                                                                          ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCa  >  1:1011303/1‑150 (MQ=255)
                                                                                                                                                  aGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                         >  1:292552/1‑71 (MQ=255)
                                                                                                                                                  aGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGt                                                                         <  2:292552/71‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCA  >  NC_002947/4899964‑4900251

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