Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,321 C→T intergenic (‑271/‑274) trxB ← / → lrp thioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,3210CT95.0% 68.5 / ‑4.0 20intergenic (‑271/‑274)trxB/lrpthioredoxin reductase, FAD/NAD(P)‑binding/leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (6/13);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAC  >  NC_000913/932187‑932468
                                                                                                                                      |                                                                                                                                                   
tCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaa                                                                                                                                     >  1:161954/1‑151 (MQ=255)
   ccTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGaaaaaaaa                                                                                                                                  >  1:18001/1‑151 (MQ=255)
            cTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCt                                                                                                                         <  1:235488/151‑1 (MQ=255)
                  tttCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATg                                                                                                                   <  1:226476/151‑1 (MQ=255)
                          tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                           <  1:153711/151‑1 (MQ=255)
                          tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                           <  1:201192/151‑1 (MQ=255)
                            tataAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                         <  2:3761/151‑1 (MQ=255)
                               aaGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACcatgc                                                                                                      <  2:162310/151‑1 (MQ=255)
                                         aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                            <  1:219840/151‑1 (MQ=255)
                                         aaTGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                            <  2:138685/151‑1 (MQ=255)
                                         aaTGATGAATAAACGCCCCTGTTAATGAACATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATAc                                                                                            <  2:17182/151‑1 (MQ=255)
                                                  tAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                    >  1:250704/1‑150 (MQ=255)
                                                             gTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                        >  2:64091/1‑151 (MQ=255)
                                                             gTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGa                                                                        >  1:48364/1‑151 (MQ=255)
                                                                                 ttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGAc                                                    <  1:164928/151‑1 (MQ=255)
                                                                                                aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTg                                      <  1:147737/150‑1 (MQ=255)
                                                                                                  gTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATg                                   >  1:242523/1‑151 (MQ=255)
                                                                                                         tttGACAATCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAg                            <  2:83171/151‑1 (MQ=255)
                                                                                                                aTCCCCTGGTGTTTTGCGAAAATATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCca                      >  2:220686/1‑150 (MQ=255)
                                                                                                                                   aaaTATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAc  <  1:177454/151‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                   
TCACCTCTGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTCAC  >  NC_000913/932187‑932468

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