Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,030,449 A→G 8.7% F74L (TTT→CTT)  tusE ← mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,030,4490AG 8.7% 121.2 / 3.5 46F74L (TTT→CTT) tusEmnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base G (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

ATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCAATCACCACAGCCAGCGGCTCA  >  NC_000913/1030302‑1030597
                                                                                                                                                   |                                                                                                                                                    
atcatcCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAAt                                                                                                                                                     <  2:553251/149‑1 (MQ=255)
      ccTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCTAATTTATTc                                                                                                                                               >  2:229757/1‑149 (MQ=255)
       cTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCg                                                                                                                                              >  1:415574/1‑149 (MQ=255)
       cTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCg                                                                                                                                              >  2:1008232/1‑149 (MQ=255)
         cGGGGGAATTAACCGATTTTCGTATCCGTTAAATACAGTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAAGTGCGGGTGTTACCTTTTTCTTCGCCAAAATTATTCgac                                                                                                                                            >  2:353563/1‑147 (MQ=255)
            gggAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTAGCTTTTTCTTCGCCAAATTTATTCGCCATc                                                                                                                                         >  2:188833/1‑149 (MQ=255)
                 ttAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCttt                                                                                                                                    >  2:938052/1‑149 (MQ=255)
                     ccGATTTTCGTTTCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGCCCTTTCGGATACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTTcc                                                                                                                                >  2:425641/1‑149 (MQ=255)
                         ttttCGAATCCGTTAAATACATTTTACCGGTTTTGGCCGGCCAGCAATTATGGTGGCTTGCATAGCCGGACCTTTAGGAAACAGACGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCAATCGCATTTCCCAGCa                                                                                                                            >  2:427958/1‑149 (MQ=255)
                             cGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAAAAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACg                                                                                                                        >  2:496624/1‑149 (MQ=255)
                                     ttAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTATCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCg                                                                                                                <  2:477957/149‑1 (MQ=255)
                                        aaTACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAGTTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGgag                                                                                                             >  1:294519/1‑149 (MQ=255)
                                        aaTACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAGTTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGgag                                                                                                             >  1:967276/1‑149 (MQ=255)
                                          tACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAGTTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGaa                                                                                                           <  2:967276/149‑1 (MQ=255)
                                          tACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAGTTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGaa                                                                                                           <  2:294519/149‑1 (MQ=255)
                                            cATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGATGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGt                                                                                                         <  1:867602/149‑1 (MQ=255)
                                                     ggTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTc                                                                                                <  1:321304/149‑1 (MQ=255)
                                                       ttttGGCAGGCCAGCATTTTTGGTGGCTTGCTTAGCAGGACCTTTCGGAAACAGTCGGTATAAATAGCTGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGATGTATTGAATTcaa                                                                                              >  1:862828/1‑147 (MQ=255)
                                                       ttttGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCa                                                                                              <  1:1004452/149‑1 (MQ=255)
                                                         ttGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGa                                                                                            <  2:725568/149‑1 (MQ=255)
                                                           ggCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATa                                                                                          <  2:595733/149‑1 (MQ=255)
                                                                      aaTTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCAGTAGAAGTATTGAATTCCAGATAGAAATCACGca                                                                               >  2:246860/1‑149 (MQ=255)
                                                                      aaTTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGACATCGCTTTTACCAGCATACGAATCGCCGGAGTAGTATTGAATTCCAGATAGAAATCCCGca                                                                               >  2:243070/1‑149 (MQ=255)
                                                                              tggCTTGCTTAGGCGGACGTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCATATAGAAATCACGCACAAAACGc                                                                       <  2:47215/149‑1 (MQ=255)
                                                                               tgcttgcttAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCAAACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGca                                                                      <  1:505535/148‑1 (MQ=255)
                                                                                gcttgcttAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGcac                                                                     >  1:644354/1‑149 (MQ=255)
                                                                                                   tCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGt                                                  <  2:800847/149‑1 (MQ=255)
                                                                                                     ggAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCGCCACTTCCCAGTGTTCCGGTGa                                                >  1:290968/1‑149 (MQ=255)
                                                                                                           aGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGa                                          <  1:188833/149‑1 (MQ=255)
                                                                                                            gTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGaa                                         <  2:26412/149‑1 (MQ=255)
                                                                                                             tCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGaaa                                        <  1:511876/149‑1 (MQ=255)
                                                                                                                gTATAAATAGCGGCTGTTAGCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACTGCGAAATcc                                     >  1:63111/1‑149 (MQ=255)
                                                                                                                gTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATAACGCACAAAACGCACCACTTACCAGAGATCCGGTGAAAGCGAAATcc                                     >  2:190844/1‑149 (MQ=255)
                                                                                                                  ataAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCt                                   >  1:332119/1‑149 (MQ=255)
                                                                                                                      aTAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGt                               >  1:880226/1‑149 (MQ=255)
                                                                                                                           ggCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTctct                          <  1:938052/149‑1 (MQ=255)
                                                                                                                           ggCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTctct                          <  2:725424/149‑1 (MQ=255)
                                                                                                                           ggCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTAACCAGCATACGAATCGCCGTAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTctct                          <  1:496624/149‑1 (MQ=255)
                                                                                                                            gCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTg                         <  2:644354/149‑1 (MQ=255)
                                                                                                                            gCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTg                         >  2:740496/1‑149 (MQ=255)
                                                                                                                               gTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCaa                      <  1:229757/149‑1 (MQ=255)
                                                                                                                               gTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCaa                      >  2:889913/1‑149 (MQ=255)
                                                                                                                               gTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCaa                      <  2:319863/149‑1 (MQ=255)
                                                                                                                                               ccAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGAAAGAAATCACGCACAAAACGCACCCCTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGATCGCTGCAATCACCACAGCCCGCgg      >  2:31453/1‑149 (MQ=255)
                                                                                                                                                cAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCAATCACCACAGCCAGCGGc     >  1:857026/1‑149 (MQ=255)
                                                                                                                                                cAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAAAGCAACACTTCCCAGTGTTCCGGTGACAGCGAAATCACTTCGTTCTCTGCAATCACCACAGCCAGCGGc     >  1:228529/1‑149 (MQ=255)
                                                                                                                                                   aTTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAATTCCCTTCGTTCTCTGCAATAACAACAGCCAGCGGCTCa  >  2:941964/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCATCCCTCGGGGGAATTAACCGATTTTCGTATCCGTTAAATACATTTTACCGGTTTTGGCAGGCCAGCAATTTTGGTGGCTTGCTTAGCCGGACCTTTCGGAAACAGTCGGTATAAATAGCGGCTGTTACCTTTTTCTTCGCCAAATTTATTCGCCATCGCTTTTACCAGCATACGAATCGCCGGAGAAGTATTGAATTCCAGATAGAAATCACGCACAAAACGCACCACTTCCCAGTGTTCCGGTGACAGCGAAATCCCTTCGTTCTCTGCAATCACCACAGCCAGCGGCTCA  >  NC_000913/1030302‑1030597

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