| 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 | ||
| * | CP000730 | 1,318,237 | 0 | T | A | 99.5% | 73.3 / ‑6.9 | 27 | F124Y (TTT→TAT) | miaB | 2‑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: ‑ |