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% 187.3 / ‑4.4 61C146Y (TGT→TAT) ygaDABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/22);  minor base A (0/1);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-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.

CACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATG  >  CP000730/2010229‑2010493
                                                                                                                                      |                                                                                                                                  
cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTa                                                                                                                              >  2:21528/1‑141 (MQ=255)
cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTa                                                                                                                              >  2:636663/1‑141 (MQ=255)
         aTTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAatat                                                                                                                     >  1:218111/1‑141 (MQ=255)
             tCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAAAATTATGGATAGTGCAATAATAATTGTTATCTAATCTAACCAAATATTCAt                                                                                                                 >  2:246054/1‑141 (MQ=255)
             tCTTAAGCTTCCAAAGAAAACGCACACAGGTAAAATATAAAATGGGATGATAAACAGTGTTGCTAAAGTCAATCTCACATCTAAAAGGAACATGATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCAt                                                                                                                 <  2:185149/141‑1 (MQ=255)
                tAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTaa                                                                                                              <  1:332819/141‑1 (MQ=255)
                  aTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAcc                                                                                                            >  1:387412/1‑141 (MQ=255)
                   tCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAccc                                                                                                           <  1:765188/141‑1 (MQ=255)
                      tCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAAACAAATATTCAATAACCCGGt                                                                                                        >  2:327743/1‑141 (MQ=255)
                         aaagaaagCGTACACCGTTAGAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                                     >  1:511289/1‑141 (MQ=255)
                         aaaGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTaa                                                                                                     >  1:268694/1‑141 (MQ=255)
                              aaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACACCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGa                                                                                                >  1:634791/1‑141 (MQ=255)
                              aaaCGTACACCGTTAAAATGTAAAATGGGAAGAAAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGACAGTGCAATAATAAATGTTAAATACTCAGTTCATATTTTCGTTAACCCGGTTAAAAtaa                                                                                                >  2:329241/1‑139 (MQ=255)
                                  gTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAtt                                                                                            >  1:373217/1‑140 (MQ=255)
                                  gTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCCGCTAAAGTCATTTACACATCTATACAGAACATTATGGTTGGTGCAATAATAATTGTTATATTATCTAATGATATTTTCTTTGACCCGGTTAAAATGAAATc                                                                                            >  1:318982/1‑141 (MQ=255)
                                   tACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCt                                                                                           >  1:670014/1‑141 (MQ=255)
                                   tACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCt                                                                                           >  2:537921/1‑141 (MQ=255)
                                      acCGTTAATATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                                                                        <  2:686094/141‑1 (MQ=255)
                                      acCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATGTTCGTTAACCCGGTTAAAATAAAATAtttt                                                                                        >  2:423833/1‑141 (MQ=255)
                                         gTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgtt                                                                                     >  2:461755/1‑141 (MQ=255)
                                         gTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATATtttgtt                                                                                     >  2:684867/1‑141 (MQ=255)
                                         gTTAAAATGTAAAATGAGACGATTAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTtttgtt                                                                                     >  2:574802/1‑141 (MQ=255)
                                                    aaatgGGTAGGTAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatc                                                                          <  2:584107/141‑1 (MQ=255)
                                                    aaatgGGAAGATTAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatc                                                                          <  2:808156/141‑1 (MQ=255)
                                                      atgGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAcatcat                                                                        <  2:634791/141‑1 (MQ=255)
                                                        gggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTa                                                                      <  1:633157/141‑1 (MQ=255)
                                                         ggAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTaa                                                                     >  1:331861/1‑141 (MQ=255)
                                                               tAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTATTCActc                                                               >  1:482414/1‑141 (MQ=255)
                                                                aaaCAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActct                                                              <  2:619993/141‑1 (MQ=255)
                                                                     gttcggcgaaaGTCAATTTCACATCTAAAAAGCACATTATGGATAGTGCAATAATAATTGGTATATAATCTAACCAAATATTCATTAACACGGTTAAAATGAAATCTTTTGTTTGTTGACAATAATTATTACCTCTAGatt                                                         >  2:8208/9‑140 (MQ=255)
                                                                     gTGCTGCTAAAGTCCATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata                                                         >  1:442851/1‑141 (MQ=255)
                                                                     gTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCACCATCATTAATCACTCTAGatt                                                         >  2:65133/1‑140 (MQ=255)
                                                                     gTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata                                                         >  2:65107/1‑141 (MQ=255)
                                                                      tgctgcTAAAGTCAATTTTACATCTAAAAAGAACATTATGGATAGTGCAAAAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGatat                                                        <  1:107300/141‑1 (MQ=255)
                                                                        ctgctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAAAATCATTAATAAATCTAAATAttc                                                      >  2:180489/1‑140 (MQ=255)
                                                                          gctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACt                                                    <  1:202367/141‑1 (MQ=255)
                                                                          gctAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACt                                                    >  1:234250/1‑141 (MQ=255)
                                                                                 tCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTa                                             <  2:522543/141‑1 (MQ=255)
                                                                                 gcAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTa                                             <  1:10459/140‑1 (MQ=255)
                                                                                  cAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACTTAc                                            >  2:363814/1‑141 (MQ=255)
                                                                                  cAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCtca                                            >  2:198479/1‑139 (MQ=255)
                                                                                   aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                           <  2:511289/141‑1 (MQ=255)
                                                                                   aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                           <  2:268694/141‑1 (MQ=255)
                                                                                    aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACtt                                          <  1:240883/141‑1 (MQ=255)
                                                                                       tCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAACTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTATTTGACCTACTTGat                                       <  1:423833/141‑1 (MQ=255)
                                                                                         acaTCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGatta                                     >  1:510756/1‑141 (MQ=255)
                                                                                              tAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGc                                <  2:755415/141‑1 (MQ=255)
                                                                                              tAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGc                                <  2:571464/141‑1 (MQ=255)
                                                                                                  aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                            <  2:387412/141‑1 (MQ=255)
                                                                                                  aaGAACAGTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTCCAACCTAATTAGAGAAGCTAGGTAGTACTTGACCTACTTGATTATTAGCATaa                            >  2:567747/1‑141 (MQ=255)
                                                                                                    gAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTCTTAGCATaaaa                          >  2:429622/1‑141 (MQ=255)
                                                                                                     aaCATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAt                         >  1:631817/1‑141 (MQ=255)
                                                                                                        aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAAtctc                      >  2:341924/1‑141 (MQ=255)
                                                                                                        aTTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCACTAATCACTCTAGATATTACTTCACCTACTAGATTATTAGAATAAAAcctt                      <  2:594847/141‑2 (MQ=255)
                                                                                                         ttATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCg                     >  1:38324/1‑141 (MQ=255)
                                                                                                              aaaaGTACAATAATAATTTTTATATAATCTAACCAAATATTCATTAACCCTGTTAAAATGAAATCTTTTGTTCAATCAACAACATAAAACACTCTAAATATTACTTGACCTACTTGATTACTAGCATAAAATCTCGCACtt                <  2:530086/138‑1 (MQ=255)
                                                                                                                    gCAATAATAATTGTTTTATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCt          >  1:601539/1‑141 (MQ=255)
                                                                                                                    gCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACCTAAAGCt          >  1:291760/1‑141 (MQ=255)
                                                                                                                      aataataatTGTTATAATATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTg        <  2:321298/141‑1 (MQ=255)
                                                                                                                            aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTCCTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATg  >  1:159721/1‑141 (MQ=255)
                                                                                                                            aatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACATACTTGATTATTAGCATAAAATCTCGCACTTGAAGCTTGTAAATg  >  1:433225/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTGTAAATG  >  CP000730/2010229‑2010493

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