Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,010,363 C→T 100% C146Y (TGT→TAT)  ygaD ← ABC superfamily ATP binding cassette transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,010,3630CT98.4% 204.3 / ‑3.9 61C146Y (TGT→TAT) ygaDABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/33);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-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.

TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTA  >  CP000730/2010227‑2010499
                                                                                                                                        |                                                                                                                                        
ttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTACAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATCGTGCAATAATAATTGTTATATAATc                                                                                                                                      >  1:168710/1‑141 (MQ=255)
ttCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATc                                                                                                                                      >  2:257403/1‑141 (MQ=255)
 tCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCt                                                                                                                                     <  1:749309/141‑1 (MQ=255)
 tCACGTGTCAATATTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGGTATATAATCt                                                                                                                                     >  2:34252/1‑141 (MQ=255)
  cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTa                                                                                                                                    >  1:356835/1‑141 (MQ=255)
   aCGTGTCAATTTTCTTAAAAAAAAAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTaa                                                                                                                                   <  2:416099/141‑1 (MQ=255)
    cGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTACAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAc                                                                                                                                  >  2:697427/1‑141 (MQ=255)
     gtgtCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAcc                                                                                                                                 >  2:761545/1‑141 (MQ=255)
        tCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAGCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCaaa                                                                                                                              >  2:266413/1‑141 (MQ=255)
           aTTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAatat                                                                                                                           >  2:391153/1‑141 (MQ=255)
              ttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCa                                                                                                                        <  2:618450/141‑1 (MQ=255)
               tCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCAt                                                                                                                       <  1:627425/141‑1 (MQ=255)
                   aaTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                                                                   >  1:557414/1‑141 (MQ=255)
                             aGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaaaa                                                                                                         >  2:644614/1‑141 (MQ=255)
                              gAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAAt                                                                                                        <  1:61696/141‑1 (MQ=255)
                               aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                                       <  2:126315/141‑1 (MQ=255)
                                 aaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaa                                                                                                     >  2:65653/1‑141 (MQ=255)
                                  aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                                    <  2:543535/141‑1 (MQ=255)
                                   cGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAt                                                                                                   >  1:441021/1‑141 (MQ=255)
                                    gTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGACATc                                                                                                  >  2:587244/1‑141 (MQ=255)
                                      acacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtt                                                                                                <  2:745775/141‑1 (MQ=255)
                                      acacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtt                                                                                                <  1:463503/141‑1 (MQ=255)
                                       cacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCttt                                                                                               <  2:770289/141‑1 (MQ=255)
                                       cacCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCttt                                                                                               <  2:597565/141‑1 (MQ=255)
                                         ccGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttg                                                                                             <  1:599766/141‑1 (MQ=255)
                                             taaaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgtttg                                                                                         >  2:382221/1‑141 (MQ=255)
                                              aaaatgtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTttgtttgt                                                                                        <  2:389418/141‑1 (MQ=255)
                                                   gtaaaatgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAc                                                                                   >  2:638994/1‑141 (MQ=255)
                                                          gggAAGATAAACATGGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTa                                                                            <  1:45622/141‑1 (MQ=255)
                                                           ggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTaa                                                                           <  2:719354/141‑1 (MQ=255)
                                                               gATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                                       >  1:229319/1‑141 (MQ=255)
                                                                 tAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAGCCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                                                                     <  1:376899/141‑1 (MQ=255)
                                                                  aaaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCAAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                                    <  2:715048/141‑1 (MQ=255)
                                                                     cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa                                                                 >  1:603611/1‑141 (MQ=255)
                                                                          ctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTa                                                            >  1:684581/1‑141 (MQ=255)
                                                                            gctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTTTTTTTTCAACATAATTAATCACTCTAGATATTACt                                                          >  2:520219/1‑141 (MQ=255)
                                                                               aaaGTCAATTTCACATGTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                       <  1:246161/141‑1 (MQ=255)
                                                                               aaaGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGa                                                       <  2:618508/141‑1 (MQ=255)
                                                                                     aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                                 <  1:724973/141‑1 (MQ=255)
                                                                                     aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                                 <  1:724972/141‑1 (MQ=255)
                                                                                      aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACtt                                                <  1:134662/141‑1 (MQ=255)
                                                                                      aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACtt                                                <  1:380338/141‑1 (MQ=255)
                                                                                        ttCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGa                                              >  2:365971/1‑141 (MQ=255)
                                                                                                  aaaaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTTCTTGACCTACTTGATTATTAGCAt                                    <  1:382221/141‑1 (MQ=255)
                                                                                                    aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                                  <  1:177536/141‑1 (MQ=255)
                                                                                                    aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                                  <  2:437776/141‑1 (MQ=255)
                                                                                                    aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                                  <  2:67832/141‑1 (MQ=255)
                                                                                                          aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtctc                            <  2:555434/141‑1 (MQ=255)
                                                                                                           ttATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCg                           >  1:39607/1‑141 (MQ=255)
                                                                                                             aTGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCa                         >  2:164356/1‑141 (MQ=255)
                                                                                                             aTGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCa                         >  1:671178/1‑141 (MQ=255)
                                                                                                              tGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCAc                        <  1:483629/141‑1 (MQ=255)
                                                                                                               ggATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACt                       <  1:518018/141‑1 (MQ=255)
                                                                                                                    gTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAg                  >  2:506090/1‑141 (MQ=255)
                                                                                                                        aataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTg              <  2:709201/141‑1 (MQ=255)
                                                                                                                         ataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGt             >  1:193361/1‑141 (MQ=255)
                                                                                                                              aatTGTTCTATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATg        >  1:622662/1‑141 (MQ=255)
                                                                                                                              aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                                    <  1:515654/81‑1 (MQ=255)
                                                                                                                              aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                                    >  2:515654/1‑81 (MQ=255)
                                                                                                                              aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATg        <  2:168710/141‑1 (MQ=255)
                                                                                                                                ttGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGt      >  1:45519/1‑141 (MQ=255)
                                                                                                                                    tatataATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGta  <  1:482378/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TTCACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATGGTTGTA  >  CP000730/2010227‑2010499

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