Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,005,859 G→A 20.4% G147E (GGG→GAG)  uup → putative transporter subunit of ABC superfamily: ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,005,8590GA20.4% 115.3 / 26.7 54G147E (GGG→GAG) uupputative transporter subunit of ABC superfamily: ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (7/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

AAACGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGG  >  CP009273/1005714‑1006001
                                                                                                                                                 |                                                                                                                                              
aaaCGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTgggg                                                                                                                                               <  2:11958/147‑1 (MQ=255)
aaaCGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGtt                                                                                                                                             <  1:1754/149‑1 (MQ=255)
aaaCGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTa                                                                                                                                            <  2:284340/150‑1 (MQ=255)
aaaCGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTa                                                                                                                                            <  1:631617/150‑1 (MQ=255)
  aCGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGt                                                                                                                                              >  1:448843/1‑146 (MQ=255)
   cGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGAt                                                                                                                                         <  1:854064/150‑1 (MQ=255)
    gCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATc                                                                                                                                        <  1:425331/150‑1 (MQ=255)
          aCGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAAc                                                                                                                                   <  1:228468/149‑1 (MQ=255)
          aCGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACg                                                                                                                                  <  1:692994/150‑1 (MQ=255)
           cggtatTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGTGCAACTGGAGTTAGATCCTAACGt                                                                                                                                 <  1:388129/147‑1 (MQ=255)
           cGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTAGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGt                                                                                                                                 >  1:381559/1‑150 (MQ=255)
             atatTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTg                                                                                                                               >  2:816777/1‑150 (MQ=255)
                     cgcCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTgtcg                                                                                                                       <  2:845077/150‑1 (MQ=255)
                       ccTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTgtcgtc                                                                                                                     <  2:836785/150‑1 (MQ=255)
                       ccTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTgtcgtc                                                                                                                     >  1:902770/1‑150 (MQ=255)
                        cTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGtcgtcg                                                                                                                    >  2:118904/1‑150 (MQ=255)
                         tGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGc                                                                                                                   >  2:888114/1‑150 (MQ=255)
                           gTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCtt                                                                                                                 >  1:925246/1‑150 (MQ=255)
                              atgaACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTcc                                                                                                              <  2:485101/150‑1 (MQ=255)
                                   cgccccGAACGAGACAAATCTCAACGAACTGGCTAAGGTTCAGGAACAGCTGGATCACCACAACCTATGGCAACTGGAAAACCGCATCAACGAATTACTGACCGAACTACGGCTAGTTCATATCTTTGCGCTTTCGTCGCTTTCcgtccg                                                                                                         >  2:623564/4‑146 (MQ=255)
                                   cGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCcggcgg                                                                                                         <  2:694057/150‑1 (MQ=255)
                                      cccGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTg                                                                                                      <  1:622409/150‑1 (MQ=255)
                                                aaaaaTCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTa                                                                                              >  1:856352/1‑148 (MQ=255)
                                                   aaTCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCg                                                                                         >  2:254195/1‑150 (MQ=255)
                                                       tcAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCAt                                                                                     >  2:966402/1‑150 (MQ=255)
                                                         aaCGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCAtt                                                                                    <  1:863483/149‑1 (MQ=255)
                                                            gAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGa                                                                                >  1:960340/1‑150 (MQ=255)
                                                                  gCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGAcgcg                                                                            <  2:309220/148‑1 (MQ=255)
                                                                       ggTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGt                                                                     >  2:102369/1‑150 (MQ=255)
                                                                         ttCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGt                                                                     >  1:265747/1‑148 (MQ=255)
                                                                                 cAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCcg                                                           >  2:159801/1‑150 (MQ=255)
                                                                                    cTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCcgcgc                                                        >  1:719890/1‑150 (MQ=255)
                                                                                         tCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCt                                                   <  1:102369/150‑1 (MQ=255)
                                                                                             cacaACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTg                                               <  1:888114/150‑1 (MQ=255)
                                                                                                   cTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTtgat                                         <  1:854485/150‑1 (MQ=255)
                                                                                                     gtgGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTtgatg                                        >  2:733144/1‑149 (MQ=255)
                                                                                                      tgGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTtgatga                                       <  1:276854/149‑1 (MQ=255)
                                                                                                         cAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAAcc                                    >  2:185744/1‑149 (MQ=255)
                                                                                                          aGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGa                                  >  1:821233/1‑150 (MQ=255)
                                                                                                          aGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGa                                  >  1:471243/1‑150 (MQ=255)
                                                                                                           gCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGa                                  >  1:821580/1‑149 (MQ=255)
                                                                                                           gCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGa                                  >  2:256531/1‑149 (MQ=255)
                                                                                                            cTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACa                                <  2:319435/150‑1 (MQ=255)
                                                                                                              ggAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACaaa                              >  2:259246/1‑150 (MQ=255)
                                                                                                               gAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAac                             <  2:799010/150‑1 (MQ=255)
                                                                                                               gAAAACCGCATCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAac                             >  2:207894/1‑150 (MQ=255)
                                                                                                                   aCCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAaccacc                         >  1:405304/1‑150 (MQ=255)
                                                                                                                        aTCAACGAAGTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGa                     >  1:382333/1‑149 (MQ=255)
                                                                                                                              gAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGaa              <  1:146443/150‑1 (MQ=255)
                                                                                                                                 gTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAAc            >  1:922259/1‑149 (MQ=255)
                                                                                                                                 gTGCTGGCGCAACTGGAGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAAcc           <  2:448843/150‑1 (MQ=255)
                                                                                                                                  tGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCa          <  2:922259/150‑1 (MQ=255)
                                                                                                                                      ggCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGa      >  1:834311/1‑150 (MQ=255)
                                                                                                                                       gcgcAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCg       <  1:809784/148‑1 (MQ=255)
                                                                                                                                          cAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGActgg  <  2:378307/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AAACGCTATCACGATATTTCGCGCCTGGTGATGAACGACCCGAGCGAGAAAAATCTCAACGAACTGGCGAAGGTTCAGGAACAGCTGGATCACCACAACCTGTGGCAGCTGGAAAACCGCATCAACGAAGTGCTGGCGCAACTGGGGTTAGATCCTAACGTTGCGCTGTCGTCGCTTTCCGGCGGCTGGTTGCGTAAAGCGGCATTAGGACGCGCGCTGGTGAGTAATCCGCGCGTGCTGTTGCTTGATGAACCGACAAACCACCTGGATATTGAAACCATCGACTGG  >  CP009273/1005714‑1006001

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