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.C94.6% 126.5 / 1.0 37intergenic (‑105/+45)rluA/minE23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (20/15);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑G‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCT  >  AE015451/1932424‑1932867
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
cAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCa                                                                                                                                                                                                                   >  1:636533/1‑236 (MQ=255)
     cTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGc                                                                                                                                                                                                                 <  2:793925/233‑1 (MQ=255)
                 ggCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                                                                                                                  <  1:195782/236‑1 (MQ=255)
                                     gcagcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                                                                                                                >  1:456286/1‑234 (MQ=255)
                                                           cACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAg                                                                                                                                                          >  2:791757/1‑234 (MQ=255)
                                                               aGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                                                                                                                    >  1:505298/1‑236 (MQ=255)
                                                                               cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGaa                                                                                                                                                                                                                                >  2:772117/1‑144 (MQ=255)
                                                                                      gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTCC‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                                                                                                      >  1:862918/1‑163 (MQ=255)
                                                                                                    aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGGGATGTTCAGCGCCAGAAACGAGCAGCTGCCCTGGtt                                                                                                                                     >  1:281962/1‑214 (MQ=255)
                                                                                                     cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCg                                                                                                                                                                                     >  2:780263/1‑165 (MQ=255)
                                                                                                     cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                                                                                                                 >  2:610663/1‑201 (MQ=255)
                                                                                                       gCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTg                                                                                                                                             <  2:456283/203‑1 (MQ=255)
                                                                                                          gCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                                                                                                         <  1:295512/236‑1 (MQ=255)
                                                                                                              tGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATTCACATCATCGTTGCCgc                                                                                                     >  2:784951/1‑235 (MQ=255)
                                                                                                                      tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGc                                                                                                                                                                                                                 <  2:197177/120‑1 (MQ=255)
                                                                                                                           aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGa                                                                                                                                                                       >  2:871421/1‑157 (MQ=255)
                                                                                                                                     tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGGATCGagcagc                                                                                                                                               >  1:826424/1‑171 (MQ=255)
                                                                                                                                         cGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTCCCAGCTCGATATGCACATCATCGTTg                                                                                                         <  2:894020/205‑1 (MQ=255)
                                                                                                                                                       tGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                                                                                                           >  2:150468/1‑125 (MQ=255)
                                                                                                                                                       tGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTg                                                                                                                                             <  1:6147/155‑1 (MQ=255)
                                                                                                                                                                tCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGc                                                   <  1:610666/236‑1 (MQ=255)
                                                                                                                                                                 cgcgGGCAAGCCAGCTCCTACAACGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCAAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                                                                                                                     <  1:232907/121‑1 (MQ=38)
                                                                                                                                                                                  cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat                                                                                                                 >  2:886707/1‑156 (MQ=255)
                                                                                                                                                                                  cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat                                                                                                                 <  1:886710/156‑1 (MQ=255)
                                                                                                                                                                                   tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAgcg                                                        >  1:703895/1‑212 (MQ=255)
                                                                                                                                                                                        aGGATATGCGTCCCCTGTAGGCGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCt                                  <  1:117507/229‑1 (MQ=255)
                                                                                                                                                                                            tatGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGc                          <  2:636530/233‑1 (MQ=255)
                                                                                                                                                                                              tGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGt                                                                     <  2:703892/188‑1 (MQ=255)
                                                                                                                                                                                                 gTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAg                                                                                                                                                                                                                      >  1:660349/1‑40 (MQ=14)
                                                                                                                                                                                                 gTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTc                   <  1:321701/235‑1 (MQ=255)
                                                                                                                                                                                                  tCCCCTGTCGGGGCGGGCTT‑G‑CCGGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCca                                                                          <  2:120869/179‑1 (MQ=255)
                                                                                                                                                                                                        gTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCCGCTCCAGAATCGAGCAGCTGCCCTGGCTTTCCAGCTCGATATGCCCATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCACCAGTTCCTTCTGc                                                            >  2:36014/1‑187 (MQ=255)
                                                                                                                                                                                                                  gCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                                                                                                                            <  1:801994/113‑1 (MQ=255)
                                                                                                                                                                                                                  gCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                                                                                                                            >  2:801991/1‑113 (MQ=255)
                                                                                                                                                                                                                  gCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCg                                                                                                      >  2:58423/1‑135 (MQ=255)
                                                                                                                                                                                                                  gCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCt                                  >  2:378648/1‑203 (MQ=255)
                                                                                                                                                                                                                  gCTT‑GCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGCTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCt  >  2:192346/1‑235 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑G‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCT  >  AE015451/1932424‑1932867

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