Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,548 G→A 6.9% R246C (CGT→TGT)  feaR ← transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,5480GA 6.9% 179.8 / 2.6 58R246C (CGT→TGT) feaRtranscriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACT  >  NC_000913/1446402‑1446694
                                                                                                                                                  |                                                                                                                                                  
gCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGaa                                                                                                                                                  <  8:32542/149‑1 (MQ=255)
gCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGaa                                                                                                                                                  >  7:32542/1‑149 (MQ=255)
gCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGaa                                                                                                                                                  >  4:335582/1‑149 (MQ=255)
gCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGaa                                                                                                                                                  <  3:335582/149‑1 (MQ=255)
     gCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAatat                                                                                                                                               <  1:118754/147‑1 (MQ=255)
             ccAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGCCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGAc                                                                                                                                     >  5:321726/1‑149 (MQ=255)
              cAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGCCGAATCGCATCTGCACAAAAATCGAGACGACGCTTACGAATATATTGCGCGAc                                                                                                                                     >  6:145701/1‑148 (MQ=255)
                         ttAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCttt                                                                                                                         >  5:329125/1‑149 (MQ=255)
                         ttAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCtt                                                                                                                          >  2:112517/1‑148 (MQ=255)
                          tAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTa                                                                                                                        >  5:204014/1‑149 (MQ=255)
                           aaaTACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTa                                                                                                                        <  2:238976/148‑1 (MQ=255)
                            aaTACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATc                                                                                                                      <  6:204014/149‑1 (MQ=255)
                                ccGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAatat                                                                                                                                               <  6:365655/120‑1 (MQ=255)
                                ccGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAatat                                                                                                                                               >  5:365655/1‑120 (MQ=255)
                                ccGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGc                                                                                                                   <  4:341881/148‑1 (MQ=255)
                                       aaaaTGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTc                                                                                                           <  2:357015/149‑1 (MQ=255)
                                               tctGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTa                                                                                                       >  6:360881/1‑145 (MQ=255)
                                               tctGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTa                                                                                                       <  5:360881/145‑1 (MQ=255)
                                                 tGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACaaa                                                                                                   <  8:298219/147‑1 (MQ=255)
                                                   gTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAAc                                                                                                  <  6:329125/146‑1 (MQ=255)
                                                   gTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACt                                                                                                 <  1:112517/147‑1 (MQ=255)
                                                   gTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGACACGACGGTTACGAATATATTGCGCGCCTACCAAACCTTTATCGGCAAACATTCGGTACAAACTAc                                                                                               <  5:145701/149‑1 (MQ=255)
                                                    tCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACg                                                                                              >  3:143883/1‑149 (MQ=255)
                                                        aaaaTCCCCAATGAAAGCCGCTGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCTAGACGACGGATACGAATATATTGCGCGACTCCCAAACCTTTATCGGCAAACATTCGGTACAAACtacttact                                                                                          >  6:234353/1‑149 (MQ=255)
                                                         aaaTACCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTg                                                                                         <  3:319231/149‑1 (MQ=255)
                                                          aaTCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGa                                                                                        <  8:131725/149‑1 (MQ=255)
                                                              cccAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGCCTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATa                                                                                    >  6:351180/1‑149 (MQ=255)
                                                              cccAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGa                                                                                        <  2:293854/145‑1 (MQ=255)
                                                              cccAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGa                                                                                        >  1:293854/1‑145 (MQ=255)
                                                                cAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAAc                                                                                                                             >  8:278534/1‑106 (MQ=255)
                                                                cAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAAc                                                                                                                             <  7:278534/106‑1 (MQ=255)
                                                                cAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATAcc                                                                                  >  5:133671/1‑149 (MQ=255)
                                                                 aaTGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACAt                                                                                     >  1:326703/1‑145 (MQ=255)
                                                                 aaTGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACAt                                                                                     <  2:326703/145‑1 (MQ=255)
                                                                    gAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGtc                                                                              >  1:296380/1‑149 (MQ=255)
                                                                    gAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACAAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGtc                                                                              >  5:208297/1‑149 (MQ=255)
                                                                    gAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACAAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGtc                                                                              >  5:202867/1‑149 (MQ=255)
                                                                                ctgccAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATc                                                                  <  6:321726/149‑1 (MQ=255)
                                                                                        ttttCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCgg                                                          >  2:185973/1‑149 (MQ=255)
                                                                                                     cGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTATTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATctct                                             <  6:309716/149‑1 (MQ=255)
                                                                                                            gCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTGGCGa                                       <  6:133671/148‑1 (MQ=255)
                                                                                                              gAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAata                                    >  4:149862/1‑149 (MQ=255)
                                                                                                                        tGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTa                          >  5:49847/1‑149 (MQ=255)
                                                                                                                          cacaAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTc                            <  4:143883/145‑1 (MQ=255)
                                                                                                                          cacaAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCt                           <  3:149862/146‑1 (MQ=255)
                                                                                                                          cacaAAAATCGAGACGACGGTTACAAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCt                           <  6:202867/146‑1 (MQ=255)
                                                                                                                          cacaAAAATCGAGACGACGGTTACAAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCt                           <  6:208297/146‑1 (MQ=255)
                                                                                                                           acaAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATc                        >  1:58505/1‑148 (MQ=255)
                                                                                                                             aaaaaTCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATCCCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAAc                     <  4:238802/149‑1 (MQ=255)
                                                                                                                             aaaaaTCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTa                          <  6:78836/144‑1 (MQ=255)
                                                                                                                             aaaaaTCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTa                          >  5:78836/1‑144 (MQ=255)
                                                                                                                                      gacgacGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACt            >  7:93831/1‑149 (MQ=255)
                                                                                                                                         gacgGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACtttt         >  1:90225/1‑149 (MQ=255)
                                                                                                                                         gacgGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACtttt         <  6:49847/149‑1 (MQ=255)
                                                                                                                                          acgGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACttttt        <  1:185973/149‑1 (MQ=255)
                                                                                                                                              ttACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGaaa    >  4:7726/1‑149 (MQ=255)
                                                                                                                                                 cGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACt  >  4:179806/1‑148 (MQ=255)
                                                                                                                                                 cGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACt  <  2:58505/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCGCGACTACCAAACCTTTATCGGCAAACATTCGGTACAAACTACGTACTGACATACCTGTCTCTCCGGCTATCCACTCCGGGCGTAATATCTCTTCGCGAATATTATCGTCTATCAACGTAACCACTTTTTGAAACT  >  NC_000913/1446402‑1446694

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