Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 443,604 C→T S13N (AGT→AAT) 
*304* (TAG→TAA
yajL ←
panE ←
oxidative‑stress‑resistance chaperone
2‑dehydropantoate reductase, NADPH‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,6040CT97.2% 105.3 / ‑4.0 36S13N (AGT→AAT) 
*304* (TAG→TAA
yajL
panE
oxidative‑stress‑resistance chaperone
2‑dehydropantoate reductase, NADPH‑specific
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCACTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAC  >  NC_000913/443485‑443731
                                                                                                                       |                                                                                                                               
gcTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAAc                                                                                                              <  1:57804/139‑1 (MQ=255)
         cgcgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCACTACCAGGGGCGAGGCAAACCAGTGCCGa                                                                                                     <  1:561841/139‑1 (MQ=255)
          gcgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGAt                                                                                                    <  2:685945/139‑1 (MQ=255)
           cgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   <  1:198780/139‑1 (MQ=255)
           cgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   <  2:147258/139‑1 (MQ=255)
                 cAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCa                                                                                             <  2:164497/139‑1 (MQ=255)
                             aGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   >  2:454586/1‑121 (MQ=255)
                             aGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   <  1:454586/121‑1 (MQ=255)
                               gTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCctt                                                                               >  1:393899/1‑139 (MQ=255)
                                         cTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                 >  1:551435/1‑111 (MQ=255)
                                         cTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                 <  2:551435/111‑1 (MQ=255)
                                             cGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCAttt                                                                 >  2:459471/1‑139 (MQ=255)
                                                 gCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCaaa                                                             >  2:268676/1‑139 (MQ=255)
                                                         tggtgACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                <  1:393036/112‑1 (MQ=255)
                                                         tggtgACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                >  2:393036/1‑112 (MQ=255)
                                                                tttGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTc                                              >  2:281999/1‑139 (MQ=255)
                                                                   gATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGTCGATGCGCTCATATTc                                                                                        <  2:321669/94‑1 (MQ=255)
                                                                   gATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGTCGATGCGCTCATATTc                                                                                        >  1:321669/1‑94 (MQ=255)
                                                                   gATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCgg                                           >  2:751463/1‑139 (MQ=255)
                                                                    aTACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGt                                          >  1:759655/1‑139 (MQ=255)
                                                                         gccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGc                                     <  1:268676/139‑1 (MQ=255)
                                                                          ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   <  2:964440/76‑1 (MQ=255)
                                                                          ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                   >  1:964440/1‑76 (MQ=255)
                                                                               aaCCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCGCTCTCCTTTCTATTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATccc                               <  2:759655/139‑1 (MQ=255)
                                                                                            aTAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCgcgcgggcg                  <  1:751463/139‑1 (MQ=255)
                                                                                               gtggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATg                            <  2:199908/126‑1 (MQ=255)
                                                                                               gtggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATg                            >  1:199908/1‑126 (MQ=255)
                                                                                                  gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                >  2:1058515/1‑71 (MQ=255)
                                                                                                  gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                <  1:1058515/71‑1 (MQ=255)
                                                                                                  gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCgcgcgggcg                  <  2:686630/133‑1 (MQ=255)
                                                                                                  gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCgcgcgggcg                  >  1:686630/1‑133 (MQ=255)
                                                                                                  gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGt            >  2:185209/1‑139 (MQ=255)
                                                                                                       ggCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAagag       >  2:1002123/1‑139 (MQ=255)
                                                                                                           tCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAc  >  1:1061505/1‑139 (MQ=255)
                                                                                                                   ttCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGaa     <  2:499437/129‑1 (MQ=255)
                                                                                                                   ttCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGaa     >  1:499437/1‑129 (MQ=255)
                                                                                                                       |                                                                                                                               
GCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCACTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAC  >  NC_000913/443485‑443731

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