Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,405,453 C→T P251P (CCG→CCA ttcA ← tRNA 2‑thiocytidine(32) synthetase TtcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,405,4530CT84.8% 122.0 / 12.4 46P251P (CCG→CCAttcAtRNA 2‑thiocytidine(32) synthetase TtcA
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base T (20/19);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

CTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATC  >  NZ_CP009273/1405320‑1405588
                                                                                                                                     |                                                                                                                                       
cttcttctGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGcac                                                                                                                                    >  2:499390/1‑139 (MQ=255)
   cttctGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCGCCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGcacc                                                                                                                                   <  2:117191/137‑1 (MQ=255)
             ccAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCa                                                                                                                         <  1:1025015/137‑1 (MQ=255)
             ccAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATc                                                                                                                       >  2:888654/1‑139 (MQ=255)
                                 ggATCTCTTCGCGATCAAACGCTAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGAt                                                                                                   <  1:1178629/139‑1 (MQ=255)
                                 ggATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGAt                                                                                                   <  1:132997/139‑1 (MQ=255)
                                      tctTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGcc                                                                                              >  2:414404/1‑139 (MQ=255)
                                         tCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGATATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTg                                                                                           >  2:25757/1‑139 (MQ=255)
                                          cgcgATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTg                                                                                           >  1:1070438/1‑138 (MQ=255)
                                                                  cGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCaa                                                                  <  1:64946/139‑1 (MQ=255)
                                                                       aCCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGt                                                             <  2:1035264/139‑1 (MQ=255)
                                                                        ccACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGca                                                                                                                                     <  2:131231/66‑1 (MQ=255)
                                                                        ccACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGca                                                                                                                                     >  1:131231/1‑66 (MQ=255)
                                                                                  aCCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCAcc                                                  <  1:876818/139‑1 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGcacc                                                                                                                                   <  2:829658/54‑1 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGcacc                                                                                                                                   >  1:829658/1‑54 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGcc                                                                                              >  1:515147/1‑91 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGcc                                                                                              <  2:515147/91‑1 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGtt                                                                                  >  2:239847/1‑103 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGtt                                                                                  <  1:239847/103‑1 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGAc                                              >  2:1275177/1‑139 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGAc                                              >  2:776244/1‑139 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCa                                                        >  2:1140467/1‑129 (MQ=255)
                                                                                      tGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCa                                                        <  1:1140467/129‑1 (MQ=255)
                                                                                         gTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGtt                                           <  2:1011965/139‑1 (MQ=255)
                                                                                             tGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCa                                                                       <  2:1245537/107‑1 (MQ=255)
                                                                                             tGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCa                                                                       >  1:1245537/1‑107 (MQ=255)
                                                                                               ccTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAgg                                      <  2:938662/138‑1 (MQ=255)
                                                                                                        tCGAACAGGTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCtgtg                            <  2:1070438/139‑1 (MQ=255)
                                                                                                                gTTGGTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGcac                    <  1:291900/139‑1 (MQ=255)
                                                                                                                    gTATCGCACAGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAgg                 <  1:888654/138‑1 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGt                                                                                                         >  2:712095/1‑44 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGt                                                                                                         <  1:712095/44‑1 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGAt                                                                                        >  1:921508/1‑61 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGAt                                                                                        <  2:921508/61‑1 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGc                                                                    <  1:476818/81‑1 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGc                                                                    >  2:476818/1‑81 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACg          >  1:319822/1‑139 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACg          >  2:192051/1‑139 (MQ=255)
                                                                                                                          cacaGATGCGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACg          >  2:1080095/1‑139 (MQ=255)
                                                                                                                          cacaGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCa            <  1:1057682/137‑1 (MQ=255)
                                                                                                                          cacaGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCa            >  2:1057682/1‑137 (MQ=255)
                                                                                                                               aTGCGATGGCACCACATTCTGCATCGTGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGaat      >  2:9720/1‑138 (MQ=255)
                                                                                                                               aTGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGAt                                                                                                   <  2:1245701/45‑1 (MQ=255)
                                                                                                                               aTGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGAt                                                                                                   >  1:1245701/1‑45 (MQ=255)
                                                                                                                                  cGATGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATc  >  1:1278048/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CTTCTTCTGGCTGCCAGCACGCCGGTTGTAGTGGGATCTCTTCGCGATCAAACGCCAGATCACCCCCGTTAACCACTTCAGAACCGTGGGTGATGCCTTTGAAATCGAACAGGTTGGTATCGCACAGATGCGACGGCACCACATTCTGCATCGCGCTGAACATCGTCTCGATACGCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATC  >  NZ_CP009273/1405320‑1405588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: