Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,039 A→T 8.9% Q231L (CAG→CTG)  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,0390AT 8.9% 118.7 / 5.1 45Q231L (CAG→CTG) queF7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01

TGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAA  >  NC_000913/2925901‑2926179
                                                                                                                                           |                                                                                                                                            
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)
                                                                                                                                           |                                                                                                                                            
TGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAA  >  NC_000913/2925901‑2926179

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: ‑