Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,608,471 T→A L674Q (CTG→CAG)  zntA → zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,4710TA97.7% 150.8 / ‑5.0 44L674Q (CTG→CAG) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (26/17);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGG  >  NC_000913/3608380‑3608567
                                                                                           |                                                                                                
gCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCAcgcg                                                                                          >  1:531225/1‑100 (MQ=255)
gCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCAcgcg                                                                                          >  1:512942/1‑100 (MQ=255)
  aaTGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGcc                                                                                        >  1:1283657/1‑100 (MQ=255)
      ggTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTc                                                                                    >  2:1081429/1‑100 (MQ=255)
          gCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGc                                                                                >  1:19007/1‑100 (MQ=255)
          gCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGAGCCAATGATTGAACAGGCACGCGCCACTCACGc                                                                                >  2:904676/1‑100 (MQ=255)
           cGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGcc                                                                               <  2:1251329/100‑1 (MQ=255)
            ggCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCa                                                                              >  2:966503/1‑100 (MQ=255)
              cacaGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAat                                                                            <  1:216779/100‑1 (MQ=255)
                   aCGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCg                                                                       <  1:1255742/100‑1 (MQ=255)
                        gcgcTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGa                                                                  >  1:735520/1‑100 (MQ=255)
                           cTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACa                                                               <  1:158789/100‑1 (MQ=255)
                            tGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACAt                                                              <  1:633147/100‑1 (MQ=255)
                             ggAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATc                                                             >  2:1052066/1‑100 (MQ=255)
                                   cgccgACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTAtt                                                       >  1:1302140/1‑100 (MQ=255)
                                          gcagcaTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCGCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTgg                                                <  2:725588/100‑1 (MQ=255)
                                               aTTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTg                                           >  2:1048304/1‑100 (MQ=255)
                                                ttAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGa                                          <  1:638304/100‑1 (MQ=255)
                                                 tAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTgga                                         <  1:881586/100‑3 (MQ=255)
                                                 tAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGaa                                         <  2:662289/100‑1 (MQ=255)
                                                   aCCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAg                                       <  2:512938/100‑1 (MQ=255)
                                                    cccATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                      >  2:37844/1‑100 (MQ=255)
                                                    cccATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                      >  2:318307/1‑100 (MQ=255)
                                                    cccATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                      >  1:1035132/1‑100 (MQ=255)
                                                    cccATAACAACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGACAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAgg                                      <  2:1110382/100‑1 (MQ=255)
                                                      cATAACCGCCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACGCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGa                                    <  2:218642/100‑1 (MQ=255)
                                                      cATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGa                                    <  1:1029568/100‑1 (MQ=255)
                                                         aaCCACCTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCt                                 >  1:286822/1‑100 (MQ=255)
                                                          accaccTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCtt                                >  1:601564/1‑100 (MQ=255)
                                                            caccTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATGTCCGCCAGAACTTCACTATTGCGCTGGGGCTAAAAGGGATCTTcc                              <  2:536765/100‑1 (MQ=255)
                                                             accTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATGGCGCTGGGGCGGAGAGGGACCTTcc                              >  1:191935/1‑99 (MQ=255)
                                                               cTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCg                           <  2:1098977/100‑1 (MQ=255)
                                                               cTGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCg                           >  1:366193/1‑100 (MQ=255)
                                                                tGCGCGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGt                          >  2:840330/1‑100 (MQ=255)
                                                                 gcgcGGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTc                         >  1:259311/1‑100 (MQ=255)
                                                                 gcgcGACCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTc                         >  1:1304024/1‑100 (MQ=255)
                                                                    cgGCCTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTcacc                      >  2:693482/1‑100 (MQ=255)
                                                                        cTGGTGCAAATGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGc                  >  1:552086/1‑100 (MQ=255)
                                                                                aaTGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCCCTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAggg          <  2:566614/100‑1 (MQ=255)
                                                                                aaTGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAggg          <  2:775505/100‑1 (MQ=255)
                                                                                  tGATTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGAt        <  2:1283646/100‑1 (MQ=255)
                                                                                    aTTGAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGa      >  2:653236/1‑100 (MQ=255)
                                                                                       gAACAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGAATCTTCCTCGTCACCACGCTGTTAGGGATGAc     >  2:273615/1‑98 (MQ=255)
                                                                                        aaCAGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCgg  >  2:939612/1‑100 (MQ=255)
                                                                                           |                                                                                                
GCAATGGGTAGCGGCACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGG  >  NC_000913/3608380‑3608567

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