Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,042 G→T 8.7% C232F (TGC→TTC)  queF → 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,0420GT 8.7% 111.6 / 2.9 46C232F (TGC→TTC) queF7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base T (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

TGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGAC  >  NC_000913/2925901‑2926182
                                                                                                                                              |                                                                                                                                            
tgctgAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGc                                                                                                                                             >  2:297091/1‑143 (MQ=255)
  ctgAAATCAAACTGCCTGAGCCCCCATCAACCAGATTGGGGTTCGGTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGt                                                                                                                                           <  2:882820/143‑1 (MQ=255)
      aaTCAAACTGCCTGATCACCCATCAACCAGATTGGTGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGaa                                                                                                                                       <  2:604988/143‑1 (MQ=255)
      aaTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGTACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCCTCGAAACGAGCTCC‑CGAACAGTGCGTGGAAc                                                                                                                                      >  1:402748/1‑143 (MQ=255)
      aaTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGaa                                                                                                                                       <  2:521919/143‑1 (MQ=255)
          aaaCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACATAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCa                                                                                                                                   >  2:451300/1‑143 (MQ=255)
            tcTGCCTGATAACCGATCAAGCAGATTGGAGTTCGAACCAAAGTAAGAATCGTCAACGGCAAATTGACAGACAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGAActcgtt                                                                                                                                 <  1:34006/142‑6 (MQ=255)
              tGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCtt                                                                                                                               >  1:505776/1‑143 (MQ=255)
                 cTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTaa                                                                                                                            <  1:826358/143‑1 (MQ=255)
                     tCACACATCAACCAGATGGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGAc                                                                                                                        <  2:147844/143‑1 (MQ=255)
                       aCCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCCTCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCt                                                                                                                      >  2:13034/1‑143 (MQ=255)
                       aCCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCAGACATCAAAACGAGTTCCACGAACAGTGCGTGCAACGCATCTTTAATGcccc                                                                                                                      >  2:344174/1‑142 (MQ=255)
                         ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑GACCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCAACGAACAGCGCGTGGAACGCATAGTTAATGACCTGt                                                                                                                    >  2:499968/1‑143 (MQ=255)
                         ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                                    <  1:634656/143‑1 (MQ=255)
                         ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                                    >  1:407238/1‑143 (MQ=255)
                           aTCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGGTa                                                                                                                  >  2:7619/1‑143 (MQ=255)
                                ccAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCtt                                                                                                             >  1:792557/1‑143 (MQ=255)
                                ccAGATTGCGGTTCCCTCCACACTCAGTATCGTGGACGCCAAATTGA‑GGAGAACGGTCGGCTGTGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCtt                                                                                                             <  1:663832/143‑1 (MQ=255)
                                 cAGATTGGGGTTCGCTCAAAATCCAGTATCGTGGACGCCAACTTGACAGAGAA‑ACACTGCTGCGTTGCTTTGTCTCATTTCGTGCTCTCAACCAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                            <  2:403982/143‑1 (MQ=255)
                                  aGATTGTGGTTCGTTCTAAATTCAGTATCGTGGTTGCCAAATTGCCGGAGAA‑AAACTGCTGCGTTACCTGGTCTCATGCCGTCATCACTACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCt                                                                                                           <  2:334517/143‑1 (MQ=255)
                                  aGATTGGGGTTCGCTACAAATTCAGTATCGTCGACGCCAAATTGACAGAGAA‑AAACTGCCGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCt                                                                                                           <  2:447988/143‑1 (MQ=255)
                                  aGATGGGGGTTCGCTCCACAGTCAGTAGCGTGGACGCCAAATGGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTGCCACGAACAGTGCGTGGTACGCATCTTGAATGACCTGTTACGCTTCt                                                                                                           <  1:111445/143‑1 (MQ=255)
                                        gggTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCCTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgccagc                                                                                                     <  1:126382/143‑1 (MQ=255)
                                                 ccAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAAt                                                                                            >  2:701957/1‑143 (MQ=255)
                                                         aGTACCGTGGACGCCAAATTGACAGAGAA‑AAACTGCCGCGTGACCTGGGCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGtt                                                                                    <  1:567237/143‑1 (MQ=255)
                                                                tggctgCCAATTTGACAGAGTA‑ACTCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCAc                                                                             <  2:53320/138‑1 (MQ=255)
                                                                       cAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATAACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATAc                                                                      <  1:79136/143‑1 (MQ=255)
                                                                                agagAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGACAAATTGAGCGTTTACGCACGTTATACCCGTCGTgg                                                             >  1:577439/1‑143 (MQ=255)
                                                                                    aa‑aaaCTGCTGCGTTACCTGGTATCATTACGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATATTTAATGACCTGTTACGCTTCTGCCCGCCAGACAAAATGACAGTATACGGACGCTATACGCCACGTGGCgct                                                         >  2:348569/1‑141 (MQ=255)
                                                                                       aaaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTCACTGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCCCGTTATACCCGCCGTGGCggccc                                                       >  2:909871/1‑140 (MQ=255)
                                                                                             ctgcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACAt                                                 <  1:334962/143‑1 (MQ=255)
                                                                                             ctgcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCAGCTTTAATGACCTGTTGCGCTTCTGCCAGCCAGAAAAATTCAGCGTTTACGCACGTTCTACCCGTCGTGGCGGGCGGGACAt                                                 >  2:644128/1‑143 (MQ=255)
                                                                                                  ttACCTGGTCTCATTCCGTCATCACAACGAGTGCCACGAACAGGGCGTGGAACGCATCTTTAATCACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAAcc                                            <  1:323544/143‑1 (MQ=255)
                                                                                                      cgcctctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTg                                        <  1:244997/139‑1 (MQ=255)
                                                                                                         gTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGgcg                                     <  1:24416/143‑1 (MQ=255)
                                                                                                              aTTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTa                                >  1:676087/1‑143 (MQ=255)
                                                                                                               ttCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTAACGCTTCTGCCAGCCAGAAAAATTCAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCGGGGCATTAACCCGTGGCGCAGTaa                               >  2:450568/1‑143 (MQ=255)
                                                                                                                  cGTCATCACAACGGGTTCCACGAACAGTGCGTGGAACGCATCTTTACTGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAg                            <  2:664280/143‑1 (MQ=255)
                                                                                                                          caACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACACACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTg                    <  1:297091/143‑1 (MQ=255)
                                                                                                                                gTTCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGATAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCAt              >  1:802806/1‑143 (MQ=255)
                                                                                                                                 ttCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGCGCGCTTATAAAAGATTGAGGTGTCGTGCCAGCCAAGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCAtg             >  2:190389/1‑142 (MQ=255)
                                                                                                                                 ttCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCAGTGGCGCAGTCATCGAGGTTTTGTcccccg             >  2:560186/1‑140 (MQ=255)
                                                                                                                                    cACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCCGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATATTGTCGCCTCGAc          >  2:531706/1‑143 (MQ=255)
                                                                                                                                     aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGAcc         >  2:251666/1‑143 (MQ=255)
                                                                                                                                         aCAGTGCGTGGGACGCATCGTTAATGCCCTGTTACGCTTCTGCCAGCCAGTAAAACTGCGCGCTTCCGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACaa     <  2:889152/143‑1 (MQ=255)
                                                                                                                                            gTGCGTGGCGCGCATCTGTGATGACCTGTTACGCGTCTGCCAGCCAGAGAGATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGAc  <  2:511668/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
TGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGAC  >  NC_000913/2925901‑2926182

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