Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,513,723 C→T 86.2% P336L (CCC→CTC)  mraY → Phospho‑N‑acetylmuramoyl‑pentapeptide‑transferas e

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,513,7230CT86.2% 128.4 / 12.5 51P336L (CCC→CTC) mraYPhospho‑N‑acetylmuramoyl‑pentapeptide‑transferas e
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base T (24/19);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCATCTTCGTGGTGGAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCGGCAAGAG  >  ACB143/1513600‑1513858
                                                                                                                           |                                                                                                                                       
gcATCTTCGTGGTGGAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGCCCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACATTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGt                                                                                                                             >  2:153303/1‑136 (MQ=255)
  aTCTTCGTGGTGGAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGTCGCCGATTCATCACCACTTTGAACTCAATGGCTGGCCAGAGCTCAGGGTGATCGTcc                                                                                                                           >  1:9612/1‑136 (MQ=255)
     ttCGTGGTGGAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGtt                                                                                                                        >  4:215803/1‑136 (MQ=255)
            tggAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGat                                                                                                                 <  6:70893/136‑1 (MQ=255)
              gAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGatca                                                                                                               >  7:185894/1‑136 (MQ=255)
                  cccTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCAc                                                                                                           <  3:286239/136‑1 (MQ=255)
                  cccTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCAc                                                                                                           <  1:237020/136‑1 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATc                                                                                                    >  7:233413/1‑136 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATc                                                                                                    >  2:166526/1‑136 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGGGATc                                                                                                    >  8:191129/1‑136 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGAGCGTCCGTTTCTGGATCATAACCGTGATc                                                                                                    >  5:161835/1‑136 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGAGCGTGTTCCGCATGGCGCCGATTCATCACCACGTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGTTGGTCCGTTTCTGGATCATCACCGTGATc                                                                                                    >  8:195076/1‑136 (MQ=255)
                         ggtggtGATCCAGGTCGCCTCCTTCAAGCTGACCGACAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAATGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATc                                                                                                    >  7:113061/1‑136 (MQ=255)
                                aTCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTAc                                                                                             >  5:111200/1‑136 (MQ=255)
                                   cAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGa                                                                                          >  2:195290/1‑136 (MQ=255)
                                   cAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGCGCTGGCCAGAGCTCAGGGTGAATGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGa                                                                                          >  6:231814/1‑136 (MQ=255)
                                    aGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGAt                                                                                         <  7:37295/136‑1 (MQ=255)
                                        cGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGc                                                                                     >  3:58993/1‑136 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCGGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCGAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  1:28428/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  5:115985/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  4:30969/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  2:252986/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  1:121422/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  3:39986/136‑1 (MQ=255)
                                          cctcctTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCt                                                                                   <  1:279400/136‑1 (MQ=255)
                                                 cAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAAc                                                                             <  6:111200/135‑1 (MQ=255)
                                                     cTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCGAGAGCGCAGTGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGa                                                                        <  8:13809/136‑1 (MQ=255)
                                                           ggCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGTGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCAc                                                                                                           <  3:264511/95‑1 (MQ=255)
                                                           ggCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCAc                                                                                                           >  4:264511/1‑95 (MQ=255)
                                                           ggCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACt                                                                    >  2:137921/1‑134 (MQ=255)
                                                           ggCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACt                                                                    >  6:44523/1‑134 (MQ=255)
                                                           ggCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCATGAAACt                                                                    >  6:254640/1‑134 (MQ=255)
                                                                              tGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCttt                                               <  4:107019/136‑1 (MQ=255)
                                                                              tGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCttt                                               <  6:56775/136‑1 (MQ=255)
                                                                                  gcCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATc                                           >  6:203605/1‑136 (MQ=255)
                                                                                   ccGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCg                                          <  5:162513/136‑1 (MQ=255)
                                                                                         catcaCCACTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTcc                                     >  1:74914/1‑135 (MQ=255)
                                                                                             accacTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATc                                                                                                                               <  2:197620/41‑1 (MQ=255)
                                                                                             accacTTTGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATc                                                                                                                               >  1:197620/1‑41 (MQ=255)
                                                                                                    tGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTcc                           <  6:256299/134‑1 (MQ=255)
                                                                                                    tGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTcc                           <  5:44523/134‑1 (MQ=255)
                                                                                                    tGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTcc                           <  2:177388/134‑1 (MQ=255)
                                                                                                    tGAACTCAAGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCg                          <  3:59155/135‑1 (MQ=255)
                                                                                                     gAACTCAAGGGATGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCa                        >  6:224924/1‑136 (MQ=255)
                                                                                                           aaGGGCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGt                    <  5:155439/134‑1 (MQ=255)
                                                                                                              ggCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCACCCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCg               >  2:49613/1‑136 (MQ=255)
                                                                                                               gCTGGCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCgg              >  5:59349/1‑136 (MQ=255)
                                                                                                                   gCCAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTc          >  6:194493/1‑136 (MQ=255)
                                                                                                                    ccAGAGCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCg         <  7:43107/136‑1 (MQ=255)
                                                                                                                      agagCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCgg        >  8:231486/1‑135 (MQ=255)
                                                                                                                        agCTCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCGGCaa     <  1:153303/136‑1 (MQ=255)
                                                                                                                           ccAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCGGCAagag  >  1:230715/1‑136 (MQ=255)
                                                                                                                           |                                                                                                                                       
GCATCTTCGTGGTGGAAACCCTGTCGGTGGTGATCCAGGTCGCCTCCTTCAAGCTGACCGGCAAGCGCGTGTTCCGCATGGCGCCGATTCATCACCACTTTGAACTCAAGGGCTGGCCAGAGCCCAGGGTGATCGTCCGTTTCTGGATCATCACCGTGATCCTGGTACTGATTGGCCTTGCAACCCTGAAACTGAGGTAGATGAGCGTGTCTTTGATCGCTTCCGACCAATTCCGCATCGTTGTCGGCCTCGGCAAGAG  >  ACB143/1513600‑1513858

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