Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,804,319 T→C 12.1% intergenic (+14/‑149) yniC → / → ydjM 2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,804,3190TC12.1% 75.1 / 3.2 33intergenic (+14/‑149)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

GCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTAT  >  CP009273/1804186‑1804466
                                                                                                                                     |                                                                                                                                                   
gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAACATAAACCCCGTTGTGTATTAGCAGACGCCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCGCGGTCAGTGACCAGGGACAGTGCCCTCGCt                                                                                                                                            >  2:300656/1‑143 (MQ=255)
gcCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                                                                      >  1:510735/1‑149 (MQ=255)
  ccGCATGCGTTCCATCGTCGTTCCTGCGCCACAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                                                                    <  1:434466/149‑1 (MQ=255)
   cGCATGCCTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGt                                                                                                                                   <  1:283752/149‑1 (MQ=255)
               aTCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCCGCCATCTCGCTGCCCCCGGTTCTTTa                                                                                                                             >  2:208116/1‑143 (MQ=255)
                      ttCCTCCGCCAGAAGCGCAAAATGATCCACGTCCTGGATTAGCAGGCGTCAAACTTTCATCGCTTACAAAACTAACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGtcattaaccgagaa                                                                                                                      <  2:306352/143‑14 (MQ=21)
                          tGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                                                                  >  2:477975/1‑143 (MQ=255)
                             gcCAGAAGCGCAAAATGATCGACGTTTTGTATTAGCAGACGTCGAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGACATTGCAtt                                                                                                               >  2:39714/1‑143 (MQ=255)
                                     cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCCGGTTCTTTATCTGAAATTGCATTCAACTGat                                                                                                       <  2:346904/143‑2 (MQ=255)
                                              aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                                                                              >  2:180528/1‑143 (MQ=255)
                                               tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTGATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTGAACTGACGGATTAATCg                                                                                             >  2:294502/1‑143 (MQ=255)
                                                      tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAgaga                                                                                >  1:352212/1‑149 (MQ=255)
                                                      tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTAGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCTTCAATTTAAgaga                                                                                >  1:352176/1‑149 (MQ=255)
                                                       ttGGATTTGCAGACGTTAAACTTGCATGGCCGACAGAACTCCCCGCAAAAGACCTCCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTTTTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                                                                     <  2:190359/143‑1 (MQ=255)
                                                          tATTAGCAGACGTCAAACTTTCATCACTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                                                            >  1:264822/1‑149 (MQ=255)
                                                             tAGCAGTCGTCAAACTTTCCTCGCTGACAGAACTCACCGCAAATGACCTTCTCGGTTAATCACCAGGGGCAGTGACCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt                                                                         <  1:178674/149‑1 (MQ=255)
                                                              cgCAGACGTCAAACTTTCACCGCCGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa                                                                              <  2:384408/142‑1 (MQ=255)
                                                               gcggACGGCAAACTTTCATCGGTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGGCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCTATTTAAGAGAAag                                                                             <  2:288053/140‑1 (MQ=255)
                                                                  gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGCTTAATCGTCAATTTAAGAGAAAGAGTTACAcc                                                                    >  1:436866/1‑149 (MQ=255)
                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                                                                    >  2:309962/1‑143 (MQ=255)
                                                                        aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcccc                                                              >  1:438594/1‑149 (MQ=255)
                                                                          aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGt                                                                  <  2:459364/143‑1 (MQ=255)
                                                                                 tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt                                                           >  2:269136/1‑143 (MQ=255)
                                                                                    cTGACAGAACTCACCGCAAAAGACCTCCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAc                                                  <  1:510616/149‑1 (MQ=255)
                                                                                      gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGt                                                      >  2:200984/1‑143 (MQ=255)
                                                                                      gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACAGGGGGCAGTGATATCTTGGCCCCTGGTTCCTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTg                                                <  1:101149/149‑1 (MQ=255)
                                                                                        cAGAACTCCCCGCAACAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCCATTCACAATAAAGAGTTACACCGTCACCACTTCCGTGCACTGta                                              <  1:258398/149‑1 (MQ=255)
                                                                                                     aaaaGACCTTCCCGGCTATTTACCAGTGGCAGTGATCTCGCTGCCCCCGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCtata                                 <  1:200405/149‑1 (MQ=255)
                                                                                                          aCCTTCGCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTa                            <  1:76416/149‑1 (MQ=255)
                                                                                                                 agataaagaaCCAGGGACAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAGTTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTCTAAAAACCCTATACTGtct                           >  2:455909/10‑141 (MQ=37)
                                                                                                                   gTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGAc                   >  1:414288/1‑149 (MQ=255)
                                                                                                                                ggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTACCGACAGTTTAg            >  2:162444/1‑143 (MQ=255)
                                                                                                                                    gTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTAt  <  1:208143/149‑1 (MQ=255)
                                                                                                                                     tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACGGTCCGTATCGACAGTTTAGTGAGt       >  2:321126/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
GCCCGCATGCGTTCCATCGTCGTTCCTGCGCCAGAAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTATAAAAATCCTATACTGTACGTATCGACAGTTTAGTGAGTTTTAT  >  CP009273/1804186‑1804466

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: