Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,770 C→T intergenic (+94/‑146) rho → / → rfe transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,7700CT96.5% 86.2 / ‑4.0 30intergenic (+94/‑146)rho/rfetranscription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/17);  total (11/19)
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.61e-01

CACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTT  >  NC_000913/3967529‑3968013
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCGCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTg                                                                                                                                                                                                                                             <  2:407545/251‑1 (MQ=255)
cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAAttt                                                                                                                                                                                                                                              <  2:305544/250‑1 (MQ=255)
   gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTg                                                                                                                                                                                                                                          >  1:34807/1‑251 (MQ=255)
    aCTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTgg                                                                                                                                                                                                                                         <  2:439103/251‑1 (MQ=255)
         ggAATAACTGCTGAAAATGTGGATCCTGCTCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTTATTAATAAACTTGCAATGACCAAGACCAATGACGATGTCTTCGAAATGATGAAACGCTCATAAATTTTTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTTTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGt                                                                                                                                                                                                                                    <  1:260643/251‑1 (MQ=255)
         ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGt                                                                                                                                                                                                                                    <  1:415048/251‑1 (MQ=255)
                  gCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCaa                                                                                                                                                                                                                             <  2:285816/249‑1 (MQ=255)
                                        aaaaCCATTCACCCTATGGGCAAAATCAATTGATGGGTATTCCTCATTAATAAACTGGCAATGACTAAGACCAATGTCGATGTCTTCGAAATGATGAAAAGCGCATAAATTTTACTTATGCCAAAAACTCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAAACTTAATGCCGCGTTGGCGAGGTTAGGAAAATTCGTGGAATTTGCTGGCATGTCAGACCATTTGCATTTCAAATGGTTAAtttt                                                                                                                                                                                                     <  1:351239/251‑1 (MQ=255)
                                                                                 ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGa                                                                                                                                                            <  1:327331/251‑1 (MQ=255)
                                                                                        aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTAGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAgg                                                                                                                                                     <  1:128725/251‑1 (MQ=255)
                                                                                              cTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCg                                                                                                                                               >  1:57003/1‑251 (MQ=255)
                                                                                                  cAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGt                                                                                                                                           <  1:219236/251‑1 (MQ=255)
                                                                                                                   gACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGttt                                                                                                                                                                                 <  2:402982/196‑1 (MQ=255)
                                                                                                                   gACGATTTCTTCGAAATGATGAAACGCCCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGttt                                                                                                                                                                                 >  1:402982/1‑196 (MQ=255)
                                                                                                                     cGATATCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTCAGTGGCGTTTTTCTTATATATATGTAATCTTAATGCCACGCAGGCGCTGTTAGGAAAATTCTTGGTATTTGATGGCATGTTATGCAATTTGACTCTCAAATGGATAATTTTTGCAAAGGACTGGTGGGTTTGGACCGGACTTTCCCTTCTGACTAAAGGACTTCGTGGTTATACTTCTTCTAATAAtt                                                                                                                        >  1:466015/1‑251 (MQ=255)
                                                                                                                              cGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTgag                                                                                                               <  1:381195/251‑1 (MQ=255)
                                                                                                                               gAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTGGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATACTTTtctatt                                                                                                                 >  1:401927/1‑245 (MQ=255)
                                                                                                                                aaaTGTTGAAACGCTCATAAAGTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGTCTTTCCCTTCTGAATAAAGGTCTGCGTGGTTATACTTCTGCTAATAATTTTCTCTGagag                                                                                                             <  2:79199/251‑1 (MQ=255)
                                                                                                                                             cTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGTCTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGa                                                                                                 >  2:337612/1‑250 (MQ=255)
                                                                                                                                                                    aaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGc                                                                                                                                                                                                                               <  2:452098/101‑1 (MQ=255)
                                                                                                                                                                    aaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGc                                                                                                                                                                                                                               >  1:452098/1‑101 (MQ=255)
                                                                                                                                                                          cACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGCATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCtcatc                                                                    >  2:471191/1‑250 (MQ=255)
                                                                                                                                                                                            tttGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTACTAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACg                                                  >  2:250675/1‑250 (MQ=255)
                                                                                                                                                                                                    tataTCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACActgtt                                         <  1:271606/250‑1 (MQ=255)
                                                                                                                                                                                                         cTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTc                                                                                                                                                               >  2:472101/1‑128 (MQ=255)
                                                                                                                                                                                                         cTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTc                                                                                                                                                               <  1:472101/128‑1 (MQ=255)
                                                                                                                                                                                                                 ttAATTCCGCTCTGGCGATGTTAGTAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAATTGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGCCTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATACTTTTCTCTGAGAGCCTGCATTGTGAATTTAATGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTTGTTTTGTTTTTTTccc                           >  1:149084/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tAGGAAAATTCTTGGAAGTTGCGGGAATGTTATGCAATTTGCATATCAAATGGTTAATTTGTGCACAGGACTGGTGGGTTTGGAACGGGCTTTCCCGTCTGAATAAAGGGCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGATCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCTTAAGGGGGCATAATAAGTc       <  2:34807/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAAAATTCTTGGAAGTTTCTGGCATGGTATGCAATTTGCATATCAAATGGTTAATTTTTGCACATGACTGGTGTGTGTGTATCTGACTTTCCCTTCTTAATAAAGGTCTTCGTGGTCATACTTCTGCGAATAATTTTCTCTGAGAGCATGCATGGTGAATTTACTGACATGGAGTACTGATCGCATCAGTAGTTTTTTATTCTCGACACTGTTGCGGTTGTTGGCCTGTTAGGGGGCAAAAAATGTCGGt    <  2:293939/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaTTCTTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGttt  <  1:461683/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
CACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCATATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTT  >  NC_000913/3967529‑3968013

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