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.C97.2% 123.8 / ‑3.1 36intergenic (‑105/+45)rluA/minE23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (20/15);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
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.

CAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCA  >  AE015451/1932408‑1932858
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                             
cAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑ccc                                                                                                                                                                                                                            <  1:764463/234‑2 (MQ=255)
                       tgttgcCCTCGGCACGGCAAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAg                                                                                                                                                                                                             <  1:113595/226‑1 (MQ=255)
                                cGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCa                                                                                                                                                                                                          >  2:527810/1‑220 (MQ=255)
                                    aCGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGCTATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAggg                                                                                                                                                                                                                    >  1:167669/1‑206 (MQ=255)
                                           ggTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑cccgcga                                                                                                                                                                                                                        >  2:318622/1‑190 (MQ=255)
                                            gTACCGACAGCGGCGGGGTCGGTTTGTTGGTCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACAGTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑ccc                                                                                                                                                                                                                            <  2:168487/190‑2 (MQ=255)
                                                         cagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcgcg                                                                                                                                                                                    >  2:54487/1‑217 (MQ=255)
                                                                      ttgATCACCAGGATGGCGGCGTCGTCGACGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                                                                                                              <  1:321873/194‑1 (MQ=255)
                                                                       tgATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAGGGGACGCATATCCTTTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAggg                                                                                                                                                                 <  2:345380/222‑1 (MQ=37)
                                                                                  aTGGGGGCGTCTGCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                                                                                                            <  1:157494/216‑1 (MQ=255)
                                                                                  aTGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGc                                                                                                                                                                                                        <  2:66894/172‑1 (MQ=255)
                                                                                       ggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                                                                                                        >  1:668732/1‑215 (MQ=255)
                                                                                         cgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTACGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCa                                                                                                                                                                                                    >  2:164036/1‑169 (MQ=255)
                                                                                                                 tCGACAGCGGCATTGCAGGTCTCACGAAAAGAGCAGTAAAGGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGTCAGGGTGAt                                                                                                                                                             <  2:395270/184‑1 (MQ=255)
                                                                                                                    aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGAc                                                                                                                                                                                          >  2:363992/1‑152 (MQ=255)
                                                                                                                                 aGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTTCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGa                                                                                                                                                                                           >  2:655491/1‑138 (MQ=255)
                                                                                                                                                     tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCAc                                                                      <  1:54487/235‑1 (MQ=255)
                                                                                                                                                       aaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCAc                                                                      >  2:639045/1‑233 (MQ=255)
                                                                                                                                                                gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACata                                                                                  >  2:439008/1‑212 (MQ=255)
                                                                                                                                                                              cTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGc                                                                             <  1:61305/203‑1 (MQ=255)
                                                                                                                                                                                 cgcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCca                                                                 >  1:154175/1‑212 (MQ=255)
                                                                                                                                                                                 cgcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAgc                                                >  2:637941/1‑229 (MQ=255)
                                                                                                                                                                                    gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGCGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGa                                                                          >  2:100019/1‑200 (MQ=255)
                                                                                                                                                                                    gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCCCATATTTg                                                                              >  1:65926/1‑196 (MQ=255)
                                                                                                                                                                                        gaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGaa                                                                                                                                               <  2:104239/126‑1 (MQ=255)
                                                                                                                                                                                               cTCCTACAAAGGATATGCGTCCCCTGTAGGCGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt                                                                                           <  1:253561/172‑1 (MQ=255)
                                                                                                                                                                                                             atGCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTTCg                                                                            >  2:571160/1‑173 (MQ=255)
                                                                                                                                                                                                               gCGTCCCCTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCt           >  1:315958/1‑236 (MQ=255)
                                                                                                                                                                                                                  tCCCCTGTAGGAGCGGGCGTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCtgc         <  1:194234/235‑3 (MQ=255)
                                                                                                                                                                                                                    cccTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGt        <  2:710660/234‑1 (MQ=255)
                                                                                                                                                                                                                     ccTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGt                              >  1:755746/1‑211 (MQ=255)
                                                                                                                                                                                                                     ccTGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGt                              <  2:755745/211‑1 (MQ=255)
                                                                                                                                                                                                                       tGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcg                                                                                                                                                                                      >  2:143241/1‑57 (MQ=18)
                                                                                                                                                                                                                       tGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCg                      >  2:244430/1‑217 (MQ=255)
                                                                                                                                                                                                                       tGTAGGAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTgcgc    >  1:594039/1‑235 (MQ=255)
                                                                                                                                                                                                                           ggAGCGGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCa  <  2:315957/233‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCcagc                                                               <  2:646298/165‑1 (MQ=255)
                                                                                                                                                                                                                                  gCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCcagc                                                               >  1:646299/1‑165 (MQ=255)
                                                                                                                                                                                                                                   cTTG‑CCAGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCACTGGTTTTCCAGCTCGat                                                                                                                <  2:76742/115‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                             
CAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCA  >  AE015451/1932408‑1932858

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: ‑