Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,010,363 C→T 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,3630CT97.1% 107.9 / NA 35C146Y (TGT→TAT) ygaDABC superfamily ATP binding cassette transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/19);  minor base G (1/0);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTG  >  CP000730/2010229‑2010487
                                                                                                                                      |                                                                                                                            
cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGACGATAAACAGTGCTGCTAAATTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATGTa                                                                                                                        >  1:47739/1‑141 (MQ=255)
cACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTa                                                                                                                        >  1:52340/1‑141 (MQ=255)
 aCGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTaa                                                                                                                       <  1:16839/141‑1 (MQ=255)
  cgttttatatttCTTAATCTTCTTAAGAACACGTACTCCGTTAAAATGTAAAATGGGAAGATATACATTGCTACTAAAGGCAGTTTCACATCTAAAAAAAACAGTATGGAAAGTGCAATAATAATAGTTATATAATCTAAc                                                                                                                      <  1:121883/132‑1 (MQ=255)
           tttCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTc                                                                                                             <  2:191826/141‑1 (MQ=255)
              ctaactCTTCCAAAGAAAACGTACACCGTTAAAATTGAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTGCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAACTAACCAAATATTCAtt                                                                                                          >  2:122252/6‑141 (MQ=255)
               ttAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAAGAGAACCAAATATTCATTa                                                                                                         >  1:2589/1‑141 (MQ=255)
                 aaTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAc                                                                                                       >  1:187061/1‑141 (MQ=255)
                  aTCTTCCAAAGCAAACGTACACCGTTAAAATGTAAAAGGGGAAGATAAACAGTGCTGATAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAcc                                                                                                      >  2:172266/1‑141 (MQ=255)
                  aTCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAAcc                                                                                                      <  1:140596/141‑1 (MQ=255)
                       ccAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGtt                                                                                                 >  2:140211/1‑141 (MQ=255)
                       ccAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGtt                                                                                                 <  1:639/141‑1 (MQ=255)
                             aaaaCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATg                                                                                           <  2:83041/141‑1 (MQ=255)
                               aaCGTACACCGTTAAAATGTAAAATGGGCAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATACTTGTTATATACTCTAACCAAATATTCATTAACCCGGTTAAAATGaa                                                                                         <  2:111885/141‑1 (MQ=255)
                                aCGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGaaa                                                                                        <  2:31580/141‑1 (MQ=255)
                                 cGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAAt                                                                                       >  1:243501/1‑141 (MQ=255)
                                  gTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATc                                                                                      >  1:123219/1‑141 (MQ=255)
                                  gTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATc                                                                                      >  1:85138/1‑141 (MQ=255)
                                      acCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCtttt                                                                                  <  2:103428/141‑1 (MQ=255)
                                                          gAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAAt                                                              <  1:86687/141‑1 (MQ=255)
                                                            aGATAAACAGTGATGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCa                                                            <  2:47739/141‑1 (MQ=255)
                                                             gATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCAc                                                           >  1:171228/1‑141 (MQ=255)
                                                               tAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCActc                                                         <  2:187061/141‑1 (MQ=255)
                                                                   cAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGa                                                     >  2:111975/1‑141 (MQ=255)
                                                                     gTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata                                                   <  1:229338/141‑1 (MQ=255)
                                                                     gTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGata                                                   <  1:140211/141‑1 (MQ=255)
                                                                                   aaTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACt                                     <  2:107689/141‑1 (MQ=255)
                                                                                    aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATCTAATCTAACCAAATATTCTTTTACCTGGTTAAAATGATATCTTGTGTTTGTTCATCATCATTAATCACTCTAGATATTACTTGACCTACtt                                    >  1:212042/1‑141 (MQ=255)
                                                                                    aTTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACtt                                    >  2:215236/1‑141 (MQ=255)
                                                                                          cgtgtAAAAAGAACATTATGGAAAGTGCAATAATAATTGTTATATAATCTAACCAGATATTCATTAACCCGCTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCGCTCTCGATATTACTTGACCTACTTGattat                              <  2:52340/137‑1 (MQ=255)
                                                                                                  aaGAACATTATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATaa                      >  2:167725/1‑141 (MQ=255)
                                                                                                          tATGGATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGc              <  2:286499/141‑1 (MQ=255)
                                                                                                              gATAGTGCAATAATAATTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACtt          <  2:140657/141‑1 (MQ=255)
                                                                                                                  gTGCAATAATAATTGTGATAGAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAg      >  2:137270/1‑141 (MQ=255)
                                                                                                                      aataataatTGTTATATAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTg  <  2:171228/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                            
CACGTGTCAATTTTCTTAATCTTCCAAAGAAAACGTACACCGTTAAAATGTAAAATGGGAAGATAAACAGTGCTGCTAAAGTCAATTTCACATCTAAAAAGAACATTATGGATAGTGCAATAATAATTGTTATACAATCTAACCAAATATTCATTAACCCGGTTAAAATGAAATCTTTTGTTTGTTCAACATCATTAATCACTCTAGATATTACTTGACCTACTTGATTATTAGCATAAAATCTCGCACTTAAAGCTTG  >  CP000730/2010229‑2010487

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: ‑