Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,993,572 G→A 24.5% V289M (GTG→ATG)  nupG → nucleoside transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,993,5720GA24.5% 101.9 / 30.4 53V289M (GTG→ATG) nupGnucleoside transporter
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base A (9/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGG  >  minE/1993503‑1993641
                                                                     |                                                                     
cTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgt                                                                      >  1:1053670/1‑71 (MQ=255)
cTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgt                                                                      >  1:135429/1‑71 (MQ=255)
  gTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTgg                                                                    >  1:888908/1‑71 (MQ=255)
       ttCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTg                                                                >  1:1992019/1‑70 (MQ=255)
        tCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              >  1:1285117/1‑71 (MQ=255)
        tCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              >  1:1685703/1‑71 (MQ=255)
         cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCg                                                                  >  1:863574/1‑66 (MQ=255)
         cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              <  1:595192/70‑1 (MQ=255)
         cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              <  1:2160796/70‑1 (MQ=255)
              cATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTg                                                         <  1:123800/70‑1 (MQ=255)
              cATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTg                                                         <  1:2196343/70‑1 (MQ=255)
                tCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGt                                                      >  1:1855538/1‑71 (MQ=255)
                    gTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGa                                                              <  1:2014886/59‑1 (MQ=255)
                      tcttctTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTgc                                                 >  1:980779/1‑70 (MQ=255)
                          cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCt                                              <  1:813331/69‑1 (MQ=255)
                          cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGt                                            >  1:2016985/1‑71 (MQ=255)
                          cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCATGCGTTTTGCGCt                                              <  1:182077/69‑1 (MQ=255)
                          cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCt                                              <  1:219214/69‑1 (MQ=255)
                            tAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTTTTGTGGCg                                                                  <  1:1589937/47‑1 (MQ=255)
                            tAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGttt                                          >  1:885331/1‑71 (MQ=255)
                             aaGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGt                                            <  1:920758/68‑1 (MQ=255)
                              aGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                                         <  1:1947564/70‑1 (MQ=255)
                               gCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCt                                       <  1:1587855/71‑1 (MQ=255)
                                  gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:614656/71‑1 (MQ=255)
                                  gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:1477584/71‑1 (MQ=255)
                                     aCGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 >  1:55835/1‑71 (MQ=255)
                                     aCGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 >  1:1483545/1‑71 (MQ=255)
                                      cGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 <  1:97382/70‑1 (MQ=255)
                                          aTTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:365416/63‑1 (MQ=255)
                                           ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGc                                        <  1:1751601/58‑1 (MQ=255)
                                           ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                            <  1:1646911/70‑1 (MQ=255)
                                           ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                            <  1:1399079/70‑1 (MQ=255)
                                             aaGAACTTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:285040/1‑60 (MQ=255)
                                             aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTg                                             >  1:721024/1‑51 (MQ=255)
                                             aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:511178/1‑60 (MQ=255)
                                             aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:98346/1‑60 (MQ=255)
                                             aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:501167/1‑60 (MQ=255)
                                              aGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                         <  1:1511725/70‑1 (MQ=255)
                                                   gTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   >  1:1102717/1‑71 (MQ=255)
                                                   gTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   >  1:1334015/1‑71 (MQ=255)
                                                    tAATGATGATCAGTATTGTGGCGTGGATCCTGCGttt                                                    >  1:1767181/1‑37 (MQ=255)
                                                    tAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                  <  1:2152439/71‑1 (MQ=255)
                                                      atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:180751/71‑1 (MQ=255)
                                                      atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:1245958/71‑1 (MQ=255)
                                                           gatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactg           >  1:685203/1‑71 (MQ=255)
                                                            atCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactg           <  1:1922535/70‑1 (MQ=255)
                                                              cAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAACCGACTCCGTTc                  >  1:1653750/1‑61 (MQ=255)
                                                               aGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact       >  1:954082/1‑71 (MQ=255)
                                                                 tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:580458/60‑1 (MQ=255)
                                                                  aTTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg                 >  1:655349/1‑58 (MQ=255)
                                                                   ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   <  1:906456/55‑1 (MQ=255)
                                                                   ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                >  1:785270/1‑58 (MQ=255)
                                                                   ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                >  1:1968654/1‑58 (MQ=255)
                                                                   ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgctg   >  1:167178/1‑71 (MQ=255)
                                                                   ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgctg   >  1:1369352/1‑71 (MQ=255)
                                                                     gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTgg  >  1:2114310/1‑70 (MQ=255)
                                                                     |                                                                     
CTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGG  >  minE/1993503‑1993641

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: