Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 209,716 T→C I300T (ATC→ACC)  tilS → tRNA lysidine(34) synthetase TilS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273209,7160TC98.1% 164.8 / ‑3.7 52I300T (ATC→ACC) tilStRNA lysidine(34) synthetase TilS
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/28);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

TTCCCGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCA  >  NZ_CP009273/209627‑209800
                                                                                         |                                                                                    
ttCCCGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAAt                                                                                      >  2:754453/1‑90 (MQ=255)
     gcgACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGAc                                                                                 <  1:946311/90‑1 (MQ=255)
            gTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCa                                                                          <  1:721702/90‑1 (MQ=255)
               ggTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTc                                                                       <  2:315607/90‑1 (MQ=255)
                    ggATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAAc                                                                  <  1:1114121/90‑1 (MQ=255)
                    ggATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAAc                                                                  <  2:628159/90‑1 (MQ=255)
                         tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  2:1576155/1‑90 (MQ=255)
                         tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  2:712387/1‑90 (MQ=255)
                         tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  1:389166/1‑90 (MQ=255)
                         tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  1:375616/1‑90 (MQ=255)
                          ggCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                            <  2:1002985/90‑1 (MQ=255)
                          ggCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                            <  2:389166/90‑1 (MQ=255)
                          ggCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                            <  2:1308912/90‑1 (MQ=255)
                                        ctggcgcGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTc                                              >  2:1457309/1‑90 (MQ=255)
                                                ggAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGTGCa                                      <  2:1224315/90‑1 (MQ=255)
                                                  aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  1:1576155/90‑1 (MQ=255)
                                                  aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  2:739669/90‑1 (MQ=255)
                                                  aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  2:1444698/90‑1 (MQ=255)
                                                  aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  2:52522/90‑1 (MQ=255)
                                                    gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGAc                                                                                 <  1:418757/43‑1 (MQ=255)
                                                    gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGAc                                                                                 >  2:418757/1‑43 (MQ=255)
                                                    gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                            >  1:238308/1‑64 (MQ=255)
                                                    gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                            <  2:238308/64‑1 (MQ=255)
                                                          tCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCa                                      >  2:1310189/1‑80 (MQ=255)
                                                          tCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCa                                      <  1:1310189/80‑1 (MQ=255)
                                                          tCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAAc                            >  1:1179309/1‑90 (MQ=255)
                                                           cACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACa                           >  1:431258/1‑90 (MQ=255)
                                                                 gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATcc                                                 >  1:784534/1‑62 (MQ=255)
                                                                 gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATcc                                                 <  2:784534/62‑1 (MQ=255)
                                                                 gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGc                     <  2:657435/90‑1 (MQ=255)
                                                                      cGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTgg                <  1:857654/90‑1 (MQ=255)
                                                                       gTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGc               >  1:1315846/1‑90 (MQ=255)
                                                                        tttGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCa              <  1:1397974/90‑1 (MQ=255)
                                                                        tttGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCa              <  2:848675/90‑1 (MQ=255)
                                                                         ttGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAg             >  1:547185/1‑90 (MQ=255)
                                                                         ttGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAg             <  2:1315846/90‑1 (MQ=255)
                                                                           gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  2:1137960/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTa                                                      <  2:1303077/47‑1 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTa                                                      >  1:1303077/1‑47 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCAcc                                           <  2:523649/58‑1 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCAcc                                           >  1:523649/1‑58 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  2:327592/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:846817/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:37089/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:169705/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:1341183/1‑90 (MQ=255)
                                                                           gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:1268677/1‑90 (MQ=255)
                                                                            ggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACg          <  1:1457309/90‑1 (MQ=255)
                                                                               gcgTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               >  1:879688/1‑66 (MQ=255)
                                                                               gcgTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               <  2:879688/66‑1 (MQ=255)
                                                                                cgTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGc      >  2:1363067/1‑90 (MQ=255)
                                                                                    tGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCa  <  2:811573/90‑1 (MQ=255)
                                                                                    tGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCa  <  2:933808/90‑1 (MQ=255)
                                                                                         |                                                                                    
TTCCCGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCA  >  NZ_CP009273/209627‑209800

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