Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,572,316 T→C 9.3% *472* (TAA→TGA)  gltX ← glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,572,3160TC 9.3% 107.2 / 4.2 43*472* (TAA→TGA) gltXglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (2/2);  total (22/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 = 1.00e+00

TTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  minE/1572247‑1572382
                                                                     |                                                                  
ttGACTAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAtt                                                                    <  1:559727/70‑1 (MQ=255)
  gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATCAct                                                                 >  1:1061075/1‑71 (MQ=255)
     gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                              >  1:1307796/1‑71 (MQ=255)
                aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATCACTGCTGAtttt                                                        <  1:1030235/66‑1 (MQ=255)
                 ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                   >  1:1184102/1‑70 (MQ=255)
                  taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                 <  1:460578/71‑1 (MQ=255)
                      aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATCACTGCTGATTTTCGCGTTCAGCa                                             <  1:159846/71‑1 (MQ=255)
                           gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        <  1:954328/71‑1 (MQ=255)
                           gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        <  1:443852/71‑1 (MQ=255)
                           gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        <  1:1932366/71‑1 (MQ=255)
                                atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                   >  1:1039830/1‑39 (MQ=255)
                                  attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                 >  1:1080874/1‑71 (MQ=255)
                                   ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                                <  1:613531/71‑1 (MQ=255)
                                     atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                  >  1:777740/1‑67 (MQ=255)
                                     atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                              >  1:1270234/1‑71 (MQ=255)
                                          cTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                     >  1:357707/1‑59 (MQ=255)
                                          cTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  1:1791316/71‑1 (MQ=255)
                                             ccGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       <  1:541194/70‑1 (MQ=255)
                                             ccGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       <  1:1815318/70‑1 (MQ=255)
                                             ccGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       <  1:2043706/70‑1 (MQ=255)
                                               gTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        <  1:1921482/51‑1 (MQ=255)
                                                tttATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                         <  1:322736/49‑1 (MQ=255)
                                                tttATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                         <  1:104034/49‑1 (MQ=255)
                                                tttATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                     <  1:868032/53‑1 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                       >  1:1510049/1‑49 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                         <  1:181223/47‑1 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                         <  1:906691/47‑1 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                     >  1:778116/1‑51 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                     >  1:35185/1‑51 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                     >  1:175428/1‑51 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                   >  1:1097675/1‑53 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  1:508992/1‑65 (MQ=255)
                                                  tATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                 >  1:2202546/1‑71 (MQ=255)
                                                    tCTTTTTACACGCTAATCACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  1:2108457/1‑63 (MQ=255)
                                                      tttttACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                             >  1:267363/1‑39 (MQ=255)
                                                      tttttACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                            <  1:406746/56‑1 (MQ=255)
                                                        tttACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             <  1:668799/69‑1 (MQ=255)
                                                        tttACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             <  1:1135232/69‑1 (MQ=255)
                                                        tttACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             <  1:981247/69‑1 (MQ=255)
                                                        tttACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc           >  1:1170090/1‑71 (MQ=255)
                                                           acacGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         >  1:773204/1‑54 (MQ=255)
                                                                  aaTTACTGCTGATTTTCGCGTTCAGCAATAAAATCGAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  >  1:29731/1‑70 (MQ=255)
                                                                  aaTTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                     >  1:275202/1‑51 (MQ=255)
                                                                  aaTTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  >  1:1419984/1‑70 (MQ=255)
                                                                     |                                                                  
TTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  minE/1572247‑1572382

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