Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,786 C→T 8.5% intergenic (‑91/‑71) argI ← / → yjgD ornithine carbamoyltransferase 1/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,7860CT 8.5% 441.8 / 19.2 153intergenic (‑91/‑71)argI/yjgDornithine carbamoyltransferase 1/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (84/56);  new base T (7/6);  total (91/62)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

GCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGC  >  minE/2919716‑2919855
                                                                      |                                                                     
     aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                                                    >  1:3737516/1‑69 (MQ=255)
     aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                                  <  1:725143/71‑1 (MQ=255)
       tGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                                                    <  1:4754993/67‑1 (MQ=255)
        gAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGt                                                               >  1:2808177/1‑71 (MQ=255)
         aaTTAAAATTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                                              <  1:2476002/71‑1 (MQ=255)
          aTTAAAATTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                                              <  1:552584/70‑1 (MQ=255)
           ttAAAATTCAATTTTTATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTgg                                                             <  1:3479434/70‑1 (MQ=255)
             aaaaTTCACTTTATATGTGTAATTATTCATTTGCACCCCCATTTCACAATTCTTTCTTACAAAGGTGGAgg                                                          >  1:4112108/1‑71 (MQ=255)
              aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAgg                                                          <  1:4147958/70‑1 (MQ=255)
              aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                         <  1:4741111/71‑1 (MQ=255)
               aaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                         <  1:1688589/70‑1 (MQ=255)
                tttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                                                  <  1:5648880/59‑1 (MQ=255)
                aTTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGAGGCaa                                                       <  1:2416299/71‑1 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                       >  1:346257/1‑70 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                       >  1:3072625/1‑70 (MQ=255)
                 ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAAt                                                      >  1:5313446/1‑71 (MQ=255)
                  tCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:4259186/1‑60 (MQ=255)
                  tCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                       >  1:5297771/1‑69 (MQ=255)
                   cAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:5261381/1‑59 (MQ=255)
                     aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                                                  <  1:496189/71‑1 (MQ=255)
                     aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                                                  <  1:667106/71‑1 (MQ=255)
                        tataTGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                                               >  1:1443233/1‑71 (MQ=255)
                         atatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtgg                                              >  1:3991974/1‑71 (MQ=255)
                            tGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                                       >  1:5444318/1‑43 (MQ=255)
                            tGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                                            <  1:1991408/70‑1 (MQ=255)
                              gatgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:2127078/1‑48 (MQ=255)
                              gatgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:1215417/1‑48 (MQ=255)
                              gatgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                         >  1:3009926/1‑71 (MQ=255)
                               atgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                                           >  1:1835171/1‑52 (MQ=255)
                                 gatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:3771574/1‑45 (MQ=255)
                                 gatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                                                >  1:1964015/1‑45 (MQ=255)
                                  attattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                      >  1:2177611/1‑70 (MQ=255)
                                   ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                                       >  1:1234272/1‑52 (MQ=255)
                                   ttattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                                              <  1:3937787/45‑1 (MQ=255)
                                    tattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                                                   <  1:3648732/39‑1 (MQ=255)
                                      ttCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                                   >  1:2149397/1‑69 (MQ=255)
                                      ttCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATaa                                 <  1:3577876/71‑1 (MQ=255)
                                       tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                                   >  1:1000203/1‑68 (MQ=255)
                                         aTTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              >  1:851574/1‑71 (MQ=255)
                                         aTTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              >  1:3890303/1‑71 (MQ=255)
                                          tttGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                                                  <  1:1173437/50‑1 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                                >  1:4966426/1‑50 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                                                >  1:2665336/1‑50 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                            <  1:518456/70‑1 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt                           >  1:5344715/1‑71 (MQ=255)
                                            tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAAAAGAGTAt                           >  1:4699726/1‑71 (MQ=255)
                                             gCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                                              <  1:2487863/35‑1 (MQ=255)
                                              cAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       <  1:2657415/57‑1 (MQ=255)
                                               aaGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                                         <  1:3225505/54‑1 (MQ=255)
                                               aaGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCgg                        <  1:2293791/71‑1 (MQ=255)
                                                aGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc                       <  1:5530268/71‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                         <  1:1558451/36‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                                         <  1:4891607/36‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAg                                                 <  1:3053669/44‑1 (MQ=255)
                                                 gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                      <  1:1005353/55‑1 (MQ=255)
                                                  tCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAg                     <  1:4546108/71‑1 (MQ=255)
                                                   cTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAg                     >  1:2026998/1‑70 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:3682534/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:3378416/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:2877922/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:299969/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:1570164/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:3528147/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:5600443/1‑71 (MQ=255)
                                                    tAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   >  1:4540978/1‑71 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   <  1:4001974/69‑1 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGa                   <  1:3125529/69‑1 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGAc                  <  1:4459273/70‑1 (MQ=255)
                                                      aaGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACa                 >  1:5501872/1‑71 (MQ=255)
                                                       aGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACAt                <  1:4154781/71‑1 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt                                   <  1:1136269/48‑1 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              >  1:4209560/1‑53 (MQ=255)
                                                           tAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              >  1:4202811/1‑53 (MQ=255)
                                                            aaaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       <  1:4790631/43‑1 (MQ=255)
                                                            aaaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       <  1:4883095/43‑1 (MQ=255)
                                                            aaaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGAAGGACATTaa             <  1:3168621/69‑1 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                                            >  1:4848204/1‑37 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                                            >  1:2773370/1‑37 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa             >  1:615444/1‑68 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag          <  1:1123610/71‑1 (MQ=255)
                                                             aaTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag          <  1:4443757/71‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                                      >  1:5053565/1‑42 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag          >  1:4933423/1‑70 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAagag          >  1:1193032/1‑70 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg         <  1:1490840/71‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg         <  1:1880655/71‑1 (MQ=255)
                                                              aTCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAgg         <  1:1066473/71‑1 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       >  1:1345603/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       >  1:5594556/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGa                                       >  1:3165966/1‑39 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACa                                     >  1:351258/1‑41 (MQ=255)
                                                                cTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGaa       >  1:1812894/1‑71 (MQ=255)
                                                                 tttGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              <  1:1536260/47‑1 (MQ=255)
                                                                 tttGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              <  1:2573747/47‑1 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag                              >  1:730605/1‑45 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                            >  1:3804617/1‑47 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa                            >  1:4772066/1‑47 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTaa             <  1:5169190/62‑1 (MQ=255)
                                                                   tGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAAtt     <  1:3418197/70‑2 (MQ=255)
                                                                     tCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGc  <  1:5083344/71‑1 (MQ=255)
                                                                      |                                                                     
GCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGCAGGACATTAAGAGGAATGAGC  >  minE/2919716‑2919855

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