Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,263,009 A→G 41.5% I58I (ATT→ATC cutC ← copper homeostasis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,263,0090AG41.5% 32.4 / 69.2 94I58I (ATT→ATCcutCcopper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base A (21/34);  new base G (27/12);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 3.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 101 total aligned reads displayed.

CCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAA  >  minE/1262951‑1263077
                                                          |                                                                    
ccTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:2358533/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:3043824/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:2473933/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:746065/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:1120823/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:2528118/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:3297298/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:2826755/69‑1 (MQ=255)
 cTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:3107011/69‑1 (MQ=255)
   cAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           >  1:3196843/1‑67 (MQ=255)
   cAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           >  1:1661394/1‑67 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGcac                                                        >  1:2954382/1‑65 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                    >  1:2397994/1‑69 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                    >  1:486537/1‑69 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                    >  1:2130393/1‑69 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                    >  1:264642/1‑69 (MQ=255)
        aTGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAggg                                                    >  1:3063282/1‑69 (MQ=255)
           gcggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATc                                                 >  1:146948/1‑69 (MQ=255)
           gcggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATc                                                 >  1:1884725/1‑69 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:1428726/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:2473248/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:3270037/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:2388614/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:1675369/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:1195172/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:1565699/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:3120631/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:1555962/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:3348459/69‑1 (MQ=255)
             ggcAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:3206347/69‑1 (MQ=255)
              acAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGGATAATCGGATGCACTGGGATCGt                                               <  1:3266171/67‑1 (MQ=255)
               cAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATg                                                           <  1:3269009/55‑1 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGa                                                                      >  1:769692/1‑40 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGa                                                                      >  1:640303/1‑40 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGa                                                                      >  1:1657404/1‑40 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGa                                                                      >  1:3387292/1‑40 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                         >  1:316262/1‑69 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                         >  1:201631/1‑69 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                         >  1:2659133/1‑69 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                         >  1:3261893/1‑69 (MQ=255)
                   cTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCg                                         >  1:428486/1‑69 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGt                                               >  1:857997/1‑61 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGc                                        >  1:3093077/1‑68 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       <  1:2468957/69‑1 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       >  1:1231177/1‑69 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       >  1:1715144/1‑69 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       >  1:2616729/1‑69 (MQ=255)
                     cACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       >  1:3428501/1‑69 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:334498/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:2767168/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:2684187/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:1316378/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:2632542/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:3079299/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:3287914/68‑1 (MQ=255)
                              tGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:1924900/68‑1 (MQ=255)
                               gTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:1985351/67‑1 (MQ=255)
                               gTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACg                               <  1:2318884/67‑1 (MQ=255)
                                     aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGTATTTCAg                       <  1:617351/69‑1 (MQ=255)
                                     aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                       <  1:2045688/69‑1 (MQ=255)
                                     aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                       <  1:609495/69‑1 (MQ=255)
                                     aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                       <  1:1742397/69‑1 (MQ=255)
                                     aaaaTCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                       <  1:1918423/69‑1 (MQ=255)
                                     aaaaTCACCACCGCCTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAg                       <  1:1098094/69‑1 (MQ=255)
                                             caccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACAccc                <  1:950027/68‑1 (MQ=255)
                                             caccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACAccc                <  1:1705933/68‑1 (MQ=255)
                                             caccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACAccc                <  1:1099401/68‑1 (MQ=255)
                                             caccGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACAccc                <  1:2549490/68‑1 (MQ=255)
                                             caccGCGTGGGCGAATTATCGGATGCACAGGGATCGACACCCGCTGGCGCACGGATTTCAGTACAccc                <  1:1135016/68‑1 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:300793/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:3377573/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:1557285/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:182610/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:2299263/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:700336/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:2535260/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:71049/1‑69 (MQ=255)
                                                 gcgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACCCCCAacga           >  1:2827386/1‑69 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCCGTac                    >  1:2123630/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:1131182/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:1447723/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:359758/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:1523692/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:1024513/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:153691/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:101692/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:1538839/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:189637/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:2138905/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:3256305/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:2517021/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:3107792/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:2740363/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:2931817/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTac                    >  1:2966206/1‑59 (MQ=255)
                                                  cgTGGGCGGATAATCGGATGCACTGGGATCGTCACCCGCTGGCGCACGGATTTCAGTa                     >  1:1092417/1‑58 (MQ=255)
                                                      ggCGAATTATCGGATGCACAGGGATCGTCACCCGCt                                       <  1:1901463/36‑1 (MQ=255)
                                                          aaTTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTaa  <  1:1497717/69‑1 (MQ=255)
                                                          aaTTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTaa  <  1:1155394/69‑1 (MQ=255)
                                                          aaTTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTaa  <  1:2277666/69‑1 (MQ=255)
                                                          |                                                                    
CCTCAAGAATGGCGGCAAACTCACCGTCGCTGTAGCAAAAATCACCACCGCGTGGGCGAATTATCGGATGCACAGGGATCGTCACCCGCTGGCGCACGGATTTCAGTACACCCAACGACGGCGTTAA  >  minE/1262951‑1263077

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