Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 383,259 C→T R226Q (CGG→CAG)  yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,2590CT98.2% 178.8 / NA 57R226Q (CGG→CAG) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/29);  minor base G (0/1);  total (27/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

GTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAAT  >  NC_000913/383143‑383383
                                                                                                                    |                                                                                                                            
gTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                        >  1:85309/1‑139 (MQ=255)
 tCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTATTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacc                                                                                                       <  1:1051098/139‑1 (MQ=255)
  cTTCTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                      <  1:1050434/139‑1 (MQ=255)
   ttttATGTAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCATCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCaa                                                                                                     >  1:1057122/1‑139 (MQ=255)
       aTGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                 <  1:854054/139‑1 (MQ=255)
           aGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                             <  2:85309/139‑1 (MQ=255)
             gcgcATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTg                                                                                           >  2:925583/1‑139 (MQ=255)
                       aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                 >  2:90877/1‑139 (MQ=255)
                       aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                 >  2:730951/1‑139 (MQ=255)
                          gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                              <  1:90877/139‑1 (MQ=255)
                           cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                             >  1:517682/1‑139 (MQ=255)
                           cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                             <  1:4046/139‑1 (MQ=255)
                             tATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGa                                                                           <  1:925583/139‑1 (MQ=255)
                               tCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGCATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                         <  1:1019019/139‑1 (MQ=255)
                                     aaTCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCTAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                   <  2:1057122/139‑1 (MQ=255)
                                     aaTCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCt                                                                               <  1:1061696/127‑1 (MQ=255)
                                     aaTCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCt                                                                               >  2:1061696/1‑127 (MQ=255)
                                     aaTCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                   <  1:1378646/139‑1 (MQ=255)
                                          aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATc                                                              >  1:844523/1‑139 (MQ=255)
                                               cgccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGTGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCgaga                                                         >  2:634372/1‑139 (MQ=255)
                                               cgccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACTATCCgaga                                                         >  1:490711/1‑139 (MQ=255)
                                                  cAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAAc                                                      <  2:844523/139‑1 (MQ=255)
                                                       aCGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                 <  2:464489/139‑1 (MQ=255)
                                                          ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                              <  2:331958/139‑1 (MQ=255)
                                                            aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGTGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                            <  1:634372/139‑1 (MQ=255)
                                                              tcctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                            <  1:359363/137‑1 (MQ=255)
                                                              tcctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                            >  2:359363/1‑137 (MQ=255)
                                                                    cgcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                    <  1:1471182/139‑1 (MQ=255)
                                                                                     cTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCg                   >  2:1261616/1‑139 (MQ=255)
                                                                                      tgtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGc                  <  2:490711/139‑1 (MQ=255)
                                                                                       gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCa                 >  2:791221/1‑139 (MQ=255)
                                                                                             ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                      <  2:1009152/48‑1 (MQ=255)
                                                                                             ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                      >  1:1009152/1‑48 (MQ=255)
                                                                                               aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                        <  2:1233592/44‑1 (MQ=255)
                                                                                               aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                        >  1:1233592/1‑44 (MQ=255)
                                                                                               aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCGGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAGATCAGTACGGAACATCCCGCAGGCACCGc         >  2:1090644/1‑139 (MQ=255)
                                                                                               aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                             <  1:1034257/55‑1 (MQ=255)
                                                                                               aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                             >  2:1034257/1‑55 (MQ=255)
                                                                                                 gATGCAGCGATTCGCCTGCGGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg       <  1:1215860/139‑1 (MQ=255)
                                                                                                   tGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGat     >  2:1547734/1‑139 (MQ=255)
                                                                                                    gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGata    >  2:779113/1‑139 (MQ=255)
                                                                                                     cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                             <  2:1179506/65‑1 (MQ=255)
                                                                                                     cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                             >  1:1179506/1‑65 (MQ=255)
                                                                                                      aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg       <  2:558197/134‑1 (MQ=255)
                                                                                                      aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg       >  1:558197/1‑134 (MQ=255)
                                                                                                        cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                   <  1:662168/72‑1 (MQ=255)
                                                                                                        cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                   >  2:662168/1‑72 (MQ=255)
                                                                                                        cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                         >  1:1502280/1‑98 (MQ=255)
                                                                                                        cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                         <  2:1502280/98‑1 (MQ=255)
                                                                                                          aTTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCt                                                                               <  1:369245/58‑1 (MQ=255)
                                                                                                          aTTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGATGAGATCAAGGTCt                                                                               >  2:369245/1‑58 (MQ=255)
                                                                                                             cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                 <  1:292564/85‑1 (MQ=255)
                                                                                                             cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                 >  2:292564/1‑85 (MQ=255)
                                                                                                             cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataat  >  1:617163/1‑132 (MQ=255)
                                                                                                             cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataat  <  2:617163/132‑1 (MQ=255)
                                                                                                                 tGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                    >  2:357781/1‑94 (MQ=255)
                                                                                                                 tGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                    <  1:357781/94‑1 (MQ=255)
                                                                                                                    |                                                                                                                            
GTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAAT  >  NC_000913/383143‑383383

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: