Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,765 G→C 10.8% intergenic (‑70/‑92) 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,7650GC10.8% 329.3 / 24.3 120intergenic (‑70/‑92)argI/yjgDornithine carbamoyltransferase 1/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (65/42);  new base C (7/6);  total (72/48)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 124 total aligned reads displayed.

AGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGC  >  minE/2919696‑2919834
                                                                     |                                                                     
aGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGt                                                                      >  1:123024/1‑71 (MQ=255)
   tCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTa                                                                   <  1:5315923/71‑1 (MQ=255)
   tCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTa                                                                   <  1:4160445/71‑1 (MQ=255)
     cctcctGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTaaa                                                                 <  1:68730/71‑1 (MQ=255)
     cctcctGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTaaa                                                                 <  1:2922780/71‑1 (MQ=255)
     cctcctGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTaaa                                                                 <  1:1618768/71‑1 (MQ=255)
     cctcctGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTaa                                                                  <  1:2596135/70‑1 (MQ=255)
         ctGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCa                                                              >  1:3638909/1‑70 (MQ=255)
         atgtgGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCa                                                              >  1:765822/2‑70 (MQ=255)
           gtgGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATa                                                            >  1:3742826/1‑70 (MQ=255)
           gtgGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATa                                                            >  1:5183378/1‑70 (MQ=255)
            tgGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATaaa                                                          >  1:1274731/1‑71 (MQ=255)
              gCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATaa                                                           >  1:3774842/1‑68 (MQ=255)
              gCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATc                                                        >  1:4369160/1‑71 (MQ=255)
               cTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATGAATCt                                                       <  1:4456394/71‑1 (MQ=255)
                    gCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                  <  1:3988443/71‑1 (MQ=255)
                       tCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCt                                                                    <  1:4184464/50‑1 (MQ=255)
                         aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATc                                                        >  1:4190779/1‑60 (MQ=255)
                         aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATc                                                        >  1:2637125/1‑60 (MQ=255)
                         aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                               >  1:3737516/1‑69 (MQ=255)
                         aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                               >  1:2115807/1‑69 (MQ=255)
                         aaTGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaaa                                             <  1:725143/71‑1 (MQ=255)
                           tGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATc                                                        <  1:4783280/58‑1 (MQ=255)
                           tGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTcaca                                               <  1:4754993/67‑1 (MQ=255)
                            gAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGt                                          >  1:2808177/1‑71 (MQ=255)
                             aaTTAATATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                         >  1:572657/1‑71 (MQ=255)
                             aaTTAAAATTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                         <  1:2476002/71‑1 (MQ=255)
                             aaTTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAg                                            >  1:1130404/1‑68 (MQ=255)
                              aTTAAAATTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                        <  1:1938506/71‑1 (MQ=255)
                               ttAAAATTCAATTTTTATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTgg                                        <  1:3479434/70‑1 (MQ=255)
                               ttAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCt                                                       <  1:3286011/55‑1 (MQ=255)
                                tAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAg                                      >  1:4909873/1‑71 (MQ=255)
                                 aaaaTTCACTTTATATGTGTAATTATTCATTTGCACCCCCATTTCACAATTCTTTCTTACAAAGGTGGAgg                                     >  1:4112108/1‑71 (MQ=255)
                                 aaaaTTCACTTTATATGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGAgg                                     >  1:5335653/1‑71 (MQ=255)
                                 aaaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTaaa                                                                 >  1:153737/1‑43 (MQ=255)
                                  aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAAt                                                         >  1:518330/1‑50 (MQ=255)
                                  aaaTTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                              >  1:3431875/1‑61 (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:3072625/1‑70 (MQ=255)
                                     ttCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                  >  1:346257/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)
                                        aaTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTc                              <  1:2953688/71‑1 (MQ=255)
                                         aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGt                                                   <  1:4381462/49‑1 (MQ=255)
                                         aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                             <  1:667106/71‑1 (MQ=255)
                                         aTTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                             <  1:496189/71‑1 (MQ=255)
                                          tttATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                         <  1:1656882/58‑1 (MQ=255)
                                          tttATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                    >  1:4966007/1‑63 (MQ=255)
                                            tataTGTGTAATTATTCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTGGAGGCaaa                                 >  1:2163919/1‑63 (MQ=255)
                                            tataTGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                          >  1:1443233/1‑71 (MQ=255)
                                             atatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCt                                                       <  1:2773216/41‑1 (MQ=255)
                                             atatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtgg                         >  1:3991974/1‑71 (MQ=255)
                                              tatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggt                        <  1:4119257/71‑1 (MQ=255)
                                              tatGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggt                        <  1:5645027/71‑1 (MQ=255)
                                                tGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTc                                                  >  1:5444318/1‑43 (MQ=255)
                                                tGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtggtg                       <  1:1991408/70‑1 (MQ=255)
                                                  gatgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                           >  1:1215417/1‑48 (MQ=255)
                                                   atgatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgtgt                    >  1:525556/1‑70 (MQ=255)
                                                     gatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                           >  1:3771574/1‑45 (MQ=255)
                                                     gatTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAgg                                           >  1:1964015/1‑45 (MQ=255)
                                                      attattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                 >  1:2177611/1‑70 (MQ=255)
                                                      attattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                        >  1:5287353/1‑47 (MQ=255)
                                                      attattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                        >  1:4554212/1‑47 (MQ=255)
                                                      attattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTgg                                        >  1:3582467/1‑47 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                              <  1:1363484/40‑1 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                              <  1:1044524/40‑1 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                  >  1:1234272/1‑52 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                  >  1:5496535/1‑52 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                  >  1:1820770/1‑52 (MQ=255)
                                                       ttattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCaa                                  >  1:2205611/1‑52 (MQ=255)
                                                       ttattCATTTGCAACCCCATTTCACAATTCTTTCTTACAAAGGTg                                         <  1:3937787/45‑1 (MQ=255)
                                                        tattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACaa                                              <  1:3648732/39‑1 (MQ=255)
                                                        tattCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                    <  1:2569352/49‑1 (MQ=255)
                                                          ttCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt              >  1:2149397/1‑69 (MQ=255)
                                                          ttCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATaa            <  1:3577876/71‑1 (MQ=255)
                                                           tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgtg                          >  1:3078145/1‑56 (MQ=255)
                                                           tCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAAt              >  1:1000203/1‑68 (MQ=255)
                                                             aTTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAagag         >  1:3890303/1‑71 (MQ=255)
                                                              tttGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCa                             <  1:1173437/50‑1 (MQ=255)
                                                               ttGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTg                                         >  1:3220703/1‑37 (MQ=255)
                                                               ttGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGgtgt                      >  1:422535/1‑56 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                           >  1:4966426/1‑50 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                           >  1:2665336/1‑50 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAgt                           >  1:3778628/1‑50 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTa       <  1:518456/70‑1 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt      >  1:5344715/1‑71 (MQ=255)
                                                                tGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTAt      >  1:2416517/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)
                                                                    aGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGc  <  1:5530268/71‑1 (MQ=255)
                                                                     gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                    <  1:4891607/36‑1 (MQ=255)
                                                                     gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGc                                    <  1:1558451/36‑1 (MQ=255)
                                                                     gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAg                            <  1:3053669/44‑1 (MQ=255)
                                                                     gTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGAc                 <  1:1005353/55‑1 (MQ=255)
                                                                     |                                                                     
AGATCCCTCCTGTGGCTAACGCCTCAATGAATTAAAATTCAATTTATATGGATGATTATTCATTTGCAAGTCTAAAGCATAAATCTTTGTCACAAAGGTGGAGGCAATGTCAGTGGTGTGTGACAATAAGAGTATCGGC  >  minE/2919696‑2919834

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