Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,853 T→C 16.5% N129D (AAT→GAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,8530TC16.5% 179.5 / 26.2 78N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/38);  new base C (3/10);  total (30/48)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATAT  >  NC_000913/1197754‑1197950
                                                                                                   |                                                                                                 
tCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGAtt                                                                                                   <  1:997529/100‑1 (MQ=255)
  gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGc                                                                                                 <  1:29586/100‑1 (MQ=255)
    aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTc                                                                                               <  1:371969/100‑1 (MQ=255)
      aTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCAt                                                                                             <  2:1314004/100‑1 (MQ=255)
         cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTAc                                                                                          <  2:1606056/100‑1 (MQ=255)
         cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTAc                                                                                          <  1:1667504/100‑1 (MQ=255)
           aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAt                                                                                        <  1:631188/100‑1 (MQ=255)
             aaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATtg                                                                                      <  1:941650/100‑1 (MQ=255)
               atatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATtgtg                                                                                    <  2:136423/100‑1 (MQ=255)
               aaaaaaaaaaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATtgtg                                                                                    <  2:1127823/94‑1 (MQ=255)
                tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGa                                                                                   <  1:1067066/100‑1 (MQ=255)
                 atCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAg                                                                                  <  2:781605/100‑1 (MQ=255)
                 atCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAg                                                                                  <  1:1021881/100‑1 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCtt                                                                               >  1:1615095/1‑100 (MQ=255)
                        aaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATc                                                                           <  1:1373337/100‑1 (MQ=255)
                            gTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCttta                                                                       <  2:994169/100‑1 (MQ=255)
                             tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttat                                                                      <  2:1447070/100‑1 (MQ=255)
                              ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                     <  1:1191735/100‑1 (MQ=255)
                                 cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttattta                                                                   >  1:157170/1‑99 (MQ=255)
                                    tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                               <  1:2054149/100‑1 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                           >  2:1254138/1‑100 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                           <  1:149937/100‑1 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                          <  2:2037523/100‑1 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCGTTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                          >  1:1828727/1‑100 (MQ=255)
                                          agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                         <  2:1808061/100‑1 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                        >  1:369342/1‑100 (MQ=255)
                                             agagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGAc                                                      >  1:796640/1‑100 (MQ=255)
                                             agagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGAc                                                      <  1:63439/100‑1 (MQ=255)
                                             agagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCTCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                        <  2:2460237/98‑1 (MQ=255)
                                               aattCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTc                                                    <  1:342810/98‑1 (MQ=255)
                                                gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                   <  1:2238817/100‑1 (MQ=255)
                                                gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                   >  1:1997194/1‑100 (MQ=255)
                                                gTTCAATATAATGTTCGTCTTCAATTTTTCTTACTTCAGAGTAATATAGATTGCTAATTACATTGTGAGCTTCATATTTATTTAATTTTCTTTTGACTcc                                                   >  2:1809320/1‑100 (MQ=255)
                                                  tCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                                 >  1:81297/1‑100 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                                >  1:1089404/1‑100 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                                >  1:730241/1‑100 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                                <  1:1283877/100‑1 (MQ=255)
                                                    aaTATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGct                                               >  2:2487034/1‑100 (MQ=255)
                                                     atatAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                                <  2:1574759/98‑1 (MQ=255)
                                                       ataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCtctc                                            <  1:84718/100‑1 (MQ=255)
                                                       ataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAACTTTCTGTTGACTCCAGCtctc                                            >  2:225073/1‑100 (MQ=255)
                                                        taATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTcc                                           >  1:1298957/1‑100 (MQ=255)
                                                         aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                          <  1:1630464/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCCCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                         <  1:1928801/100‑1 (MQ=255)
                                                           tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                        <  2:157168/100‑1 (MQ=255)
                                                            gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATTGATTGCTCATTACATTGTGAGATTCATCTTTATTTAATTTTCTGTTGACTCCAGTTCTCCGt                                         <  2:998999/98‑1 (MQ=255)
                                                               tGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATaa                                    <  1:2256014/100‑1 (MQ=255)
                                                               tGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATaa                                    <  2:1221495/100‑1 (MQ=255)
                                                                  cTTCAATTTTTTTTTTTTTTTGTTAATATAGATTGCTCATTACATTTTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGAt                                      >  2:2420014/1‑95 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTTAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                 <  1:382345/100‑1 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                 >  1:2214132/1‑100 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                 <  1:387902/100‑1 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                 >  1:1876212/1‑100 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                 >  1:294960/1‑100 (MQ=255)
                                                                   ttCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                                >  2:1765206/1‑100 (MQ=255)
                                                                    tCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtt                               >  1:335806/1‑100 (MQ=255)
                                                                      aaTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtttt                             <  2:2435939/100‑1 (MQ=255)
                                                                       aTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTta                            <  1:2582157/100‑1 (MQ=255)
                                                                         ttttCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtata                          >  2:1649547/1‑100 (MQ=255)
                                                                          tttCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATaa                         <  1:2431538/100‑1 (MQ=255)
                                                                          tttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATaa                         <  2:2361335/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtt                       >  1:238290/1‑99 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                      <  1:2157999/100‑1 (MQ=255)
                                                                                 cTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                  >  1:2364050/1‑100 (MQ=255)
                                                                                 cTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                  >  1:1249702/1‑100 (MQ=255)
                                                                                  ttCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                 <  2:2262207/100‑1 (MQ=255)
                                                                                     aGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa              <  2:369338/100‑1 (MQ=255)
                                                                                     aGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa              <  1:32810/100‑1 (MQ=255)
                                                                                     aGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa              <  2:347064/100‑1 (MQ=255)
                                                                                     aGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa              <  1:565385/100‑1 (MQ=255)
                                                                                      gggTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTAt             >  1:372279/1‑100 (MQ=255)
                                                                                      gggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTATATGCTTAt             >  2:628879/1‑100 (MQ=255)
                                                                                       ggTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATc            <  1:2147012/100‑1 (MQ=255)
                                                                                       ggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATc            <  2:1933897/100‑1 (MQ=255)
                                                                                           atatAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa        >  1:1941331/1‑100 (MQ=255)
                                                                                             ataGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg      <  1:48134/100‑1 (MQ=255)
                                                                                               aGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGat    >  1:745813/1‑100 (MQ=255)
                                                                                                gATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata   >  2:24179/1‑100 (MQ=255)
                                                                                                gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata   >  2:2334118/1‑100 (MQ=255)
                                                                                                  ttGCTCATTAGATTGTGAGCTTCATCTTTATTTAATGTTCTGTTGACTGCAGCTCTGGGTGGTAATGGTGTTGTGATGAGATGTTTTTCTCAAAGatat  >  2:2170695/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATAT  >  NC_000913/1197754‑1197950

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