Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,108,996 T→C 19.9% L165L (TTA→CTA)  potA → spermidine/putrescine 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,108,9960TC19.9% 38.1 / 2.4 20L165L (TTA→CTA) potAspermidine/putrescine ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base C (2/2);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCAGGTTATGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGATG  >  CP000730/1108876‑1109134
                                                                                                                        |                                                                                                                                          
tCAGGTTATGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCCATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACt                                                                                                                        <  2:501/141‑1 (MQ=255)
    gTTATGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGaaa                                                                                                                    >  1:93354/1‑141 (MQ=255)
    gTTATGAAAAAAGAAATATTAATGAAATGACTAGCTGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAAGATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGACCTGata                                                                                                                    >  1:77407/1‑139 (MQ=255)
       aTGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCCACGTGTTGCAAATGCACGTGCGATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGc                                                                                                                 >  2:57004/1‑141 (MQ=255)
       aTGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCGGAAATATTATTGTTAGATGAAGCTTTATCCGCACTAGATTTGAAATTGCGTACTGAAATGc                                                                                                                 >  2:109275/1‑141 (MQ=255)
       aTGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCACTAGATTTGAAATTGCGTACTGAAATGc                                                                                                                 >  2:109280/1‑141 (MQ=255)
         gAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGATATTGCGTACTGAAATGCaa                                                                                                               <  1:267522/141‑1 (MQ=255)
                   atatTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTAc                                                                                                     <  2:194392/141‑1 (MQ=255)
                     atTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATGACga                                                                                                   >  1:221982/1‑141 (MQ=255)
                                 aGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTACTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATACGCATTAGATTTGAAACTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAAACt                                                                                       >  2:103782/1‑141 (MQ=255)
                                  gTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTa                                                                                      <  2:200151/141‑1 (MQ=255)
                                         aCAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAgg                                                                               <  1:252002/141‑1 (MQ=255)
                                          cAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGCACCAGATGTATTAGTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGt                                                                              <  1:222111/141‑1 (MQ=255)
                                                    gttttGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTta                                                                    <  1:11738/138‑1 (MQ=255)
                                                                                                   ttAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGActtt                     >  1:259775/1‑141 (MQ=255)
                                                                                                                  tCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTaaaa      >  1:232223/1‑141 (MQ=255)
                                                                                                                    cGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGa    <  1:57004/141‑1 (MQ=255)
                                                                                                                    cGCACTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGa    <  1:109275/141‑1 (MQ=255)
                                                                                                                    cGCACTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGa    <  1:109280/141‑1 (MQ=255)
                                                                                                                      cATTAGATTTGAAATTGCGGACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGATg  <  1:71855/141‑1 (MQ=255)
                                                                                                                        |                                                                                                                                          
TCAGGTTATGAAAAAAGAAATATTAATGAAATGAGTGGCGGACAAAAGCAACGTGTTGCAATTGCACGTGCTATCGTAAATGAACCAGAAATATTATTGTTAGATGAATCTTTATCCGCATTAGATTTGAAATTGCGTACTGAAATGCAATATGAATTACGAGAATTGCAATCTAGATTAGGTATTACATTTATATTTGTAACACATGATCAAGAAGAAGCGTTAGCATTAAGTGACTTTCTTTTTGTATTAAAAGATG  >  CP000730/1108876‑1109134

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