Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,940:1 +A 15.4% intergenic (+17/‑45) leuS → / → PP_4795 leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,458,9401.A15.4% 77.2 / 12.3 26intergenic (+17/‑45)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (9/13);  new base A (2/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

AGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458801‑5459088
                                                                                                                                            |                                                                                                                                                    
aGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCg                                                                                                                                            <  2:17573/150‑1 (MQ=255)
 gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCgg                                                                                                                                           <  2:376875/150‑1 (MQ=255)
   cgcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGt                                                                                                                                         >  1:683409/1‑150 (MQ=255)
     cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTgg                                                                                                                                       <  1:753971/150‑1 (MQ=255)
                         cccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                                   <  2:795716/150‑1 (MQ=255)
                                   cGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAacaaca                                                                                                         <  1:702639/150‑1 (MQ=255)
                                                cgcTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGGTGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCaa                                                                                            >  1:386298/1‑150 (MQ=255)
                                                       tCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTg                                                                                     <  2:683409/150‑1 (MQ=255)
                                                                    cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                        <  1:174263/150‑1 (MQ=255)
                                                                     aTCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGc                                                                       >  1:266070/1‑150 (MQ=255)
                                                                      tCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTGGGc                                                                                                                                     >  1:252436/1‑87 (MQ=38)
                                                                      tCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCAATCGATGCGGGTGGGc                                                                                                                                     <  2:252436/87‑1 (MQ=38)
                                                                           aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGc                                                                                                                                     >  1:523347/1‑82 (MQ=255)
                                                                           aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGc                                                                                                                                     <  2:523347/82‑1 (MQ=255)
                                                                           aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCAt                                                                 <  2:787500/150‑1 (MQ=255)
                                                                                  tCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                          <  2:266070/150‑1 (MQ=255)
                                                                                          ccGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg                                                  <  1:387827/150‑1 (MQ=255)
                                                                                                                           tgaAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCt                 <  1:680457/150‑1 (MQ=255)
                                                                                                                             aaaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGa               >  2:679379/1‑150 (MQ=255)
                                                                                                                              aaaCGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAg              >  1:274634/1‑150 (MQ=255)
                                                                                                                                 cGCCCACCCGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGt           >  2:53313/1‑150 (MQ=255)
                                                                                                                                    ccACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGaa        <  1:679379/150‑1 (MQ=255)
                                                                                                                                        ccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGa         >  1:85349/1‑145 (MQ=255)
                                                                                                                                        ccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGa         <  2:85349/145‑1 (MQ=255)
                                                                                                                                        ccGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGaa    >  2:490961/1‑150 (MQ=255)
                                                                                                                                          gCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAt  >  2:467329/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
AGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458801‑5459088

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