Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,075,342 G→A 8.7% L13L (CTG→TTG)  ygaY ← transporter YgaY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,075,3420GA 8.7% 127.9 / 2.3 46L13L (CTG→TTG) ygaYtransporter YgaY
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (2/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

CAGGCCGGCGCCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGATCGATGAAC  >  NC_002947/2075196‑2075490
                                                                                                                                                  |                                                                                                                                                    
cAGGCCGGCGCCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATcagca                                                                                                                                                    <  5:189832/149‑1 (MQ=255)
   gCCGGCGCCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGa                                                                                                                                                 >  8:259886/1‑149 (MQ=255)
        cgcCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAg                                                                                                                                            <  1:34824/149‑1 (MQ=255)
         gcCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAgg                                                                                                                                           >  3:175738/1‑149 (MQ=255)
            aTAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGgcgc                                                                                                                                        >  8:265597/1‑149 (MQ=255)
                 tGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGGGGCCATCAGCAGGATCAAGGCGCGGCt                                                                                                                                    >  3:109027/1‑148 (MQ=255)
                   aGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTga                                                                                                                                  <  4:175738/148‑1 (MQ=255)
                          cgGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGATAGTgggg                                                                                                                          >  4:98894/1‑149 (MQ=255)
                              aaTCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAg                                                                                                                                            <  1:269425/127‑1 (MQ=255)
                              aaTCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAg                                                                                                                                            >  2:269425/1‑127 (MQ=255)
                                         tgcCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATtc                                                                                                           >  2:101378/1‑149 (MQ=255)
                                         tgcCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATtc                                                                                                           >  2:169098/1‑149 (MQ=255)
                                           ccGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATtctc                                                                                                         >  7:173063/1‑149 (MQ=255)
                                                 cTGGCGGTGCTCAAGCCAAACTTCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTgggg                                                                                                   <  1:144380/149‑1 (MQ=255)
                                                 cTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCATGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTgggg                                                                                                   <  5:273891/149‑1 (MQ=255)
                                                  tGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTggggg                                                                                                  >  8:343188/1‑149 (MQ=255)
                                                    gCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGgtgt                                                                                                                      >  6:240375/1‑127 (MQ=255)
                                                      ggTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAACAGGATCAAgg                                                                                                                                           >  5:212463/1‑104 (MQ=255)
                                                      ggTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAACAGGATCAAgg                                                                                                                                           >  5:207082/1‑104 (MQ=255)
                                                      ggTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAACAGGATCAAgg                                                                                                                                           <  6:207082/104‑1 (MQ=255)
                                                      ggTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAACAGGATCAAgg                                                                                                                                           <  6:212463/104‑1 (MQ=255)
                                                           tCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGgtgt                                                                                         >  6:275514/1‑149 (MQ=255)
                                                               gagataCTGCTGCGCGATGCTTTGCAGCTGTGGTTGGGCGTATTTGTTGCTGGCTTGGGCCAGGGCTGTGGCGGTGTCCATGAGCAGGATAAAGGCGCGGCTGAGAGTGGGGGTGTTCATGCAATTCTCGTGGCGGCTGTGGTGCGTAg                                                                                     <  8:339276/144‑1 (MQ=255)
                                                               gCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAg                                                                                     <  7:259886/149‑1 (MQ=255)
                                                                ccAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGt                                                                                    >  8:114888/1‑149 (MQ=255)
                                                                     ctgctgGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGtgatga                                                                               <  4:77037/149‑1 (MQ=255)
                                                                                  cTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAgt                                                                   <  6:166464/148‑1 (MQ=255)
                                                                                    gtgCAGCAGTGGCTCCGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGGGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTc                                                                <  4:256021/149‑1 (MQ=255)
                                                                                          cagTGGTTGGGCGTAATCGCTCCGGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGCCGCGGCTGAGAGGGGGGGGGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCAt                                                          <  8:173063/149‑1 (MQ=255)
                                                                                              ggTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTTATAGTGTGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGac                                                       <  5:24866/148‑1 (MQ=255)
                                                                                                  gggCGTAATAGTTGCTGGCCACGGCCAGGCCTCTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCaa                                                  >  1:245242/1‑149 (MQ=255)
                                                                                                   ggCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCATGCCTGacactta                                                  <  1:169098/148‑4 (MQ=255)
                                                                                                      gTAATAGTTGCTGGCCACGGCCAGGCCTGGGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGta                                              <  4:265613/149‑1 (MQ=255)
                                                                                                       tAATCGTTGCTGGCCCCGGCCATGCCTGTGGCGGTGGCCATCAGCAGGATCCAGGCGCGGCTGCGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGtat                                             <  5:122263/149‑1 (MQ=255)
                                                                                                       tAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCGGAGAGTGGGGGGGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGtat                                             <  4:325385/149‑1 (MQ=255)
                                                                                                               gCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGtttt                                       >  3:194646/1‑147 (MQ=255)
                                                                                                                    ccccGGCCAGGCCTGTCGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGACTGGGGGTGTGCATGCAATTCTCCTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTa                                <  4:355276/146‑1 (MQ=255)
                                                                                                                         gCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCGGAGAGTGGGCGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCAt                           <  7:265597/149‑1 (MQ=255)
                                                                                                                         gCCAGGCATGTGGCGGTGGCCATCAGCAGGATAAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCa                            >  1:259305/1‑148 (MQ=255)
                                                                                                                          ccAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGGGGGGGTGTGCATGCATTTGTCGTGGGGGCAGTGGTGGGTTGTGATGCTTATCAGTAAGTGTCGGGCAGGCCACCAATGTATAGTTTTCAGCCTATTCATc                          >  2:211814/1‑149 (MQ=255)
                                                                                                                           cAGGCCTGTGGCTGTGGCCATCAGCATGATCACGGCGCCGCTGACAGTGGGCGTGTCCATGCAATTCTCCTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCt                         <  5:275514/149‑1 (MQ=255)
                                                                                                                           cAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCt                         <  2:245242/149‑1 (MQ=255)
                                                                                                                                cggtcgCGGTCGCCATCAGCAGGATCAAGCCGCCGCTGAGAGGGGGCGCGTGCAGCCAATTCTCGTGGGCGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGaaa                    <  7:114888/144‑1 (MQ=255)
                                                                                                                                 tgtgGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAg                   >  5:2403/1‑149 (MQ=255)
                                                                                                                                             cATCAGCCGGATCCCCGCGCGGCGGAGAGTGGGCGTCGCCAAGCAACTCTCGGCGGGGCAGTGGTGTGAAGTGATGATTATCAGTAAGGGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGAGCGa       <  4:109027/149‑1 (MQ=255)
                                                                                                                                             cATCAGCAGGCTCAAGGCGCGGCTGAGAGGTGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGATCGa       <  3:98894/149‑1 (MQ=255)
                                                                                                                                             cATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGATCGa       >  3:381877/1‑149 (MQ=255)
                                                                                                                                                  gcagGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTGTCGTGGGGGCTGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGCCACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGATCGATGAAc  >  5:108361/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGGCCGGCGCCATAGCTGAGTTGCGCGGCAATCACGATGCTGCCGGCGCTGGCGGTGCTCAAGCCAAACTGCTGGGCGATGCTGTGCAGCAGTGGTTGGGCGTAATAGTTGCTGGCCACGGCCAGGCCTGTGGCGGTGGCCATCAGCAGGATCAAGGCGCGGCTGAGAGTGGGGGTGTGCATGCAATTCTCGTGGGGGCAGTGGTGTGTAGTGATGATTATCAGTAAGTGTCAGGCATGACACCAATGTATAGTTTTCAGCCTATTCATCTTGAAAGCAGATAGATCGATGAAC  >  NC_002947/2075196‑2075490

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