Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 362,252 T→A 8.4% T310S (ACC→TCC)  lacY ← lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913362,2520TA 8.4% 103.4 / 2.4 48T310S (ACC→TCC) lacYlactose permease
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base A (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

CTTAAAGAAGCAGAAACAGACCAGATAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAAT  >  NC_000913/362106‑362400
                                                                                                                                                  |                                                                                                                                                    
cTTAAAGAAGCAGAAACAGCCCAGATTAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTGCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTgg                                                                                                                                                    >  1:838190/1‑149 (MQ=255)
 ttAAAGAAGCAGAAACAGACCAGATAAATCGTCGCTGAAAAACGCACTTCAAACTGGCAGGTAATATATTTAAAGCAGCCCACCAGCAGGACCGGTACTTCAAACATATGCAGCGTTTTCAGCATTACCACTTCCAGCGCTGAGGTGGc                                                                                                                                                   >  2:922604/1‑149 (MQ=255)
   aaaGAAGCAGAAACAGACCAGATAAATCGTCGCTGAAAAACTCACTTCAAACTGGCTGGTAATATATTTAAAGCAGACCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCCGCGCTGAGGTGGcta                                                                                                                                                 >  2:172366/1‑147 (MQ=255)
     tgaagCAGAAACAGTCCAGATAAATCGTCGCTGAAAAACGCACTTCAAACTGGATGGTAATATATTTAAAGCAGGCCTCCAGCAGGAACGGTACTTCTAAGATATGCAGAGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAAc                                                                                                                                               <  2:143139/148‑1 (MQ=255)
           aGAAACAGACCAGATAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAg                                                                                                                                         <  2:959907/149‑1 (MQ=255)
                       taTAAATCGTTGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCATAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGt                                                                                                                             <  1:288266/148‑1 (MQ=255)
                        aTAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATTCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAAAACGTa                                                                                                                            >  2:306151/1‑149 (MQ=255)
                         tAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTAc                                                                                                                           <  2:118032/149‑1 (MQ=255)
                         tAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTAc                                                                                                                           <  1:799155/149‑1 (MQ=255)
                            aTCGTCGCTGAAAAACGCACTTCACACTGGCTGGTAAGATATTTAAAGCAGCCCACCAGCAGGACCGGTACTGCAAACATATCCAGCGTTTTCAGAAGACCCCCTGCCCGCGCTGAGGTGGCGAACGAGGAGCCAATACTACGGAccgc                                                                                                                        >  2:911856/1‑146 (MQ=255)
                              cgtcgCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGGAGGAACGGTACTTCAAACATATGAAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTAcacaca                                                                                                                      <  2:266474/149‑1 (MQ=255)
                                  gCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCTGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAAGAATACGTACAGACataat                                                                                                                  >  1:605573/1‑149 (MQ=255)
                                     gAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACATCAAACATATGCCGCGTTTTCAGAATAACCACTACCAGCGCTGATGTCGCGAACGATGAGCCAATAATACGTACAGACATAATAGt                                                                                                               >  1:753050/1‑149 (MQ=255)
                                        aaaCGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCATCAGCAGGAAAGGTACTTCAAACATTTGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTgcc                                                                                                            <  1:734853/149‑1 (MQ=255)
                                            gCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTgccagcc                                                                                                        >  1:986803/1‑149 (MQ=255)
                                             cACTTCAAACTGGCTGGTCATATATTTAAAGACGACCCCCATCATTAACGGTACTTCTAACATATTGAGCGTTTTCTGAAAAACCACTTCCGGCGGTTAGGTGGGGAACGATGAGCCCTTAATACGTACAGACATAATAGTgcgaggca                                                                                                       <  2:516906/149‑1 (MQ=255)
                                                ttCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGCCCGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCcagca                                                                                                    >  1:736492/1‑149 (MQ=255)
                                                  cAAACTGGCTGGTAATATATTTAAAGCAGCCCCCCAGCAGGAACGGTACTTCAAACATATGCCGCGTTTTCAGAATACCCACGTCCCGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTTCAGCCATAATAGTGCCAGCcagcagc                                                                                                  >  1:606299/1‑149 (MQ=255)
                                                       tgtctggTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCCGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGc                                                                                             <  1:907082/146‑1 (MQ=255)
                                                          gtggTAATATATTTCTAGCAGCCATCCAGCAGGTACGGTACTTCATAAATATGCATCGTTTTCAGAATAACCCCTTCCAGCGCTGAGGGGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGtt                                                                                          <  1:575164/148‑1 (MQ=255)
                                                          ctggTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGATGAGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGtt                                                                                          >  2:448267/1‑149 (MQ=255)
                                                          ctggTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGAGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGtt                                                                                          >  2:452093/1‑149 (MQ=255)
                                                              tAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTc                                                                                      >  2:616682/1‑149 (MQ=255)
                                                                  atatTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCAc                                                                                  >  2:352249/1‑149 (MQ=255)
                                                                    atTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCg                                                                                >  2:342354/1‑149 (MQ=255)
                                                                             cgaccactaGCAGGAACGGTACTTCTAACATATGCAGCTTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTa                                                                       <  2:508037/146‑1 (MQ=255)
                                                                             agcCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTa                                                                       <  1:106980/149‑1 (MQ=255)
                                                                                 caccaGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGa                                                                   >  1:152948/1‑149 (MQ=255)
                                                                                       cagGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGAc                                                                                                                       <  2:770510/91‑1 (MQ=255)
                                                                                       cagGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGAc                                                                                                                       >  1:770510/1‑91 (MQ=255)
                                                                                             cGGTACTTCAAACATATGCAGCGTTTTCAGAATATCCACTTCCAGCGCTGAGGGGGCGAACGATGAGCCAATCATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCaa                                                       <  2:87265/149‑1 (MQ=255)
                                                                                               gTACTTCAAACATATGCATCGTTTTCAGAATAACCTCTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCTGCCAGCAGCAGGGCGTTTTTCCCTCCGATGCGATTAATGATCAGTGGCGCAAAg                                                     <  2:753050/149‑1 (MQ=255)
                                                                                                    tCAAACATATGCAGGGTTTTCAGAATAACCACTTCCAGCGCTGCGGTGGCGAACGATGAGCCAATAAAACGTACATCCATAATAGTGCCATCCAGCAGGAGGTCG‑TTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATa                                               <  2:632451/149‑1 (MQ=255)
                                                                                                    tCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACAt                                                >  1:260008/1‑149 (MQ=255)
                                                                                                     cAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCAAGCCAGCAGCAGGGCGTTTTTCCCACAGATGCGATTAATGATCAGTGGCGCAAAGAACATa                                               >  2:555189/1‑149 (MQ=255)
                                                                                                        aCATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATc                                            >  2:915699/1‑149 (MQ=255)
                                                                                                        aCAAATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGGTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATc                                            <  1:483349/149‑1 (MQ=255)
                                                                                                          atatGCAGCTTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCTGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGa                                          <  2:827226/149‑1 (MQ=255)
                                                                                                           tatGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGAGTAATGATCAGTGGCGCAAAGAACATAATCGAg                                         >  2:534527/1‑149 (MQ=255)
                                                                                                                  cGTTATCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGAGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTaa                                  <  1:485412/149‑1 (MQ=255)
                                                                                                                            aaaaCCACTTCCAGCGCTGAGGTGGCGAACGAAGAGCCAATAATACGTTCAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGcc                        <  2:605573/147‑1 (MQ=255)
                                                                                                                            aTAACCACTTCCAGCGCTGAGGAGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGcc                        <  1:452093/149‑1 (MQ=255)
                                                                                                                            aTAACCACTTCCAGCGCTGAGGAGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGcc                        <  1:448267/149‑1 (MQ=255)
                                                                                                                               aCCACTTCCAGCGCTGAGGTGGATAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCAt                     <  1:306151/149‑1 (MQ=255)
                                                                                                                                 cACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTCTTTCCCACCGAGGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCTCCCATTg                   >  2:634406/1‑149 (MQ=255)
                                                                                                                                  acatCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGAGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGt                  <  1:922604/146‑1 (MQ=255)
                                                                                                                                        aGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTAc            >  2:199589/1‑149 (MQ=255)
                                                                                                                                               aGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCa     >  2:43172/1‑149 (MQ=255)
                                                                                                                                                  tGGCGAACGATGAGCCAATAATACGTACAGATATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGACTTAAGTAATTCGCCCAGTGTCGTTACGTAGCCaaaa  >  2:905369/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTAAAGAAGCAGAAACAGACCAGATAAATCGTCGCTGAAAAACGCACTTCAAACTGGCTGGTAATATATTTAAAGCAGCCCACCAGCAGGAACGGTACTTCAAACATATGCAGCGTTTTCAGAATAACCACTTCCAGCGCTGAGGTGGCGAACGATGAGCCAATAATACGTACAGACATAATAGTGCCAGCCAGCAGCAGGGCGTTTTTCCCACCGATGCGATTAATGATCAGTGGCGCAAAGAACATAATCGAGGCGTTAAGTAATTCGCCCATTGTCGTTACGTAGCCAAAT  >  NC_000913/362106‑362400

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