Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,932,638:1 +C 100% intergenic (‑105/+45) rluA ← / ← minE 23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,932,6381.C100.0% 90.5 / NA 26intergenic (‑105/+45)rluA/minE23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/14);  total (12/14)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATG  >  AE015451/1932410‑1932863
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
gACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcg                                                                                                                                                                                                                             <  1:259471/235‑5 (MQ=255)
                                  aCGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGCTATGCGTCCCCTGTAGGCGCGGGCTTGcc                                                                                                                                                                                                                                 >  1:10571/1‑197 (MQ=255)
                                               gACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                                                                                                              <  1:510195/220‑1 (MQ=255)
                                                         gggTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                                                                                                       <  2:124431/233‑1 (MQ=255)
                                                                                aTGGCGGCGTCGTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                                                                                                                   <  2:464831/182‑1 (MQ=255)
                                                                                                       cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                                                                                                            >  2:459463/1‑182 (MQ=255)
                                                                                                          gACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                                                                                                                 <  1:328204/190‑1 (MQ=255)
                                                                                                           gcctgcGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCACGCTCCTACAAAGGATATGCGTCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                                                                                                              <  1:97334/155‑1 (MQ=39)
                                                                                                                  aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                                                                                                               >  1:352207/1‑200 (MQ=255)
                                                                                                                   cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat                                                                                                          >  1:225854/1‑236 (MQ=255)
                                                                                                                       ggCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACGAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGa                                                                                                                                                     >  2:250804/1‑189 (MQ=255)
                                                                                                                       ggCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt                                                                                                      >  1:297843/1‑236 (MQ=255)
                                                                                                                                   ctcACGAAGATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                                                                                                                            <  1:406903/202‑1 (MQ=255)
                                                                                                                                   ctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                                                                                                    >  2:152659/1‑226 (MQ=255)
                                                                                                                                       ggAAAATAGCAGTAAACGCTGGGGCCGATTTGCGGCCCTTGGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCa                                                                                                                                                                                                                    <  1:11191/109‑1 (MQ=38)
                                                                                                                                         aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAtt                                                                                     >  1:119662/1‑235 (MQ=255)
                                                                                                                                                      aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                                                                                                            <  1:379071/135‑1 (MQ=255)
                                                                                                                                                         cTGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcg                                                                                                                                                                                                                             <  2:74533/82‑5 (MQ=255)
                                                                                                                                                              gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCgg                                                                                >  2:387970/1‑219 (MQ=255)
                                                                                                                                                                        ggCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGTGCTGCAAGGCAGCCCCAAATTGTCGTCTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCg                                                                                                                                                 <  2:233282/144‑1 (MQ=255)
                                                                                                                                                                                          ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGt                                   <  1:351557/236‑1 (MQ=255)
                                                                                                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                                                                                                              <  1:279906/90‑1 (MQ=255)
                                                                                                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                                                                                                              >  2:279902/1‑90 (MQ=255)
                                                                                                                                                                                                   cagaGGGTATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                                                                                                     <  2:461392/158‑1 (MQ=255)
                                                                                                                                                                                                     aaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCcgcg                                                    >  2:439338/1‑208 (MQ=255)
                                                                                                                                                                                                             gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCt                >  1:103679/1‑236 (MQ=255)
                                                                                                                                                                                                                 ccccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGt             <  1:250808/235‑1 (MQ=255)
                                                                                                                                                                                                                        aGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGCTGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGc                               >  2:149783/1‑210 (MQ=255)
                                                                                                                                                                                                                               ggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACgatg  <  1:152663/232‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
GACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATG  >  AE015451/1932410‑1932863

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