Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,263,602 C→T V210V (GTG→GTA qorA ← putative quinone oxidoreductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,263,6020CT97.6% 124.0 / ‑4.3 42V210V (GTG→GTAqorAputative quinone oxidoreductase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/20);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

TAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTCGCG  >  NC_000913/4263482‑4263737
                                                                                                                        |                                                                                                                                       
taggTTCACCCCTGTAACCGCGCCTGATGAGTTCCCAAAACTGCCCATTAAGCCGCGGCGTTACAGGCAACCCAGCGTCCGTTCCCAGGTGTTTCTGCCCACGGAATCTTACAGCACGCTTACTTACTTACCGCCGGTGa                                                                                                                      <  2:142100/137‑1 (MQ=255)
tAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGa                                                                                                                      <  2:681902/140‑1 (MQ=255)
        cacCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCGACGGAATCGTACGCCATGCGTACTTGCTTAGCGCCGGTGATCTCtttt                                                                                                              >  2:262846/1‑140 (MQ=255)
        cacCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCtttt                                                                                                              >  2:542717/1‑140 (MQ=255)
             gTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCg                                                                                                         >  1:162809/1‑140 (MQ=255)
              tAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCAAGCGCACTTTCTCACCGCCGGTCATCTCTTTTAACCGc                                                                                                        >  1:314864/1‑140 (MQ=255)
               caCCCCACCTGATGCGTTCCCAAAATGGACCTTTAGGCCGGGGCGGTGCGGGCCATCCAGCGACCCTCCCCAGGTGTCCCTTCCCCCGGCATCGTCCCACACGGGTCCTTTCCGCCCGCCGGTGATCTCTTTTAGCCgat                                                                                                       <  2:218316/139‑3 (MQ=255)
               aaCCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCt                                                                                                       <  1:626456/140‑1 (MQ=255)
               aaCCGCACCTGATGACTTGCCACAAATGACCATTAAGTCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCt                                                                                                       <  1:240287/140‑1 (MQ=255)
                   gCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGAc                                                                                                   >  1:576671/1‑140 (MQ=255)
                       cTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGa                                                                                               >  1:551772/1‑140 (MQ=255)
                              gTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCtt                                                                                        <  1:474167/140‑1 (MQ=255)
                                tGCCAAAACTGCCCATTAGGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGCCTCTGCCCACGGAATCGTCCACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCa                                                                                      <  2:508463/140‑1 (MQ=255)
                                tGCCAAAACTGACCATTGAGCCGCGGCGTTGCAGGCAATGCAGCGACCCTTCCCAGGTGTCTCTGCCCCCGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCa                                                                                      <  2:238398/140‑1 (MQ=255)
                                tGCCAAAACTGACCATTAAGCCGCGGCGTTGGAGGCAAAACAGCGACCATCCCCAGGTGTCTCTGCCCACGGAATCGTGCACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCa                                                                                      <  2:306524/140‑1 (MQ=255)
                                 gCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTGTCTTACCGCCGGTGATATCTTTTAAACGCTTGACCAGATCCTCTTCAc                                                                                     >  2:378666/1‑140 (MQ=255)
                                  ccAAAACTGACCATTGAGCCGCGGCCTTGTAGGCAATACAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACg                                                                                    <  1:378666/140‑1 (MQ=255)
                                   cAAACCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTCCTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGa                                                                                   <  2:629397/140‑1 (MQ=255)
                                   cAAAACTGACCATTACGCCGCGTCGCTGCTGGCAATCTAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGa                                                                                   <  1:262846/140‑1 (MQ=255)
                                   cAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGa                                                                                   <  1:507586/140‑1 (MQ=255)
                                    aaaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGAt                                                                                  >  1:698314/1‑140 (MQ=255)
                                    aaaaCTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGAt                                                                                  <  1:542717/140‑1 (MQ=255)
                                            ccATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTtaata                                                                          >  2:186074/1‑140 (MQ=255)
                                                       cggcgTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgc                                                               >  1:531953/1‑140 (MQ=255)
                                                       cggcgTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGAACACCACGCGTACTTTCTTACCGCCGGCGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCAcgc                                                               >  2:461982/1‑140 (MQ=255)
                                                          cgTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACCCCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGcc                                                            <  1:196978/140‑1 (MQ=255)
                                                            ttGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCg                                                          >  2:862644/1‑140 (MQ=255)
                                                               cAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCtt                                                       <  2:531953/140‑1 (MQ=255)
                                                                    aaTCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAgc                                                  <  2:647734/140‑1 (MQ=255)
                                                                         aGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGct                                             >  1:791820/1‑140 (MQ=255)
                                                                          gCGACCGTCCCGAGGTGTTTCTCCCCAAGGAATCGTACACCTGGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTAGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGctc                                            <  2:551772/140‑1 (MQ=255)
                                                                            gACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTg                                          >  1:531674/1‑140 (MQ=255)
                                                                                 ttCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCtt                                     >  1:184856/1‑140 (MQ=255)
                                                                                       ggTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcg                               <  2:531674/140‑1 (MQ=255)
                                                                                       ggTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcg                               <  1:862644/140‑1 (MQ=255)
                                                                                       ggTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcg                               <  2:576671/140‑1 (MQ=255)
                                                                                       ggTGTCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcg                               <  2:162809/140‑1 (MQ=255)
                                                                                        gtgtCTCTGCCCACGGAATCGTACCCCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCTACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcgc                              <  2:459174/140‑1 (MQ=255)
                                                                                        gtgtCTCTGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTgcgc                              <  2:114390/140‑1 (MQ=255)
                                                                                           tctcTGCCCACAGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGt                           >  1:769884/1‑140 (MQ=255)
                                                                                             tctGCCCACGGAATCGTACACCACGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTc                         >  1:453093/1‑140 (MQ=255)
                                                                                                                accacGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGttt      >  1:664003/1‑140 (MQ=255)
                                                                                                                 ccacGCGTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTc     <  2:664003/140‑1 (MQ=255)
                                                                                                                    cgcgTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGGTTTTTGCACGGCTCGTACGGTCCCGATAAGTTTcgcg  >  1:316170/1‑140 (MQ=255)
                                                                                                                    cgcgTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTcgcg  >  2:530934/1‑140 (MQ=255)
                                                                                                                    cgcgTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTcgcg  >  2:588502/1‑140 (MQ=255)
                                                                                                                    cgcgTACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTcgcg  >  2:761653/1‑140 (MQ=255)
                                                                                                                        |                                                                                                                                       
TAAGTTCACACCGGTAACCGCACCTGATGAGTTGCCAAAACTGACCATTAAGCCGCGGCGTTGCAGGCAATCCAGCGACCGTTCCCAGGTGTCTCTGCCCACGGAATCGTACACCACGCGCACTTTCTTACCGCCGGTGATCTCTTTTAACCGCTCGACCAGATCCTCTTCACGATAGTTAATAACCTGCCACGCGCCCGCTTTTAGCGCGCTCTGCGCTTTTTGCGCGGTTCCTACGGTGCCGATAAGTTTCGCG  >  NC_000913/4263482‑4263737

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