Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,441,742 G→A intergenic (‑137/‑22) fabB ← / → mnmC 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,441,7420GA94.7% 110.5 / ‑0.4 38intergenic (‑137/‑22)fabB/mnmC3‑oxoacyl‑[acyl‑carrier‑protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/15);  minor base . (1/1);  total (22/16)
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.00e-01

CCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAG  >  NC_000913/2441612‑2441874
                                                                                                                                  |                                                                                                                                    
ccTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAg                                                                                                                              <  2:1177014/139‑1 (MQ=255)
   ctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACg                                                                                                                           >  1:1093436/1‑139 (MQ=255)
      tAAGTCGCACATAGAGTAAGTTTCGAATGCATAATAGCGTGCACTGGTACGGTGAACAAGTCCGATCAGCCACTTAATAGAGAAATTTGCGCAGCGTTACACACATCGCTAAGAGCGATCCACCACCTGTAAGACGAGt                                                                                                                        >  1:1100906/1‑139 (MQ=255)
      tAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGt                                                                                                                        >  2:794850/1‑139 (MQ=255)
                aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTCACTTACgcg                                                                                                              >  1:805954/1‑137 (MQ=255)
                aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTg                                                                                                              <  2:96490/139‑1 (MQ=255)
                 tAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGa                                                                                                             <  2:1093436/139‑1 (MQ=255)
                   gagTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGaaa                                                                                                           <  1:510182/139‑1 (MQ=255)
                              gAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACactac                                                                                                     >  2:739966/1‑134 (MQ=255)
                              gAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACactac                                                                                                     <  1:739966/134‑1 (MQ=255)
                              gAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCAt                                                                                                >  2:1339521/1‑139 (MQ=255)
                                             tACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCt                                                                                 >  2:1066948/1‑139 (MQ=255)
                                                    gTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACacta                                                                                                      >  2:1274579/1‑111 (MQ=255)
                                                    gTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACacta                                                                                                      <  1:1274579/111‑1 (MQ=255)
                                                    gTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAAttt                                                                          >  2:625911/1‑139 (MQ=255)
                                                       cGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCA‑‑‑CCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCt                                                                    >  2:192224/1‑139 (MQ=255)
                                                          cGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCt                                                                    <  1:775297/139‑1 (MQ=255)
                                                             aCAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCA‑‑‑CCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGt                                                              <  1:192224/139‑1 (MQ=255)
                                                                    cGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCa                                                                                     >  2:566625/1‑112 (MQ=255)
                                                                    cGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCa                                                                                     <  1:566625/112‑1 (MQ=255)
                                                                        cAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGt                                                      >  1:436795/1‑139 (MQ=255)
                                                                                 aaTAGAGAAATTTGCGCAGCCCTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTTCGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCCGTTTcctggga                                             <  1:794850/139‑6 (MQ=255)
                                                                                     gagaAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCGGAATtttt                                         >  1:1329194/1‑139 (MQ=255)
                                                                                                              aTCGCTAATATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACacta                                                                                                      <  1:966507/53‑1 (MQ=255)
                                                                                                              aTCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACacta                                                                                                      >  2:966507/1‑53 (MQ=255)
                                                                                                              aTCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATa                >  2:200721/1‑139 (MQ=255)
                                                                                                              aTCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATa                >  1:314567/1‑139 (MQ=255)
                                                                                                              aTCGCTAAGATCGAGCCACCACCTGTAAAACGAGAAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATa                >  1:885767/1‑139 (MQ=255)
                                                                                                                cGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATg                                  <  2:1132503/119‑1 (MQ=255)
                                                                                                                cGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATg                                  >  1:1132503/1‑119 (MQ=255)
                                                                                                                    aaGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGc          <  2:436795/139‑1 (MQ=255)
                                                                                                                        tCGAGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGa      >  1:382724/1‑139 (MQ=255)
                                                                                                                           aGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGa                                                                              >  1:1241720/1‑64 (MQ=255)
                                                                                                                           aGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGa                                                                              <  2:1241720/64‑1 (MQ=255)
                                                                                                                           aGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAaga   >  1:1081122/1‑139 (MQ=255)
                                                                                                                           aGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAaga   >  1:244766/1‑139 (MQ=255)
                                                                                                                           aGCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAaga   <  1:1339521/139‑1 (MQ=255)
                                                                                                                            gCCACCACCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAagag  <  1:200721/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTGAGGGTACACCTGTTTCCCGAGATTTTGACGATGTCTATTTTTCCAACGATAACGGGCTGGAAGAG  >  NC_000913/2441612‑2441874

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