Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,416,458 A→T 100% I216N (ATT→AAT)  oppD ← oligopeptide ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,416,4580AT98.5% 226.9 / ‑3.7 65I216N (ATT→AAT) oppDoligopeptide ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (31/33);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTC  >  CP000730/1416324‑1416594
                                                                                                                                      |                                                                                                                                        
tAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACt                                                                                                                                    <  2:434891/141‑1 (MQ=255)
 aaCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTACTATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                                                                   <  2:28611/141‑1 (MQ=255)
 aaCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                                                                   >  2:662635/1‑141 (MQ=255)
 aaCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTg                                                                                                                                   >  1:23721/1‑141 (MQ=255)
  aCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCATTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGa                                                                                                                                  <  2:625168/141‑1 (MQ=255)
  aCTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGa                                                                                                                                  <  1:258956/141‑1 (MQ=255)
   cTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGAc                                                                                                                                 >  1:496582/1‑141 (MQ=255)
   cTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGAc                                                                                                                                 <  2:369153/141‑1 (MQ=255)
   cTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGAc                                                                                                                                 <  2:398686/141‑1 (MQ=255)
            aGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTGGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTCCa                                                                                                                        >  2:547311/1‑141 (MQ=255)
             gTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCaa                                                                                                                       <  2:398764/141‑1 (MQ=255)
                 tCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACa                                                                                                                   <  1:750054/141‑1 (MQ=255)
                   ccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACACTATTTCAATATATg                                                                                                                 >  2:432578/1‑141 (MQ=255)
                   ccTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATg                                                                                                                 <  2:614208/141‑1 (MQ=255)
                          aaTTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAg                                                                                                          >  2:730454/1‑141 (MQ=255)
                           aTTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                                                                         >  1:507947/1‑141 (MQ=255)
                              gtgtTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAt                                                                                                      >  1:107361/1‑141 (MQ=255)
                                   cTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgtgtg                                                                                                 >  1:316128/1‑141 (MQ=255)
                                      tatTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGata                                                                                              >  1:296140/1‑141 (MQ=255)
                                          tAATTTGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGAtaataa                                                                                          >  2:511662/1‑141 (MQ=255)
                                          tAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTATTATTACCTAGTTCAATGATTTGACCTGCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGAtaataa                                                                                          <  1:468993/141‑1 (MQ=255)
                                            aTTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                                                                        >  1:688698/1‑141 (MQ=255)
                                             ttAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaata                                                                                       >  2:650379/1‑141 (MQ=255)
                                               aGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAaataat                                                                                     <  1:470944/141‑1 (MQ=255)
                                                 aCGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGt                                                                                   >  2:517300/1‑141 (MQ=255)
                                                   gTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGttt                                                                                 >  2:635738/1‑141 (MQ=255)
                                                    ttCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTg                                                                                <  1:688261/141‑1 (MQ=255)
                                                      cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGAc                                                                              >  2:383241/1‑141 (MQ=255)
                                                       aGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGAcc                                                                             <  1:379314/141‑1 (MQ=255)
                                                         tataCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAt                                                                           <  2:311482/141‑1 (MQ=255)
                                                         tataCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAt                                                                           <  2:200441/141‑1 (MQ=255)
                                                                          tCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATAGAATATGc                                                          >  2:542294/1‑141 (MQ=255)
                                                                          tCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGc                                                          >  1:155262/1‑141 (MQ=255)
                                                                          tCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGc                                                          <  2:12856/141‑1 (MQ=255)
                                                                            aTACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTg                                                        <  2:472804/141‑1 (MQ=255)
                                                                             tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGt                                                       <  2:253083/141‑1 (MQ=255)
                                                                               cTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGtaa                                                     <  1:730454/141‑1 (MQ=255)
                                                                                   aaaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATaaa                                                 >  2:147127/1‑141 (MQ=255)
                                                                                   aaaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATaaa                                                 >  1:663663/1‑141 (MQ=255)
                                                                                     aTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATc                                               >  1:301374/1‑141 (MQ=255)
                                                                                      tGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATcc                                              <  2:496582/141‑1 (MQ=255)
                                                                                         ttAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATc                                           >  2:457741/1‑141 (MQ=255)
                                                                                         ttAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATc                                           >  1:686558/1‑141 (MQ=255)
                                                                                           aaTATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCAc                                         <  1:165420/141‑1 (MQ=255)
                                                                                                tACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGAt                                    <  2:502571/141‑1 (MQ=255)
                                                                                                 aCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGAtt                                   >  2:651733/1‑141 (MQ=255)
                                                                                                 aCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGAtt                                   <  1:529492/141‑1 (MQ=255)
                                                                                                  ccTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTa                                  <  1:255320/141‑1 (MQ=255)
                                                                                                           aTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAAGATGCTGTAATAAATCCATCACTTGATTACCGTTCaaa                         <  1:568242/141‑1 (MQ=255)
                                                                                                               tttGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCAt                     >  1:702905/1‑141 (MQ=255)
                                                                                                                   aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCaa                 >  2:630559/1‑141 (MQ=255)
                                                                                                                    ccTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAat                >  2:493999/1‑141 (MQ=255)
                                                                                                                    ccTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAat                <  1:130884/141‑1 (MQ=255)
                                                                                                                     cttcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAatg               <  1:280096/141‑1 (MQ=255)
                                                                                                                           ttAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCaa         >  2:707342/1‑141 (MQ=255)
                                                                                                                             aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACt       <  1:147127/141‑1 (MQ=255)
                                                                                                                             aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACt       <  2:409091/141‑1 (MQ=255)
                                                                                                                              aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTTTATGCATTCATTGATGCAACTg      >  2:231438/1‑141 (MQ=255)
                                                                                                                              aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTg      >  2:757392/1‑141 (MQ=255)
                                                                                                                              aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTg      >  1:434496/1‑141 (MQ=255)
                                                                                                                                aCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGt    <  1:168602/141‑1 (MQ=255)
                                                                                                                                 cATAATTGTCCTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGtt   >  1:690936/1‑141 (MQ=255)
                                                                                                                                 cATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGtt   <  1:307037/141‑1 (MQ=255)
                                                                                                                                 cATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGtt   <  1:87828/141‑1 (MQ=255)
                                                                                                                                  aTAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTc  <  2:413862/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                        
TAACTCAATCATAGTAATCCCTCTTTAATTGTGTTCTATATTTAATTAGACGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAATAAATCCATCACTTGATTACCGTTCAAAGCATCCAATGATGCAACTGGTTC  >  CP000730/1416324‑1416594

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