Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,036 A→C 8.8% E230A (GAA→GCA)  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,0360AC 8.8% 108.7 / 3.0 45E230A (GAA→GCA) queF7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base C (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

CCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACC  >  NC_000913/2925896‑2926175
                                                                                                                                             |                                                                                                                                           
ccACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCGAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAAc                                                                                                                                           >  1:768086/1‑143 (MQ=255)
     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:407238/1‑143 (MQ=255)
                              ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                             <  1:634656/143‑1 (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)
                                                                                                                                             |                                                                                                                                           
CCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAA‑AAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACC  >  NC_000913/2925896‑2926175

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