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,6040CT98.3% 181.4 / ‑6.6 58S13N (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 (28/29);  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 = 9.99e-01

GCGCATCCGCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCACTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAACCATTG  >  NC_000913/443470‑443736
                                                                                                                                      |                                                                                                                                    
gcgcATCCGCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTAcc                                                                                                                                  <  1:343319/139‑1 (MQ=255)
 cgcATCCGCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCa                                                                                                                                 >  2:911935/1‑139 (MQ=255)
        gCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             <  2:1442763/136‑1 (MQ=255)
        gCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             >  1:1442763/1‑136 (MQ=255)
          cagcagCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAgg                                                                                                                        <  1:1071404/139‑1 (MQ=255)
                cTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAAcc                                                                                                                  >  2:453982/1‑139 (MQ=255)
                 ttCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCa                                                                                                                 <  1:1293393/139‑1 (MQ=255)
                    aCGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             <  2:1138135/124‑1 (MQ=255)
                    aCGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             >  1:1138135/1‑124 (MQ=255)
                        cgcgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGa                                                                                                          >  1:1049268/1‑139 (MQ=255)
                          cgcgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGTATAACAATAGTGGTGACGGCTTCAGTCTCTTGATTACCAGGGGCGAGGCAAGCCAGAGCCGATg                                                                                                        <  2:430359/139‑1 (MQ=255)
                            cgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                      <  1:453982/139‑1 (MQ=255)
                            cgAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCACCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                      <  2:1338679/139‑1 (MQ=255)
                              aGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCt                                                                                                    <  2:1039941/139‑1 (MQ=255)
                                             ggTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                     >  1:654150/1‑139 (MQ=255)
                                             ggTTACAATCGCTGCCGACGCTGGCAGTAGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTTGTGACGGCTTCAGTCTCTGCACTACCAGGGGCGAGGCAAACAAGTGCCGATGCGCTCAAATTCACTCGCct                                                                                     <  2:128428/139‑1 (MQ=255)
                                              gTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             >  2:1400560/1‑98 (MQ=255)
                                              gTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAgggg                                                                                                                             <  1:1400560/98‑1 (MQ=255)
                                                   cATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        <  1:1250183/114‑1 (MQ=255)
                                                   cATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        >  2:1250183/1‑114 (MQ=255)
                                                        cTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                     <  2:252146/128‑1 (MQ=255)
                                                        cTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCct                                                                                     >  1:252146/1‑128 (MQ=255)
                                                         tGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCa                                                                         <  2:792578/139‑1 (MQ=255)
                                                            cGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCAttt                                                                      >  1:34899/1‑139 (MQ=255)
                                                               cGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATGCACTCTCCTTTCTTTTTACCATTTCaa                                                                   <  2:1049268/139‑1 (MQ=255)
                                                                          gtgACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGgtg                                                        >  2:140774/1‑139 (MQ=255)
                                                                                      ccgccgCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACg                                            <  2:34899/139‑1 (MQ=255)
                                                                                        gccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGc                                          >  2:730413/1‑139 (MQ=255)
                                                                                        gccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGc                                          >  2:995790/1‑139 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAAc                                                                                                                   >  1:1533639/1‑65 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAAc                                                                                                                   <  2:1533639/65‑1 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                      >  1:1450430/1‑78 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATgcg                                                                                                      <  2:1450430/78‑1 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTa                                                                            >  1:983419/1‑104 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTa                                                                            <  2:983419/104‑1 (MQ=255)
                                                                                         ccgcGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCa                                         >  1:221258/1‑139 (MQ=255)
                                                                                              aaCCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGttt                                                     >  2:1059986/1‑122 (MQ=255)
                                                                                              aaCCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGttt                                                     <  1:1059986/122‑1 (MQ=255)
                                                                                                 cagcagATCGATAGTGGTGACGGCTTCAGTCTCTTCATTACCATGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATg                                 >  2:1429938/1‑139 (MQ=255)
                                                                                                           aTAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCgcgcgggcg                       <  1:140774/139‑1 (MQ=255)
                                                                                                           aTAGTGGTGACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCATTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCgcgcgggcg                       >  2:322805/1‑139 (MQ=255)
                                                                                                               tggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCgcggc                   >  1:1053641/1‑139 (MQ=255)
                                                                                                                ggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        <  2:522936/53‑1 (MQ=255)
                                                                                                                ggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        >  1:522936/1‑53 (MQ=255)
                                                                                                                ggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCgcggcg                  <  1:995790/139‑1 (MQ=255)
                                                                                                                ggtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCgcggcg                  <  1:730413/139‑1 (MQ=255)
                                                                                                                 gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAgg                                                              >  2:1430660/1‑94 (MQ=255)
                                                                                                                 gtgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAgg                                                              <  1:1430660/94‑1 (MQ=255)
                                                                                                                  tgACGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTa                >  1:1286225/1‑139 (MQ=255)
                                                                                                                     cGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        >  2:414563/1‑48 (MQ=255)
                                                                                                                     cGGCTTCAGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATg                                                                                                        <  1:414563/48‑1 (MQ=255)
                                                                                                                            aGTCTCTTCATTACCAGGTGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGACAAcc      >  1:1012251/1‑139 (MQ=255)
                                                                                                                            aGTCTCTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAcc      <  1:911935/139‑1 (MQ=255)
                                                                                                                              tctcTTCATTACCATGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAACCAt    <  1:1429938/139‑1 (MQ=255)
                                                                                                                              tctcTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAcc      >  1:602202/1‑137 (MQ=255)
                                                                                                                              tctcTTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAAcc      <  2:602202/137‑1 (MQ=255)
                                                                                                                                tctTCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAACCATTg  >  1:945084/1‑139 (MQ=255)
                                                                                                                                   tCATTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCa                                                                                            <  1:1494097/46‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GCGCATCCGCCAGCAGCTTCACGCCGCGCGAGCAGGTAATCGCCAGGTTACCATCGCTGGCGACGCTGGCAGTGGTGACTTTGATACCGCCGCGAACCAGCAGATCGATAGTGGTGACGGCTTCAGTCTCTTCACTACCAGGGGCGAGGCAAACCAGTGCCGATGCGCTCATATTCACTCTCCTTTCTTTTTACCATTTCAAACAGGCGGGTGTTTTCCGGTACGGCAATCCCATGCGCGCGGGCGCGGCGTAAGAGAAAACCATTG  >  NC_000913/443470‑443736

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑