Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 988,184 T→C 8.5% M30T (ATG→ACG)  oppB1 → oligopeptide ABC superfamily ATP binding cassette transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730988,1840TC 8.5% 156.1 / 4.0 57M30T (ATG→ACG) oppB1oligopeptide ABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/21);  new base C (2/3);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAAGGGAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAGGCG‑ATTTTGGTAA  >  CP000730/988044‑988322
                                                                                                                                            |                                                                                                                                           
tttGAAGGGAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAAt                                                                                                                                             <  2:144366/141‑1 (MQ=255)
   gAAGGGAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGcc                                                                                                                                          <  1:474157/141‑1 (MQ=255)
     aGGGAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAg                                                                                                                                        >  1:362427/1‑141 (MQ=255)
       ggAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGt                                                                                                                                      <  1:263176/141‑1 (MQ=255)
                 tttAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTa                                                                                                                            >  1:708667/1‑141 (MQ=255)
                     gAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGa                                                                                                                        <  1:155851/141‑1 (MQ=255)
                        aTCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGc                                                                                                                     <  1:597622/141‑1 (MQ=255)
                          cGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCta                                                                                                                   >  1:623228/1‑141 (MQ=255)
                           ggAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCtaa                                                                                                                  <  1:619177/141‑1 (MQ=255)
                                ttGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCtaaattaa                                                                                                             >  1:625058/1‑141 (MQ=255)
                                          aggGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACa                                                                                                   >  1:527765/1‑141 (MQ=255)
                                                ttAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGa                                                                                             <  1:680156/141‑1 (MQ=255)
                                                           aaTATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTTTTTTTTTTTTTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATg                                                                                  >  2:229432/1‑141 (MQ=255)
                                                            atatatTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTTTTTTTTTTTTTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGa                                                                                 >  2:196409/1‑141 (MQ=255)
                                                             tatatTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGaa                                                                                >  1:554508/1‑141 (MQ=255)
                                                             tatatTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGaa                                                                                >  1:360251/1‑141 (MQ=255)
                                                              atatTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTGGCCATTTAACGATGCTAAATTAAATGCTTAACAAAAAGAAATTTTAAATGaaa                                                                               <  1:523906/141‑1 (MQ=255)
                                                              atatTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGaaa                                                                               <  1:523969/141‑1 (MQ=255)
                                                                       aaagttttatttaTATTCTTATTTCTTTATTTATTATTATTACAATTACCTTTTTCTTAATGAAATTAACGCCAGGTTAGCCATTTAGCGATGATACATTAAAAAAAAAAAAAAAAGAAATTTTAAATGAAAAATATGGAt                                                                      <  2:136020/136‑1 (MQ=255)
                                                                       aaCGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCAGTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGAt                                                                      >  1:673592/1‑141 (MQ=255)
                                                                        aCGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGAtt                                                                     >  1:597041/1‑141 (MQ=255)
                                                                          gATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAGAAATATGGATTaa                                                                   >  2:592297/1‑141 (MQ=255)
                                                                          gATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTaa                                                                   >  2:169022/1‑141 (MQ=255)
                                                                           atttatttATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTaaa                                                                  <  1:169022/141‑1 (MQ=255)
                                                                               atttatATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCt                                                                                                                    >  1:451027/1‑87 (MQ=255)
                                                                               atttatATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCt                                                                                                                    <  2:451027/87‑1 (MQ=255)
                                                                               atttatATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTCAATGAt                                                              <  2:584339/141‑1 (MQ=255)
                                                                               aattatATGCATATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTATATTAAATGCTGAATAAAAAGAAAATATAAAAAAAAAATATGGATTAAATGAt                                                              <  2:410803/139‑1 (MQ=255)
                                                                                 ttatATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATcc                                                            <  1:495440/141‑1 (MQ=255)
                                                                                   atatGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTg                                                          >  1:209692/1‑141 (MQ=255)
                                                                                    tatGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGt                                                         >  2:104404/1‑141 (MQ=255)
                                                                                          tATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTAc                                                   <  2:298348/141‑1 (MQ=255)
                                                                                            tttCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGc                                                 >  1:104753/1‑141 (MQ=255)
                                                                                            tttCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGc                                                 >  1:254107/1‑141 (MQ=255)
                                                                                            tttCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGc                                                 <  2:623228/141‑1 (MQ=255)
                                                                                                tttatttaTTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTa                                             >  1:13130/1‑141 (MQ=255)
                                                                                                          attattaCAATTACATTTTTCTTAATGTATTTATTGCCAGGTTCGCCATTTAACGATGCTAAATTTAATGCTGAACAAAAAGAAATTTTAAATAAAAAATATGAATTAAATGATCCTGTAGCTACGCAGTATTTACattat                                   >  2:385950/1‑141 (MQ=255)
                                                                                                          attattaCAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAATGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACaatat                                   >  2:439961/1‑141 (MQ=255)
                                                                                                           tttttataattaaaTTTTTCTTAATGAAATTAATGCCAGGTTTGCCATTTACCGATGCTAAATTAAATGCTGAAAAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACattatt                                  <  2:324907/138‑1 (MQ=255)
                                                                                                             attacaattacaTTTTTCTTTTTTTATTTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTa                                >  2:397872/1‑141 (MQ=255)
                                                                                                               tacaattataTTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAAAGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTaaa                              <  2:87650/141‑1 (MQ=255)
                                                                                                                   attacaTTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAAt                          <  1:576515/141‑1 (MQ=255)
                                                                                                                    ttacaTTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAAtg                         >  1:657036/1‑141 (MQ=255)
                                                                                                                      acaTTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAAtgtt                       >  2:68224/1‑141 (MQ=255)
                                                                                                                          ttttCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCtaa                                                                                                                  >  2:85459/1‑46 (MQ=255)
                                                                                                                          ttttCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCtaa                                                                                                                  <  1:85459/46‑1 (MQ=255)
                                                                                                                           tttCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTAc                  >  1:319424/1‑141 (MQ=255)
                                                                                                                             tttttttttttttAATGCCAGATTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATTGATTAAATGATGCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAg                >  2:29435/12‑141 (MQ=255)
                                                                                                                             tCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAg                >  2:553185/1‑141 (MQ=255)
                                                                                                                              cTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAgg               <  1:389793/141‑1 (MQ=255)
                                                                                                                              attAATGAAATAAATACAAAAATCGATATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATAAAAAAAAAAAAAAAAATGTTGTTACAgg               <  2:362427/140‑1 (MQ=255)
                                                                                                                                   tttttttAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAGGCG‑Att         >  2:119051/6‑141 (MQ=255)
                                                                                                                                       aTTAACGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTaaa                              >  2:646072/1‑117 (MQ=255)
                                                                                                                                       aTTAACGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTaaa                              >  2:579259/1‑117 (MQ=255)
                                                                                                                                       aTTAACGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTaaa                              <  1:646072/117‑1 (MQ=255)
                                                                                                                                       aTTAACGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAAGATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTaaa                              <  1:579259/117‑1 (MQ=255)
                                                                                                                                        ttAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAGGCG‑ATTTTGGt    >  1:517829/1‑141 (MQ=255)
                                                                                                                                           aTGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAGGCGCTTTTTGGTaa  >  1:129469/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TTTGAAGGGAATGGAAATTTAGAAATCGGAATTTGTTAGAGGAGGGGATTAGATGGGGAAATATATTTTCAAACGATTTATTTATATGCTTATTTCTTTATTTATTATTATTACAATTACATTTTTCTTAATGAAATTAATGCCAGGTTCGCCATTTAACGATGCTAAATTAAATGCTGAACAAAAAGAAATTTTAAATGAAAAATATGGATTAAATGATCCTGTAGCTACGCAGTATTTACATTATTTAAAAAATGTTGTTACAGGCG‑ATTTTGGTAA  >  CP000730/988044‑988322

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