Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,291,126 G→A 16.3% intergenic (+241/+212) yeeN → / ← asnW conserved hypothetical protein/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,291,1260GA16.3% 125.4 / 14.9 49intergenic (+241/+212)yeeN/asnWconserved hypothetical protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base A (2/6);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGA  >  minE/1291059‑1291187
                                                                   |                                                             
tACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTAc                                                            >  1:1090032/1‑71 (MQ=255)
 aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTtata                                                             >  1:510611/1‑69 (MQ=255)
 aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACg                                                           >  1:2060597/1‑71 (MQ=255)
      tAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTAc                                                       <  1:21546/70‑1 (MQ=255)
      tAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTAc                                                       >  1:32409/1‑70 (MQ=255)
        gTTTCCCGATAAAGACGCTATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACtt                                                     >  1:1810901/1‑70 (MQ=255)
        gTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACtt                                                     <  1:2041757/70‑1 (MQ=255)
        gTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACtt                                                     >  1:1736001/1‑70 (MQ=255)
         tttCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCa                                                   >  1:2212776/1‑71 (MQ=255)
                aTAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACa                                            <  1:1265837/71‑1 (MQ=255)
                          gATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTAc                                                            <  1:115715/45‑1 (MQ=255)
                          gATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTAc                                                            <  1:1865195/45‑1 (MQ=255)
                          gATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGt                                                         <  1:1019417/48‑1 (MQ=255)
                            tAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACgg                                                          >  1:1201942/1‑45 (MQ=255)
                               cGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTg                                   <  1:509118/65‑1 (MQ=255)
                                   aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACgg                                                          <  1:842077/38‑1 (MQ=255)
                                   aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACgg                                                          <  1:1899646/38‑1 (MQ=255)
                                   aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACgg                                                          <  1:2024211/38‑1 (MQ=255)
                                   aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                         <  1:743328/71‑1 (MQ=255)
                                   aTCCCGCATAAGGAATTTTTCTCTTTATCTTTGTACgg                                                          <  1:1506422/38‑1 (MQ=255)
                                    tCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCtt                              <  1:2071340/65‑1 (MQ=255)
                                    tCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGaa                                              <  1:454761/49‑1 (MQ=255)
                                      ccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGAcc                                >  1:1858906/1‑61 (MQ=255)
                                       cGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                     >  1:422522/1‑71 (MQ=255)
                                       cGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                     <  1:1870236/71‑1 (MQ=255)
                                         cATAAGTCATTTTTCTCTTTATCTTTGTACGGTACTTCa                                                   <  1:347321/39‑1 (MQ=255)
                                         cATAAGGCATTTTTCTCTTTATCTTTGTACGGTACtt                                                     <  1:1322790/37‑1 (MQ=255)
                                           tAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGaa                                              >  1:1308440/1‑42 (MQ=255)
                                            aaGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa                   >  1:1960951/1‑68 (MQ=255)
                                             aGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTc               <  1:931697/71‑1 (MQ=255)
                                               gCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg             <  1:1951452/71‑1 (MQ=255)
                                                cATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGt            >  1:506842/1‑71 (MQ=255)
                                                  tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                     <  1:1539108/60‑1 (MQ=255)
                                                  tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGt            <  1:2086261/69‑1 (MQ=255)
                                                  tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg          <  1:840186/71‑1 (MQ=255)
                                                  tttttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGtt           <  1:2075641/70‑1 (MQ=255)
                                                   ttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                     <  1:365737/59‑1 (MQ=255)
                                                   ttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                     <  1:189591/59‑1 (MQ=255)
                                                   ttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg          <  1:266453/70‑1 (MQ=255)
                                                   ttttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg         >  1:731967/1‑71 (MQ=255)
                                                     ttCTCTTTATCTTTGTACGTTACTTCATGGAACAGAGttt                                      >  1:489183/1‑40 (MQ=255)
                                                     ttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa                   >  1:1364190/1‑59 (MQ=255)
                                                          tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGa                                  >  1:512693/1‑39 (MQ=255)
                                                          tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGgggaggga  >  1:1931186/1‑71 (MQ=255)
                                                          tttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGgggagg    <  1:1325900/69‑1 (MQ=255)
                                                            tATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg        <  1:600280/63‑1 (MQ=255)
                                                            tATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg          >  1:730115/1‑61 (MQ=255)
                                                              tCTTTGTACGGTACTTCATGGAACAGAGTTTTTGCCCTTGCGAATCGTGATGTCTGTTGggga      >  1:1933247/1‑63 (MQ=255)
                                                              tCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGggga      >  1:533573/1‑63 (MQ=255)
                                                                tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGt                <  1:159227/51‑1 (MQ=255)
                                                                   |                                                             
TACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGA  >  minE/1291059‑1291187

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