Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,817,482 A→G 7.1% S10S (AGT→AGC ygaD ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,817,4820AG 7.1% 234.7 / 6.5 84S10S (AGT→AGCygaDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (4/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

CCACCGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAG  >  minE/1817414‑1817552
                                                                    |                                                                      
ccaccGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACt                                                                      >  1:1945674/1‑71 (MQ=255)
 caccGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACtt                                                                     <  1:5094713/71‑1 (MQ=255)
 caccGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACtt                                                                     <  1:729175/71‑1 (MQ=255)
  gccGGTACAAGACTCAGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCGCTTa                                                                    <  1:5081849/70‑1 (MQ=255)
  gccGGTACAAGACTCAGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCGCTTa                                                                    <  1:3611963/70‑1 (MQ=255)
  accGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACtt                                                                     >  1:2697397/1‑70 (MQ=255)
            gACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAg                                                          >  1:88515/1‑71 (MQ=255)
            gACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAg                                                          <  1:1498194/71‑1 (MQ=255)
             aCTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGt                                                         >  1:1589122/1‑71 (MQ=255)
              cTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGtt                                                        <  1:1491087/71‑1 (MQ=255)
                cGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTc                                                       >  1:4456877/1‑70 (MQ=255)
                cGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCa                                                      >  1:4585644/1‑71 (MQ=255)
                 ggCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTc                                                       <  1:3638865/69‑1 (MQ=255)
                  gCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACtt                                                                     <  1:2238843/54‑1 (MQ=255)
                    aGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgt                                                  <  1:3107937/71‑1 (MQ=255)
                    aGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgt                                                  >  1:4917766/1‑71 (MQ=255)
                      ttgttACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgt                                                  >  1:3768592/1‑69 (MQ=255)
                      ttgttACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgt                                                  >  1:4164986/1‑69 (MQ=255)
                          tACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgtcagtca                                            >  1:4634339/1‑71 (MQ=255)
                          tACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgtcagtca                                            >  1:3772296/1‑71 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACt                                                                      >  1:2191415/1‑44 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCa                                                           >  1:1362065/1‑55 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgtcagtca                                            <  1:4546497/70‑1 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgtcagtca                                            <  1:2502923/70‑1 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgtcagtca                                            <  1:3391856/70‑1 (MQ=255)
                           aCGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGtcat                                           >  1:1710665/1‑71 (MQ=255)
                            cGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGtcatc                                          <  1:3073998/71‑1 (MQ=255)
                                     cACGTGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATc                                                            <  1:4212194/44‑1 (MQ=255)
                                     cACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATc                                                            <  1:3798727/44‑1 (MQ=255)
                                     cACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATc                                                            <  1:2288039/44‑1 (MQ=255)
                                     cACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACctct                                 <  1:3377661/71‑1 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTaa                                                                   >  1:3287898/1‑35 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGc                                                               >  1:5304894/1‑39 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTgt                                                  <  1:4310777/52‑1 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACctctct                               >  1:1777/1‑69 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgt                               >  1:4535954/1‑71 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgt                               <  1:498416/71‑1 (MQ=255)
                                       cGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgt                               >  1:858791/1‑71 (MQ=255)
                                        gggCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtc                              <  1:4162979/71‑1 (MQ=255)
                                          gCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtcag                            <  1:5400440/71‑1 (MQ=255)
                                          gCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGGCAGTCATCACAACCTCTgtcag                            <  1:2653939/71‑1 (MQ=255)
                                            tttCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtcagt                           >  1:5230728/1‑70 (MQ=255)
                                              tCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtcagtca                         >  1:3506762/1‑70 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCAc                                                     >  1:97389/1‑40 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtca                             >  1:3477621/1‑64 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtca                             >  1:2652119/1‑64 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       >  1:5544554/1‑70 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       <  1:3725165/70‑1 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       >  1:3508128/1‑70 (MQ=255)
                                                aGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       <  1:2584842/70‑1 (MQ=255)
                                                    ccTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAAtt                   <  1:954435/70‑1 (MQ=255)
                                                     cTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAAtt                   <  1:2491278/69‑1 (MQ=255)
                                                       gCCCAACCTGTTCGCTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:2141382/1‑45 (MQ=255)
                                                       gCCCAACCTGTTCGCTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:1853131/1‑45 (MQ=255)
                                                       gCCCAACCTGTTCGCTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:4597250/1‑45 (MQ=255)
                                                       gCCCAACCTGTTCGCTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:1837279/1‑45 (MQ=255)
                                                       gCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATcac                                        <  1:481672/46‑1 (MQ=255)
                                                       gCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       <  1:5524255/63‑1 (MQ=255)
                                                       gCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCaaa                       <  1:1976377/63‑1 (MQ=255)
                                                        cccAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTAt                <  1:4945622/69‑1 (MQ=255)
                                                        cccAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGc              <  1:510872/71‑1 (MQ=255)
                                                         ccAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACaa                                      >  1:2891178/1‑46 (MQ=255)
                                                         ccAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACaa                                      >  1:4413301/1‑46 (MQ=255)
                                                         ccAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACaa                                      >  1:2193565/1‑46 (MQ=255)
                                                           aaCCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAAtt                   <  1:1273244/63‑1 (MQ=255)
                                                           aaCCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCg            >  1:697977/1‑70 (MQ=255)
                                                            aCCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAAcc                                    <  1:2049076/45‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:800129/1‑39 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:1586300/1‑39 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGtcatca                                         >  1:3731144/1‑39 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTg                                <  1:88287/48‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTg                                <  1:2215694/48‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTg                                <  1:1169882/48‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGAc         <  1:187840/71‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGAc         <  1:5654878/71‑1 (MQ=255)
                                                             ccTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGAc         <  1:4138990/71‑1 (MQ=255)
                                                              cTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACa        <  1:3860814/71‑1 (MQ=255)
                                                               tGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTgtc                              >  1:2870180/1‑48 (MQ=255)
                                                                gTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGAc         >  1:3876611/1‑68 (MQ=255)
                                                                 ttCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGAc         >  1:2317587/1‑67 (MQ=255)
                                                                 ttCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGa     <  1:3511314/71‑1 (MQ=255)
                                                                  tCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAAtt                   <  1:709476/56‑1 (MQ=255)
                                                                   cACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACaa       >  1:117198/1‑67 (MQ=255)
                                                                   cACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGa     >  1:3912824/1‑69 (MQ=255)
                                                                    aCTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAg  >  1:883743/1‑71 (MQ=255)
                                                                    |                                                                      
CCACCGGTACAAGACTCGGCAGTTGTTACGGTTGCGCCACGGGCTTTCAGCGCCTGCCCAACCTGTTCACTTAACTGCATCAGTTCACTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAG  >  minE/1817414‑1817552

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: