Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,706,297 A→C 5.5% intergenic (‑3/‑155) yhjY ← / → tag autotransporter beta‑domain protein/3‑methyl‑adenine DNA glycosylase I, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,706,2970AC 5.5% 2424.1 / 23.5 914intergenic (‑3/‑155)yhjY/tagautotransporter beta‑domain protein/3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base A (183/668);  new base C (3/60);  total (186/728)
Fisher's exact test for biased strand distribution p-value = 5.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 946 total aligned reads displayed.

CACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCATAGT  >  CP009273/3706231‑3706364
                                                                  |                                                                   
 aCTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAac                                                                    >  1:19284702/1‑67 (MQ=255)
 aCTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAac                                                                    >  1:28480952/1‑67 (MQ=255)
    gATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACg                                                                 >  1:3678133/1‑67 (MQ=255)
     aTTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGAc                                                               <  1:1646437/68‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGGTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACt                                                              <  1:3550521/68‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCTACCACCGCTTTTTTTTATGATCATTAACACGACt                                                              <  1:17607250/68‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATTATCATTAACACGACt                                                              <  1:14398390/68‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCTTTAACACGACt                                                              <  1:22270234/68‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGAc                                                               <  1:6917843/67‑1 (MQ=255)
      ttACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACt                                                              <  1:38348960/68‑1 (MQ=255)
      ttACCACGCGCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACt                                                              <  1:18942401/68‑1 (MQ=255)
        aCCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTcc                                                            <  1:9200690/68‑1 (MQ=255)
         ccACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCa                                                           >  1:28176149/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGGGACCACCGCTTTTTTTTGTGATCATTAACACGACTCCAgg                                                         >  1:47390382/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAgg                                                         >  1:11766747/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAgg                                                         >  1:33448118/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAgg                                                         >  1:11197831/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAgg                                                         >  1:20285163/1‑68 (MQ=255)
           aCGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTAATGATCATTAACACGACGCCAgg                                                         >  1:17114710/1‑68 (MQ=255)
             gCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGCCTCCAGGGa                                                       >  1:40420145/1‑68 (MQ=255)
             gCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGa                                                       >  1:30337418/1‑68 (MQ=255)
             gCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGa                                                       >  1:51091764/1‑68 (MQ=255)
               tCGCCAGCTGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAag                                                     <  1:15667257/68‑1 (MQ=255)
               tCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAag                                                     >  1:880962/1‑68 (MQ=255)
                       agGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGa                                               <  1:36100209/66‑1 (MQ=255)
                       agGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAAt                                             >  1:20855176/1‑68 (MQ=255)
                         gCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTACCACGACTCCAGGGAAGAGCCGAATTc                                           >  1:46250166/1‑68 (MQ=255)
                           ttAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTGACACGACTCCAGGGAAGAGCCGAATTCg                                          <  1:31193305/67‑1 (MQ=255)
                            tcgCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGc                                        <  1:2405574/66‑1 (MQ=255)
                            taccTGCCAGGGACCACCGCTTTTTTTTATGATCATAAACACGACTCCAGGGAAGAGCCGAATTCGGc                                        <  1:34599195/65‑1 (MQ=255)
                               cTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     >  1:40775614/1‑68 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGATTTCGGCAata                                    >  1:48962754/1‑68 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAata                                    <  1:19793225/68‑1 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAata                                    >  1:24532969/1‑68 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAata                                    >  1:32235289/1‑68 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     <  1:37707319/67‑1 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     <  1:41596492/67‑1 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     <  1:7824017/67‑1 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGATAATTACCACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     <  1:42002190/67‑1 (MQ=255)
                                tGCCAGCGACCACCGCTTTTTTTTATGACCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAat                                     <  1:43201768/67‑1 (MQ=255)
                                 gCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:34062619/68‑1 (MQ=255)
                                 gCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:49490072/68‑1 (MQ=255)
                                 gCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:27097089/68‑1 (MQ=255)
                                 gCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGAATCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:12716722/68‑1 (MQ=255)
                                 gCCAGCGACCACCGCTTTTTTTTATGATCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:45352809/68‑1 (MQ=255)
                                  ccAGCGACCCCCGCTTTTTTTTATGATCATAACCACGACTCCAGGGAAGAGCCGAATTCGGCAatat                                   <  1:21690564/67‑1 (MQ=255)
                                   cAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAAtatt                                  <  1:22580771/67‑1 (MQ=255)
                                     gCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                               >  1:25075349/1‑68 (MQ=255)
                                       gACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATattatt                               <  1:1002948/66‑1 (MQ=255)
                                       gACCACCGCTTTTTTTTATGATCATTAACACGACACCAGGGAAGAGCCGAATTCGGCAATATTATTGt                             <  1:12347047/68‑1 (MQ=255)
                                        accaccGCTTTTTTTTTTTATAATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTc                            <  1:6452534/68‑1 (MQ=255)
                                        accaccGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTc                            >  1:13043110/1‑68 (MQ=255)
                                        accaccGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTc                            >  1:39027254/1‑68 (MQ=255)
                                        accaccGCTTTTTTTTATGATAATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGt                             <  1:39460331/67‑1 (MQ=255)
                                          caccGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCAt                          <  1:23527839/68‑1 (MQ=255)
                                           accGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCAtt                         <  1:29545782/68‑1 (MQ=255)
                                            ccGCTTTTTTTTATGATCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTg                        <  1:29910650/68‑1 (MQ=255)
                                             ggCTTTTTTTTTTGAACTTTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:27393020/67‑1 (MQ=255)
                                             ggCTTTTTTTTTTCAACATTACCACGACTCCATGGAAGAGCCGTATTCGGCAATATTATTGTCATTGt                       <  1:33082211/67‑1 (MQ=255)
                                             cGCTTTTTTTTATTATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:42616038/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:45017512/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:47729659/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:48528684/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:9465277/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:13808048/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:1511160/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:27747746/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:14808542/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:26692718/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTACTGTCATTGt                       <  1:48270711/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:45727590/68‑1 (MQ=255)
                                             cGCTTTTTTTTATGATAATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGt                       <  1:31445638/68‑1 (MQ=255)
                                             cGCTTTTTTTAATGAGCATTACCACGACTCCAGGGAAGAGCCGAATGCGGCAATATTATTGTCATTGt                       <  1:39947692/68‑1 (MQ=255)
                                              gCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTa                      >  1:28799951/1‑68 (MQ=255)
                                               cTTTTTTTTTTTAATATAAACTGGACTCGAGGGAAGAGCCGAATTCGGCAATATTATTGTCTTTGTa                      <  1:46641358/67‑1 (MQ=255)
                                               cTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTAt                     >  1:16798046/1‑68 (MQ=255)
                                                ttttttttATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTAt                     <  1:10439752/67‑1 (MQ=255)
                                                ttttttttATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTAt                     <  1:5681138/67‑1 (MQ=255)
                                                 tttttttATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGa                   <  1:30759716/68‑1 (MQ=255)
                                                    ttttATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAgg                <  1:43954316/68‑1 (MQ=255)
                                                     ttttTGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGa               <  1:10768279/68‑1 (MQ=255)
                                                     tttATTATAATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGa               <  1:3514929/68‑1 (MQ=255)
                                                      ttttGATCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGat              <  1:44541381/65‑1 (MQ=255)
                                                      ttATGATCATTACCACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGat              <  1:18556376/68‑1 (MQ=255)
                                                      ttATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGat              <  1:50555828/68‑1 (MQ=255)
                                                      ttATGATAATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGat              <  1:6039729/68‑1 (MQ=255)
                                                        aTGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGatat            <  1:18074923/68‑1 (MQ=255)
                                                           aTCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:25728578/68‑1 (MQ=255)
                                                           aTCATAAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:50202443/68‑1 (MQ=255)
                                                            tCATTACCACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:8149988/67‑1 (MQ=255)
                                                            tCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:15885555/67‑1 (MQ=255)
                                                            tCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:23394639/67‑1 (MQ=255)
                                                            tCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCgg         <  1:3708842/67‑1 (MQ=255)
                                                             cATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGc       <  1:38633372/68‑1 (MQ=255)
                                                              aTTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCa      <  1:26257380/68‑1 (MQ=255)
                                                                tAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAta    <  1:23029998/68‑1 (MQ=255)
                                                                tAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAta    <  1:7531463/68‑1 (MQ=255)
                                                                tAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAta    <  1:26887099/68‑1 (MQ=255)
                                                                tAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAta    <  1:41117368/68‑1 (MQ=255)
                                                                tAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCAta    <  1:32779573/68‑1 (MQ=255)
                                                                  |                                                                   
CACTGATTACCACGCTCGCCAGCAGGCTTAGCTGCCAGCGACCACCGCTTTTTTTTATGATCATTAACACGACTCCAGGGAAGAGCCGAATTCGGCAATATTATTGTCATTGTATGAAGGATATCGGGCATAGT  >  CP009273/3706231‑3706364

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