Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,134,811 Δ1 bp coding (681/693 nt) yghT → putative ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,134,8110T.97.1% 133.6 / ‑2.9 34coding (681/693 nt)yghTputative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (13/20);  minor base G (1/0);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGGAAGATTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTA  >  NC_000913/3134680‑3134930
                                                                                                                                   |                                                                                                                       
gCGAGGAAGATTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCa                                                                                                                 <  2:42367/139‑1 (MQ=255)
          ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCaaaa                                                                                                       <  1:558021/139‑1 (MQ=255)
          ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCaaaa                                                                                                       <  2:20540/139‑1 (MQ=255)
           tGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCaaaaa                                                                                                      >  2:40458/1‑139 (MQ=255)
                tAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGaa                                                                                                 >  1:451904/1‑139 (MQ=255)
                       acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                          <  2:363120/139‑1 (MQ=255)
                       acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                          <  2:417006/139‑1 (MQ=255)
                             aaCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGcc                                                                                    >  2:187250/1‑139 (MQ=255)
                              aCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCa                                                                                   <  2:334424/139‑1 (MQ=255)
                                   tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAcc                                                                              >  1:561894/1‑139 (MQ=255)
                                                cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                 <  2:451904/139‑1 (MQ=255)
                                                   gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCa                                                                                                                 <  2:495663/88‑1 (MQ=255)
                                                   gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCa                                                                                                                 >  1:495663/1‑88 (MQ=255)
                                                   gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAg                                                              >  2:484696/1‑139 (MQ=255)
                                                                ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                 <  2:741472/139‑1 (MQ=255)
                                                                     cAGCCACTGTAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                            <  1:660398/139‑1 (MQ=255)
                                                                     cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                            <  1:187250/139‑1 (MQ=255)
                                                                     cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                            <  2:461377/139‑1 (MQ=255)
                                                                        ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                               <  2:242978/101‑1 (MQ=255)
                                                                        ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                               >  1:242978/1‑101 (MQ=255)
                                                                               aaCAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTa                                  <  1:703186/139‑1 (MQ=255)
                                                                                      gttggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCcca                           >  2:214563/1‑139 (MQ=255)
                                                                                         ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATgg                                                                                                   <  2:319306/64‑1 (MQ=255)
                                                                                         ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATgg                                                                                                   >  1:319306/1‑64 (MQ=255)
                                                                                              aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                   >  2:246022/1‑139 (MQ=255)
                                                                                              aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                   >  2:782730/1‑139 (MQ=255)
                                                                                              aTGCAGAAAAAGCGATTACGCCTTTTATGTCCGCGGGGGTTATCACTAATCAAAAATGGAAATGATCGATCGCCAGGACCGGGCATTTTCAGGAAGCTTAAATCAACTGCAATGCTATCCACTACAGCACCCCagaaag                   >  1:36818/1‑139 (MQ=255)
                                                                                               tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaaga                  >  1:487133/1‑139 (MQ=255)
                                                                                                        aGCGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgcc         <  2:487133/139‑1 (MQ=255)
                                                                                                          cGATTACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccgc       <  1:782730/139‑1 (MQ=255)
                                                                                                              tACGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgt   <  2:196232/139‑1 (MQ=255)
                                                                                                               aCGGATTTTATGACCGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgta  <  1:484696/139‑1 (MQ=255)
                                                                                                                            ccGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                       <  1:529464/89‑1 (MQ=255)
                                                                                                                            ccGCGCGGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                       >  2:529464/1‑89 (MQ=255)
                                                                                                                                   |                                                                                                                       
GCGAGGAAGATTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTA  >  NC_000913/3134680‑3134930

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