Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,318,237 T→A 100% F124Y (TTT→TAT)  miaB → 2‑methylthioadenine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,318,2370TA99.5% 73.3 / ‑6.9 27F124Y (TTT→TAT) miaB2‑methylthioadenine synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (18/8);  minor base G (1/0);  total (19/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGAATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT  >  CP000730/1318102‑1318373
                                                                                                                                       |                                                                                                                                        
tgaatgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTg                                                                                                                                     >  2:228737/1‑141 (MQ=255)
   atgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGaaa                                                                                                                                  <  1:102244/141‑1 (MQ=255)
   atgCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGaaa                                                                                                                                  <  1:31088/141‑1 (MQ=255)
                    gTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCAt                                                                                                                 >  1:80523/1‑141 (MQ=255)
                                ataCTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGAGTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGaa                                                                                                     >  2:5825/1‑141 (MQ=255)
                                ataCTTAAAGCATTAGGCTATCAAGCAACGACAGAAATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGAGAAATAAGTAATTTGAAGCATTTGAGAAAAGaa                                                                                                     >  2:219656/1‑141 (MQ=255)
                                  aCTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAACTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACg                                                                                                   <  1:5229/141‑1 (MQ=255)
                                    ttGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGGGATTAGAGATAATGCCGAGTACACAGTGTATAGGGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGAc                                                                                                 >  2:5874/1‑141 (MQ=255)
                                                    tCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAAGTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCgg                                                                                 <  1:196134/141‑1 (MQ=255)
                                                       aGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGgtgt                                                                              <  1:259065/141‑1 (MQ=255)
                                                          aaCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtg                                                                           >  2:35854/1‑141 (MQ=255)
                                                           aCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtgt                                                                          <  1:161565/141‑1 (MQ=255)
                                                           aCGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTtgt                                                                          >  2:14517/1‑141 (MQ=255)
                                                                 gATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGt                                                                    >  1:19288/1‑141 (MQ=255)
                                                                     ttAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTGGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGt                                                                >  2:57467/1‑141 (MQ=255)
                                                                               aGATGTTATTTTAATCAATACATGTGCGATTAGGGAAAATGCCGAGAACAAAGGGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAagaaga                                                      <  2:175430/141‑1 (MQ=255)
                                                                                gATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGGATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAagag                                                     >  1:105458/1‑141 (MQ=255)
                                                                                          taattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTgc                                           >  1:10378/1‑140 (MQ=255)
                                                                                           aattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGaa                                          >  2:18109/1‑141 (MQ=255)
                                                                                           aattaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGaa                                          >  2:156545/1‑141 (MQ=255)
                                                                                            attaatACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGCAt                                         >  2:212879/1‑141 (MQ=255)
                                                                                                   cATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGATATATGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAAttt                                  >  1:151066/1‑141 (MQ=255)
                                                                                                     tgtgCGATTAGAGAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTa                                >  2:92360/1‑141 (MQ=255)
                                                                                                              agagAAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTaaactagaag                       >  2:140676/1‑134 (MQ=255)
                                                                                                               gagaAAATGCCGAGAACAAAGTGTATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTCGTGAATAAAATTTTAAAATCGTATc                      >  2:98980/1‑141 (MQ=255)
                                                                                                                 gaAAGTGCGCAGGACAAAGTGTATAGTGAAATAGGTAATTGGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTTTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCaa                    <  2:12970/141‑1 (MQ=255)
                                                                                                                                   gtgtATAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATAttt  >  1:78712/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TGAATGCACATGACACTGAGGTCATTGCTGGTATACTTGAAGCATTAGGCTATCAAGCAACGACTGATATTAACACTGCAGATGTTATTTTAATTAATACATGTGCGATTAGAGAAAATGCCGAGAACAAAGTGTTTAGTGAAATAGGTAATTTGAAGCATTTGAAAAAAGAACGACCTGATATTTTAATCGGTGTTTGTGGTTGTATGTCACAAGAAGAGTCAGTAGTGAATAAAATTTTAAAATCGTATCAAAATGTAGATATGATATTT  >  CP000730/1318102‑1318373

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