Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,423 G→A 18.5% intergenic (+26/+65) 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,4230GA18.5% 65.2 / 7.6 33intergenic (+26/+65)deaD/ygiDATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base A (4/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGCGACAGCTTCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGC  >  NC_002947/2092274‑2092565
                                                                                                                                                     |                                                                                                                                              
ggTGCGCGACAGCTTCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCgg                                                                                                                                                <  2:103114/150‑1 (MQ=255)
                 gcttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAtttt                                                                                                                               >  2:502658/1‑150 (MQ=255)
                  cttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttt                                                                                                                              >  2:252795/1‑150 (MQ=255)
                    tGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttttt                                                                                                                            >  2:249730/1‑150 (MQ=255)
                    tGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAttttttc                                                                                                                            <  1:485757/150‑2 (MQ=255)
                                         gcatCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTcgc                                                                                                       <  1:502658/147‑1 (MQ=255)
                                          ccTCGAGAAGCTGCTGTCCCAGCTCAAGGAAACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGCCTTAGGTCGGGGGTTTTTTTTACCTGTACCGGCCCTTTcccg                                                                                                      >  2:154708/1‑150 (MQ=255)
                                            tCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCggg                                                                                                    >  2:558459/1‑150 (MQ=255)
                                                        tgTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGCTCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGaa                                                                                        <  1:610971/150‑1 (MQ=255)
                                                             aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAttttttc                                                                                                                            <  2:383770/109‑2 (MQ=38)
                                                             aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAttttttc                                                                                                                            >  1:383770/1‑108 (MQ=38)
                                                              aGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                  <  1:252795/150‑1 (MQ=255)
                                                                    aGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATAAACCACATTCGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCGTGCCCGCGAAGAGGCCCTTCAg                                                                            >  1:13404/1‑150 (MQ=255)
                                                                         acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc                                                                                                                         <  2:598960/100‑1 (MQ=255)
                                                                         acTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc                                                                                                                         >  1:598960/1‑100 (MQ=255)
                                                                           tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCa                                                                             >  1:220465/1‑142 (MQ=255)
                                                                           tCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCa                                                                             <  2:220465/142‑1 (MQ=255)
                                                                                        aaGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCtt                                                                               <  2:96018/127‑1 (MQ=255)
                                                                                        aaGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCtt                                                                               >  1:96018/1‑127 (MQ=255)
                                                                                                  aGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAgcgcg                                              <  1:249730/150‑1 (MQ=255)
                                                                                                     tGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa                                                   <  1:660894/142‑1 (MQ=255)
                                                                                                     tGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa                                                   >  2:660894/1‑142 (MQ=255)
                                                                                                             tACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATc                                                       <  1:440916/130‑1 (MQ=255)
                                                                                                             tACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATc                                                       >  2:440916/1‑130 (MQ=255)
                                                                                                                 gcgAGGATTGCTCCTGGCGTAACGAAAAAACCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCTTGACCGCGAAGAGGCCCGTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAcccc                               <  1:154708/150‑1 (MQ=255)
                                                                                                                               tGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTTAGTCGAACCGATAGATGTCCAGCGCGAGTGCGCCCAGGGTAAAACCCTGATGCACCACGcc                 >  2:552360/1‑150 (MQ=255)
                                                                                                                                  cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc                                <  2:306330/132‑1 (MQ=255)
                                                                                                                                  cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAAccc                                >  1:306330/1‑132 (MQ=255)
                                                                                                                                          aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa                                                   <  2:108466/105‑1 (MQ=38)
                                                                                                                                          aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCa                                                   >  1:108466/1‑105 (MQ=38)
                                                                                                                                           aaaaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTccc     >  1:562855/1‑150 (MQ=255)
                                                                                                                                              aTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTACCTGCAAAGTTGCCCTTAAGAGGATGCTGGAGATTGCGATCCCCAGCGCGCGCAGGGTAAACCCCTGCTGCACCACGCCGAACTTGCTCCCgcc  >  1:279522/1‑148 (MQ=255)
                                                                                                                                                 cccGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg                                         >  1:86119/1‑108 (MQ=255)
                                                                                                                                                 cccGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAg                                         <  2:86119/108‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGTGCGCGACAGCTTCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGC  >  NC_002947/2092274‑2092565

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