Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,092,423 G→A 12.8% 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,4230GA12.8% 72.7 / 2.3 31intergenic (+26/+65)deaD/ygiDATP‑dependent DEAD‑box RNA helicase DeaD/4,5‑DOPA‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base A (2/2);  total (15/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.90e-01

TCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCC  >  NC_002947/2092288‑2092567
                                                                                                                                        |                                                                                                                                                
tcagcttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCGAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAt                                                                                                                                    >  1:425853/1‑150 (MQ=255)
     ttGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGAtttttt                                                                                                                               >  1:225002/1‑150 (MQ=255)
         tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACAGTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCgg                                                                                                                                        <  1:519646/137‑1 (MQ=255)
         tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTAc                                                                                                                           <  1:504468/150‑1 (MQ=255)
         tCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAAACCCCGGCTTAGGTCgg                                                                                                                                        >  2:519646/1‑137 (MQ=255)
                  cGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAtt                                                                                                                                   <  2:84344/133‑1 (MQ=255)
                  cGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGACCTAAGCCGGGGAtt                                                                                                                                   >  1:84344/1‑133 (MQ=255)
                       gACGGCCTCGAGAAGCTGCAGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGGCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCtt                                                                                                             >  1:19150/1‑150 (MQ=255)
                           gCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGTTCCTGGCGTAACGAAAAAACCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTcgc                                                                                                         <  2:249264/150‑1 (MQ=255)
                                   aaGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGc                                                                                                 <  2:225002/150‑1 (MQ=255)
                                              cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaaggg                                                                                      >  1:396213/1‑146 (MQ=255)
                                              cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAaga                                                                                        >  1:39105/1‑148 (MQ=255)
                                              cAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAaga                                                                                        <  2:39105/148‑1 (MQ=255)
                                               aaGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagggc                                                                                     <  2:396213/150‑6 (MQ=255)
                                                     aaGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCgg                                                                                                                  >  2:576080/1‑115 (MQ=255)
                                                     aaGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCgg                                                                                                                  <  1:576080/115‑1 (MQ=255)
                                                       ggACTCTCGCGTGGCTGGCAAGCAGCTGACGATGCAGGCGTACCGCGAGGATTGATCCTTGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGt                                                                             <  2:299634/150‑1 (MQ=255)
                                                                   ggctggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCTTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGa                                                                 <  1:187707/150‑1 (MQ=255)
                                                                      tggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagg                                                                                       <  1:562710/125‑4 (MQ=38)
                                                                      tggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCTCTTCGCGGGCAAGCCCGCGaaagg                                                                                       >  2:562710/1‑122 (MQ=38)
                                                                       ggcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTc                                                             <  1:521678/150‑1 (MQ=255)
                                                                        gcAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTcc                                                            >  1:176985/1‑150 (MQ=255)
                                                                                              gCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTa                                       <  2:425853/150‑1 (MQ=255)
                                                                                                      gAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTg                               >  1:313833/1‑150 (MQ=255)
                                                                                                          aTTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGc                           >  2:579350/1‑150 (MQ=255)
                                                                                                               tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                     <  2:291540/87‑1 (MQ=255)
                                                                                                               tCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                     >  1:291540/1‑87 (MQ=255)
                                                                                                                     cGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGaa                <  1:579350/150‑1 (MQ=255)
                                                                                                                            aaaaaaTCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCACTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCCGCGCGCCCATGGTAAACCCCTTATGCACCACGCCTAACTTGCTc         >  2:76483/1‑150 (MQ=255)
                                                                                                                             aaaaaTCCCCGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTcc        >  1:209774/1‑150 (MQ=255)
                                                                                                                                   ccccGACCTAAGCCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGccc  <  2:209774/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                
TCAGCTTGGTCGAGCTGCCGGAAGACGGCCTCGAGAAGCTGCTGTCCAAGCTCAAGGACACTCGCGTGGCTGGCAAGCAGCTGAAGCTGCGTCGCTACCGCGAGGATTGATCCTGGCGTAACGAAAAAATCCCCGGCTTAGGTCGGGGATTTTTTTTACCTGTACCGGCCCTTTCGCGGGCTTGCCCGCGAAGAGGCCCTTCAGTCGAACCGATAGATGTCCATCCCCAGCGCGCCCAGGGTAAACCCCTGATGCACCACGCCGAACTTGCTCCCGGCCC  >  NC_002947/2092288‑2092567

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