Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,229 C→T *233* (TAG→TAA thiQ ← thiamine/thiamine pyrophosphate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,2290CT100.0% 125.0 / NA 40*233* (TAG→TAAthiQthiamine/thiamine pyrophosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/16);  total (24/16)

CGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72095‑72361
                                                                                                                                      |                                                                                                                                    
cGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc                                                                                                                                  <  1:537071/139‑1 (MQ=255)
   aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAGTGCATATTAACCCg                                                                                                                               <  2:605729/139‑1 (MQ=255)
          ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATccc                                                                                                                        >  2:581774/1‑139 (MQ=255)
          ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATccc                                                                                                                        >  2:1047476/1‑139 (MQ=255)
                   aagagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:424074/3‑123 (MQ=255)
                   aagagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:424074/121‑1 (MQ=255)
                     gagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             <  1:1047476/139‑1 (MQ=255)
                     gagaTAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCg                                                                                                             <  1:802294/139‑1 (MQ=255)
                         tAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             >  1:529939/1‑119 (MQ=255)
                         tAGTCTTGTGCGGGTTGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa                                                                                                                             <  2:529939/119‑1 (MQ=255)
                          aGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        >  1:575757/1‑139 (MQ=255)
                          aGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCa                                                                                                        <  1:885020/139‑1 (MQ=255)
                            tCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACt                                                                                                      <  1:581774/139‑1 (MQ=255)
                                gtgCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGcc                                                                                                  >  2:271476/1‑139 (MQ=255)
                                          ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTca                                                                                        >  2:923699/1‑139 (MQ=255)
                                                    cGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  1:633359/1‑90 (MQ=255)
                                                    cGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  2:633359/90‑1 (MQ=255)
                                                                         ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGAt                                                         >  1:351429/1‑139 (MQ=255)
                                                                         ccTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGAt                                                         >  2:806303/1‑139 (MQ=255)
                                                                            aCAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATgcg                                                      >  2:342865/1‑139 (MQ=255)
                                                                             cAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATgcgc                                                     >  2:290887/1‑139 (MQ=255)
                                                                                  cGCACTACCCGTAGGTCGGCTAAGGCGTTCACGCAGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCTGGCGATGCGCCCGTc                                                >  1:706324/1‑139 (MQ=255)
                                                                                  cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               <  1:36864/60‑1 (MQ=255)
                                                                                  cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg                                                                                                                               >  2:36864/1‑60 (MQ=255)
                                                                                  cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTc                                                >  1:874810/1‑139 (MQ=255)
                                                                                  cGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTc                                                >  2:925925/1‑139 (MQ=255)
                                                                                   gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                               <  1:923699/139‑1 (MQ=255)
                                                                                   gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                               >  1:865155/1‑139 (MQ=255)
                                                                                   gCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                               <  1:342865/139‑1 (MQ=255)
                                                                                       tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                               <  1:595110/103‑1 (MQ=255)
                                                                                       tACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCAt                                                                               >  2:595110/1‑103 (MQ=255)
                                                                                        aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                               <  1:577798/134‑1 (MQ=255)
                                                                                        aCCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCg                                               >  2:577798/1‑134 (MQ=255)
                                                                                                   ggATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAgc                               >  1:292493/1‑139 (MQ=255)
                                                                                                      tAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGcgcg                            >  2:367561/1‑139 (MQ=255)
                                                                                                               cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTGATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATcc                   >  2:755469/1‑139 (MQ=255)
                                                                                                               cACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGTGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATcc                   >  2:202245/1‑139 (MQ=255)
                                                                                                                      cATCCGACAGTGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGcgccgc            <  1:925925/139‑1 (MQ=255)
                                                                                                                                tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg  <  2:874810/139‑1 (MQ=255)
                                                                                                                                tGCATATTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACg  >  1:526376/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACG  >  NC_000913/72095‑72361

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