Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 18,192 C→A 10.4% intergenic (+31/+301) nhaR → / ← rpsT DNA‑binding transcriptional activator/30S ribosomal subunit protein S20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE18,1920CA10.4% 208.0 / 7.3 67intergenic (+31/+301)nhaR/rpsTDNA‑binding transcriptional activator/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base C (24/36);  new base A (4/3);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACA  >  minE/18131‑18260
                                                             |                                                                    
ttCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCaaa                                                              <  1:4199485/70‑1 (MQ=255)
          tttAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACcatc                                                     <  1:3612578/69‑1 (MQ=255)
                   gCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAAt                                          <  1:994125/71‑1 (MQ=255)
                   gCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAAt                                          <  1:4388896/71‑1 (MQ=255)
                       tGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACcatc                                                     <  1:1031055/56‑1 (MQ=255)
                           ttAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGt                                  <  1:909484/71‑1 (MQ=255)
                             aaTCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTg                                >  1:5665668/1‑71 (MQ=255)
                              aTCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTgg                               <  1:4354530/71‑1 (MQ=255)
                              aTCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTgg                               <  1:4007014/71‑1 (MQ=255)
                               tCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTaa                                            <  1:1640604/57‑1 (MQ=255)
                                      cTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAACCACCATCATACACTAAATCAGTAAGTTGGCAGCATCa                       >  1:1008144/1‑71 (MQ=255)
                                        cccccAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATa                                                  >  1:2704915/1‑42 (MQ=255)
                                         ccccAAAGTTAAGGTTGGAGCACCTCAAAAACACcatca                                                    <  1:4319724/39‑1 (MQ=255)
                                           ccAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctac                  >  1:255250/1‑71 (MQ=255)
                                              aaGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTc               <  1:4907359/71‑1 (MQ=255)
                                               aGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTc               >  1:1682707/1‑70 (MQ=255)
                                                      actggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCaaaaaagcgactca       <  1:5137758‑M1/69‑15 (MQ=255)
                                                      actggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCaaaaaagcgactca       >  1:2263577‑M1/3‑57 (MQ=255)
                                                      actGGAGCACCTCAAAAACACCATCATACACTAAATCAGt                                      <  1:561254/38‑1 (MQ=255)
                                                      actGGAGCACCTCAAAAACACCATCATACACTAAATCAGt                                      <  1:2849870/38‑1 (MQ=255)
                                                      actGGAGCACCTCAAAAACACCATCATACACTAAATCAGt                                      <  1:5634696/38‑1 (MQ=255)
                                                      actGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgca                             <  1:1946613/47‑1 (MQ=255)
                                                      actGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgca                             <  1:5346513/47‑1 (MQ=255)
                                                      actGGAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGtt                                 <  1:3315589/43‑1 (MQ=255)
                                                       ctGGAGCACCTCAAAAACACCATCATACACTAAATCAg                                       <  1:1532448/37‑1 (MQ=255)
                                                       ctGGAGCACCTCAAAAACACCATCATACACTAAATCAg                                       <  1:4652186/37‑1 (MQ=255)
                                                       ctGGAGCACCTCAAAAACACCATCATACACTAAATCAg                                       <  1:395446/37‑1 (MQ=255)
                                                       ctGGAGCACCTCAAAAACACCATCATACACTAAATCAGTaa                                    <  1:2163445/40‑1 (MQ=255)
                                                       ctGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGt                                  <  1:942139/42‑1 (MQ=255)
                                                       ctGGAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCAt                         <  1:2326418/51‑1 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTaa                                    >  1:2958711/1‑40 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCAt                         <  1:3124528/51‑1 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCAt                         <  1:5616841/51‑1 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCaccttt                  <  1:2594360/58‑3 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCac                      >  1:4495373/1‑54 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCac                      >  1:3415163/1‑54 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCaaaa                 <  1:2905206‑M1/59‑5 (MQ=255)
                                                        tGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCaaaa                 <  1:2093994‑M1/59‑5 (MQ=255)
                                                        tGGAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgc                              <  1:4214973/46‑1 (MQ=255)
                                                        tGGAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCcgtattgtccatgcga     >  1:5010402‑M1/1‑55 (MQ=255)
                                                        tGGAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACAcgtattgtccatgcga     >  1:5244862‑M1/1‑55 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACTCTAAATCAGTAAGTTGGCAGCATCACCaaaaaagcgactcacat    >  1:4556855‑M1/1‑54 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCa                                        <  1:528705/35‑1 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTaa                                    >  1:4744654/1‑39 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGt                                  >  1:3385409/1‑41 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTg                                >  1:1922725/1‑43 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTg                                >  1:2543590/1‑43 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgcag                            <  1:3822309/47‑1 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCAt                         >  1:1370068/1‑50 (MQ=255)
                                                         ggAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCAt                         >  1:5087899/1‑50 (MQ=255)
                                                         ggAGAACCTCAAAAACACCATCATACACTAAATCAGTa                                     >  1:4631326/1‑38 (MQ=255)
                                                         ggAGAACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCa                       >  1:1743335/1‑52 (MQ=255)
                                                          gagaACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgcagca                          >  1:4527828/1‑48 (MQ=255)
                                                          gagaACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgcagca                          >  1:3083399/1‑48 (MQ=255)
                                                          gagaACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCcgtattgtccatgcgatg   <  1:2711398‑M1/71‑19 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAACTCAGTAAGTTGGCAGCATCCCCaaaaaagcgactcaca     <  1:4669689‑M1/69‑17 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAg                                   >  1:5228150/1‑39 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAg                                   >  1:5334881/1‑39 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAg                                   >  1:5315803/1‑39 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGt                                  >  1:285678/1‑40 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTgg                               >  1:1502344/1‑43 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTgg                               >  1:1593612/1‑43 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGgca                             >  1:5284182/1‑45 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCaaaaaagcgactcaca     <  1:5532234‑M1/69‑17 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATcac   <  1:76182/71‑1 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATc     <  1:5336068/69‑1 (MQ=255)
                                                          gAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATc     <  1:4497109/69‑1 (MQ=255)
                                                          gAGCACCTCAAAAACACAATCATACACTAAATCAGTaa                                    >  1:4313716/1‑38 (MQ=255)
                                                           aGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTgg                               <  1:541411/42‑1 (MQ=255)
                                                           aGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATcaca  <  1:3657564/71‑1 (MQ=255)
                                                             cACCTCAAAAACACCATCATACACTAAATCAGTAAGtt                                 >  1:2973373/1‑38 (MQ=255)
                                                             cACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATc                        >  1:5546148/1‑47 (MQ=255)
                                                             |                                                                    
TTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACA  >  minE/18131‑18260

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: