Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,805,220 T→G 7.2% L152V (TTA→GTA)  tilS → tRNA lysidine(34) synthetase TilS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,805,2200TG 7.2% 150.3 / 2.5 52L152V (TTA→GTA) tilStRNA lysidine(34) synthetase TilS
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base G (2/2);  total (25/27)
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.95e-01

CGGCCCGGGCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTGCAGCAACGTTGGC  >  NC_002947/1805072‑1805364
                                                                                                                                                    |                                                                                                                                                
cGGCCCGGGCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGt                                                                                                                                                  >  3:315782/1‑149 (MQ=255)
    ccGGGCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACt                                                                                                                                              <  2:83508/149‑1 (MQ=255)
     cGGGCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTg                                                                                                                                             <  5:115564/149‑1 (MQ=255)
       ggCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGc                                                                                                                                           <  2:168120/149‑1 (MQ=255)
          gATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGCCCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCGGGTCAGGCCGGTACTGGCCTg                                                                                                                                        >  8:322176/1‑149 (MQ=255)
              tCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCAGTAACGGCCTgacca                                                                                                                                    <  1:302935/149‑5 (MQ=255)
               ccTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTccc                                                                                                                                   >  4:26211/1‑149 (MQ=255)
                cTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCg                                                                                                                                  >  1:293712/1‑149 (MQ=255)
                 tGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGCGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGCCCGTTACTGGCCTGTTCCCgc                                                                                                                                 >  8:321788/1‑148 (MQ=255)
                      aCCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGTGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACa                                                                                                                            >  3:295405/1‑149 (MQ=255)
                      aCCCGCCAGCACCGCGACGACCAGGCCGAGACCCAGCAGTTCCGCCTGTTGCCCCGTGCCGGCCAGCGGGGCCAGGCCGCAAAGCCGGGGCAGCCGGCGCAGCGGCAGGGCAGCCTGGTCAGGCCCTTACAGGCCTGTTGCCGGCaaaa                                                                                                                            >  2:258122/1‑147 (MQ=255)
                          gCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTACGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTg                                                                                                                        >  5:8344/1‑149 (MQ=255)
                            cAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCa                                                                                                                      <  4:189669/149‑1 (MQ=255)
                             aGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGTGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAg                                                                                                                     <  4:295405/149‑1 (MQ=255)
                                          ccAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAg                                                                                                        <  2:293712/149‑1 (MQ=255)
                                           cAGGCCGAGACCGTGCTGTTGCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGATACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGc                                                                                                       <  7:321788/149‑1 (MQ=255)
                                           cAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGc                                                                                                       >  1:166106/1‑149 (MQ=255)
                                                cgtgaCCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAg                                                                                                  <  1:155780/146‑1 (MQ=255)
                                                               ccGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATc                                                                                    >  8:156969/1‑148 (MQ=255)
                                                                 gCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATg                                                                                                              <  3:274374/120‑1 (MQ=255)
                                                                 gCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATg                                                                                                              >  4:274374/1‑120 (MQ=255)
                                                                         cgcgGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGt                                                                          <  7:152486/148‑1 (MQ=255)
                                                                          gcgGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCAGTAACGGCCTgac                                                                                                                                      >  4:72228/1‑85 (MQ=255)
                                                                          gcgGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCAGTAACGGCCTgac                                                                                                                                      <  3:72228/87‑3 (MQ=255)
                                                                                ccGGCCTGCGGGGCCTGGCGGAAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCGTGTTCCCGGCAACATTTGCTGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGt                                                                   <  1:258122/148‑1 (MQ=255)
                                                                                 cGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTc                                                                  >  8:156315/1‑148 (MQ=255)
                                                                                     cTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGaa                                                                                 >  8:204797/1‑129 (MQ=255)
                                                                                     cTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGaa                                                                                 <  7:204797/129‑1 (MQ=255)
                                                                                        cGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGacacac                                                           >  6:40119/1‑148 (MQ=255)
                                                                                            gCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTAGTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTc                                                      <  8:57241/149‑1 (MQ=255)
                                                                                            gCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTc                                                      >  1:115745/1‑149 (MQ=255)
                                                                                              cTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTc                                                      <  5:269808/147‑1 (MQ=255)
                                                                                              cTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTc                                                      >  6:269808/1‑147 (MQ=255)
                                                                                              cTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGc                                                    <  3:26211/149‑1 (MQ=255)
                                                                                               tGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGcc                                                   >  6:320550/1‑149 (MQ=255)
                                                                                                 gcggcAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGccc                                                  <  2:115745/148‑1 (MQ=255)
                                                                                                 gcggcAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGccc                                                  <  5:40119/148‑1 (MQ=255)
                                                                                                    gcAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCaa                                              <  7:156969/149‑1 (MQ=255)
                                                                                                              ggCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTgcgc                                    <  6:8344/149‑1 (MQ=255)
                                                                                                                        gcTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGAt                                                                                     >  7:95710/1‑90 (MQ=255)
                                                                                                                        gcTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGAt                                                                                     <  8:95710/90‑1 (MQ=255)
                                                                                                                          tGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATg                        >  4:248899/1‑149 (MQ=255)
                                                                                                                          tGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATg                        <  4:315782/149‑1 (MQ=255)
                                                                                                                              gcagggcagCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCg                     >  6:175137/1‑148 (MQ=255)
                                                                                                                                    cagCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTg               >  3:197112/1‑148 (MQ=255)
                                                                                                                                     agCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTgca             <  7:322176/149‑1 (MQ=255)
                                                                                                                                     agCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTgca             <  2:166106/149‑1 (MQ=255)
                                                                                                                                        cTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTgcagc           >  3:336406/1‑148 (MQ=255)
                                                                                                                                           gTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGa                                                                                                                   >  6:177841/1‑41 (MQ=255)
                                                                                                                                           gTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGa                                                                                                                   <  5:177841/41‑1 (MQ=255)
                                                                                                                                           gTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTGCAGCAACg       >  6:6549/1‑149 (MQ=255)
                                                                                                                                            tCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTGCAGCAACGt      >  2:201199/1‑149 (MQ=255)
                                                                                                                                                gCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTGCAGCAACGTTGGc  <  7:156315/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGGCCCGGGCGATATCCTGTTCACCGGCCAGCACCGCGACGACCAGGCCGAGACCCTGCTGTTCCGCCTGTTGCGCGGTGCCGGCCTGCGGGGCCTGGCGGCAATGCCGGGGCAGCGGGCGCTGGGGCAGGGCAGCCTGGTCAGGCCGTTACTGGCCTGTTCCCGGCAACATTTGCAGGAGTATGCTCAAGCCCAAGGGCTGACCTGGATCGAAGACCCGTCAAACGTCGACACACAATTCGCCCGCAACTACCTGCGCGGCGAAGTGATGCCGCACTTGCAGCAACGTTGGC  >  NC_002947/1805072‑1805364

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