Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,469 T→A 13.5% intergenic (+72/+19) deaD → / ← ygiD ATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,092,4690TA13.5% 85.4 / 8.4 37intergenic (+72/+19)deaD/ygiDATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base A (2/3);  total (16/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.15e-01

AGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCAT  >  NC_002947/2092321‑2092600
                                                                                                                                                    |                                                                                                                                   
agaAGCTGCTGTCCAAGCTCAAGGCCACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCtt                                                                                                                                    <  2:327526/150‑1 (MQ=255)
    gcttctGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGccc                                                                                                                                <  2:164334/150‑1 (MQ=255)
      tgctgTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCcgc                                                                                                                              <  2:429786/150‑1 (MQ=255)
        ctgTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTTCGCCAACGAGAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGa                                                                                                                            >  2:29671/1‑150 (MQ=255)
          gTCCAAGCTCAAGGACACTCACGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCAAGCCCGCGaaa                                                                                                                          <  2:694081/150‑2 (MQ=255)
             cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAgg                                                                                                                       >  2:564031/1‑150 (MQ=255)
                          acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGaa                                                                                                          >  1:374409/1‑150 (MQ=255)
                            tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGt                                                                                                              >  2:439721/1‑144 (MQ=255)
                            tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGt                                                                                                              <  1:439721/144‑1 (MQ=255)
                                  ggctggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGa                                                                                                  >  2:630959/1‑150 (MQ=255)
                                   gctggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGAt                                                                                                 <  1:630959/150‑1 (MQ=255)
                                     tggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACTGATAGATGt                                                                                               >  1:439781/1‑150 (MQ=255)
                                      ggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTc                                                                                              <  1:25576/150‑1 (MQ=255)
                                               cTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAg                                                                                     <  1:659484/150‑1 (MQ=255)
                                                         gTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAgg                                                                           <  1:341265/150‑1 (MQ=255)
                                                         gTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTCAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAgg                                                                           <  1:733652/150‑1 (MQ=255)
                                                            gCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTa                                                                        >  2:792711/1‑150 (MQ=255)
                                                             cTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTaa                                                                       >  2:479242/1‑150 (MQ=255)
                                                                                   ggTAACGAAAAAAGCCCCGGCTGAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGaa                                                 <  2:306877/149‑1 (MQ=255)
                                                                                       aCGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACtcg                                             >  1:58305/1‑148 (MQ=255)
                                                                                         gAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCACAGGGTAAACCCCTGATGCACCACGCCGAACTTGCt                                           >  1:136407/1‑150 (MQ=255)
                                                                                             aaaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCg                                       >  1:552349/1‑150 (MQ=255)
                                                                                              aaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCgg                                      <  1:564031/150‑1 (MQ=255)
                                                                                               aTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCCCGGCGAACTTGCTCCCGGc                                     <  1:29671/150‑1 (MQ=255)
                                                                                                      gCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGTCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAgg                              <  1:746427/150‑1 (MQ=255)
                                                                                                               ggggATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGTCCCGTCAGTCGAACCGATAGATGTCCATGCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAAc                     >  1:410152/1‑150 (MQ=255)
                                                                                                                gggATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACa                    <  1:479242/150‑1 (MQ=255)
                                                                                                                   aTTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAc                                                                     >  1:137552/1‑98 (MQ=255)
                                                                                                                   aTTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAc                                                                     <  2:137552/98‑1 (MQ=255)
                                                                                                                    ttttttttACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCgg                <  2:552349/150‑1 (MQ=255)
                                                                                                                    ttttttttACCTGTACCGGCCCTTTCCCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCACCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCgg                <  2:114846/150‑1 (MQ=255)
                                                                                                                                gTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtc    <  1:792711/150‑1 (MQ=255)
                                                                                                                                  aCCGGCCCTTTCGCGGGCTTGCCCGCGCATAGGCGCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGGAAACCCCTGATGCACCCCGCCGAACTTGCTGCCGGCCCCAAGGGCAAAGAACAGCGGCAGCaaaatctcgt  >  2:407542/1‑143 (MQ=255)
                                                                                                                                  aCCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtcat  <  2:439781/150‑1 (MQ=255)
                                                                                                                                        ctcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtca   >  1:786322/4‑143 (MQ=255)
                                                                                                                                        ctcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtca   <  2:786322/140‑1 (MQ=255)
                                                                                                                                         tcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtca   >  2:581496/3‑142 (MQ=255)
                                                                                                                                         tcttCGCGGGCAAGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCtca   <  1:581496/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                   
AGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCCCAAGGGCAAAGAACAGCGGCAGCAAATGCTCAT  >  NC_002947/2092321‑2092600

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