Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,018,582 T→A 9.8% K113* (AAA→TAA)  mgt ← transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,018,5820TA 9.8% 128.2 / 7.9 41K113* (AAA→TAA) mgttransglycosylase
Reads supporting (aligned to +/- strand):  ref base T (15/22);  new base A (2/2);  total (17/24)
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.52e-01

AATTCTTTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGT  >  CP000730/2018449‑2018721
                                                                                                                                     |                                                                                                                                           
aatTCTTTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCTAATCg                                                                                                                                      <  2:345988/141‑1 (MQ=255)
  ttCTTTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGaa                                                                                                                                    <  1:266112/141‑1 (MQ=255)
       ttACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCat                                                                                                                               >  1:104674/1‑141 (MQ=255)
       ttACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCat                                                                                                                               >  2:410994/1‑141 (MQ=255)
       ttACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCat                                                                                                                               >  2:45032/1‑141 (MQ=255)
                        gAACGATCCTTATCATAAAAATAATTTTTGACGACTTGGTGTGTAATGGTACTGCCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGt                                                                                                              <  2:364255/141‑1 (MQ=255)
                              tCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCt                                                                                                        >  1:347517/1‑141 (MQ=255)
                                              aatTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAATTCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAgg                                                                                        <  1:270213/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  1:419565/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  2:39979/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  2:100764/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  2:52920/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTACAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  2:414170/141‑1 (MQ=255)
                                                tttttGACAACTTGTTGTGTAATGGTACTACCACCTTGCACACCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCa                                                                                      <  2:438921/141‑1 (MQ=255)
                                                     gACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCttt                                                                                 >  1:163748/1‑141 (MQ=255)
                                                          cTTGTTGTGTAATGGTACTACCACCTTGCACATGTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACAt                                                                            <  2:115938/141‑1 (MQ=255)
                                                           ttgttgTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATa                                                                           <  2:221493/141‑1 (MQ=255)
                                                            tgttgtGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATAc                                                                          <  1:286841/141‑1 (MQ=255)
                                                              ttgtGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATActc                                                                        >  1:380842/1‑141 (MQ=255)
                                                              ttgtGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATActc                                                                        >  1:8761/1‑141 (MQ=255)
                                                              ttgtGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATActc                                                                        >  1:224588/1‑141 (MQ=255)
                                                                gtgtAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTg                                                                      <  1:331155/141‑1 (MQ=255)
                                                                         tactacCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTAt                                                             >  2:250818/1‑141 (MQ=255)
                                                                               caccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTg                                                       >  1:455350/1‑141 (MQ=255)
                                                                               caccTTGCACATCTCTGTCGCTAATCGATGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTg                                                       >  1:127144/1‑141 (MQ=255)
                                                                                       acaTCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTACCATACTCTGGCATGTTATCAGCTGACACaaaa                                               >  1:413459/1‑141 (MQ=255)
                                                                                       acaTCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACaaaa                                               >  1:371658/1‑141 (MQ=255)
                                                                                                           tGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAAtttt                           <  1:320011/141‑1 (MQ=255)
                                                                                                              aaaTAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCt                        <  2:3257/141‑1 (MQ=255)
                                                                                                                 tAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAg                     <  2:119073/141‑1 (MQ=255)
                                                                                                                 tAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAg                     <  2:366933/141‑1 (MQ=255)
                                                                                                                  aaaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGt                    <  2:8761/141‑1 (MQ=255)
                                                                                                                  aaaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGt                    <  2:224588/141‑1 (MQ=255)
                                                                                                                  aaaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGt                    <  2:380842/141‑1 (MQ=255)
                                                                                                                   aaGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGtt                   <  2:5836/141‑1 (MQ=255)
                                                                                                                     gCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCa                 >  2:220717/1‑141 (MQ=255)
                                                                                                                           gTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCcaca           >  1:441916/1‑141 (MQ=255)
                                                                                                                            ttGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACAt          <  2:159240/141‑1 (MQ=255)
                                                                                                                              gTACCTTACAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTa        <  2:373820/141‑1 (MQ=255)
                                                                                                                                aCCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATc      >  1:250088/1‑141 (MQ=255)
                                                                                                                                aCCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATc      >  2:287110/1‑141 (MQ=255)
                                                                                                                                   tttCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATcgcg   >  1:355419/1‑141 (MQ=255)
                                                                                                                                   tttCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATcgcg   >  1:251358/1‑141 (MQ=255)
                                                                                                                                    ttCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCgt  >  1:312759/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                           
AATTCTTTTACTTTTCTAGTAAATGAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGT  >  CP000730/2018449‑2018721

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