Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 818,656 G→A R29R (CGG→CGA moaC → cyclic pyranopterin monophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913818,6560GA95.5% 134.4 / ‑1.8 44R29R (CGG→CGAmoaCcyclic pyranopterin monophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (22/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GTACGCGTCCGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATC  >  NC_000913/818527‑818788
                                                                                                                                 |                                                                                                                                    
gTACGCGTCCGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCt                                                                                                                            <  2:183680/140‑1 (MQ=255)
gTACGCGTCCGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCt                                                                                                                            >  1:1024527/1‑140 (MQ=255)
 tACGCGTCCGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCtt                                                                                                                           >  1:98824/1‑140 (MQ=255)
        ccGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCAcc                                                                                                                    >  2:54362/1‑140 (MQ=255)
                      tccaCATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGAc                                                                                                      >  1:641207/1‑140 (MQ=255)
                       ccaCATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACg                                                                                                     <  1:54362/140‑1 (MQ=255)
                             ttGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCg                                                                                               <  2:641207/140‑1 (MQ=255)
                             ttGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCg                                                                                               <  1:368531/140‑1 (MQ=255)
                                     aTAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACg                                                                                                     <  2:702712/126‑1 (MQ=255)
                                     aTAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACg                                                                                                     >  1:702712/1‑126 (MQ=255)
                                          tATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATg                                                                                  >  1:707128/1‑140 (MQ=255)
                                          tATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATg                                                                                  <  2:222692/140‑1 (MQ=255)
                                             gTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTc                                                                               >  1:371158/1‑140 (MQ=255)
                                             gTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTc                                                                               >  2:891124/1‑140 (MQ=255)
                                             gTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTc                                                                               >  2:502679/1‑140 (MQ=255)
                                               cGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGc                                                                             >  1:139590/1‑140 (MQ=255)
                                                    cTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGccacca                                                                        >  2:317571/1‑140 (MQ=255)
                                                       aCCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACaa                                                                     >  2:788247/1‑140 (MQ=255)
                                                       aCCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACaa                                                                     >  2:321985/1‑140 (MQ=255)
                                                       aCCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACaa                                                                     >  2:149568/1‑140 (MQ=255)
                                                            tatCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCg                                                                <  2:1024527/140‑1 (MQ=255)
                                                            tatCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCg                                                                <  2:195861/140‑1 (MQ=255)
                                                            tatCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCg                                                                >  1:195861/1‑140 (MQ=255)
                                                              tCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGAc                                                              <  1:149568/140‑1 (MQ=255)
                                                              tCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGAc                                                              <  1:788247/140‑1 (MQ=255)
                                                              tCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGAc                                                              <  2:522717/140‑1 (MQ=255)
                                                                       cTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGcc                                                     <  1:502679/140‑1 (MQ=255)
                                                                        tGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCa                                                    <  2:98824/140‑1 (MQ=255)
                                                                             aaGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTgcg                                               <  1:1209436/140‑1 (MQ=255)
                                                                               gcgcACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGcgcg                                             >  1:84570/1‑140 (MQ=255)
                                                                                cgcACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGt                                            >  1:705528/1‑140 (MQ=255)
                                                                                 gcACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTa                                           <  1:891124/140‑1 (MQ=255)
                                                                                     aTGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGc                                       <  2:930784/140‑1 (MQ=255)
                                                                                               tctcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAg                             >  1:1122040/1‑140 (MQ=255)
                                                                                                 tcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAggc                           >  1:159225/1‑140 (MQ=255)
                                                                                                 tcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAggc                           >  2:239762/1‑140 (MQ=255)
                                                                                                   cGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAggcgg                         >  1:965416/1‑140 (MQ=255)
                                                                                                   cGCCAAAGCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAggcgg                         >  1:561046/1‑140 (MQ=255)
                                                                                                          gCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACg                  <  2:139590/140‑1 (MQ=255)
                                                                                                          gCGGAAACCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACg                  <  1:239762/140‑1 (MQ=255)
                                                                                                                aCCGTGCGTGAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTg            <  2:707128/140‑1 (MQ=255)
                                                                                                                    tgcgtgAAGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGAt        >  1:1086317/1‑140 (MQ=255)
                                                                                                                        tgAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGa    <  1:317571/140‑1 (MQ=255)
                                                                                                                          aaGCGCGAGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATc  <  2:705528/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GTACGCGTCCGTGTAATTTCCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATC  >  NC_000913/818527‑818788

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