Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,023,368 A→G R172R (CGT→CGC yccS ← putative transporter YccS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,023,3680AG97.4% 120.6 / ‑4.2 39R172R (CGT→CGCyccSputative transporter YccS
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/17);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTGAT  >  NC_000913/1023238‑1023494
                                                                                                                                  |                                                                                                                              
gATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGcg                                                                                                                              >  1:46929/1‑133 (MQ=255)
  ttCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGcgcg                                                                                                                            <  2:46929/133‑1 (MQ=255)
           gCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGt                                                                                                                   >  1:363602/1‑133 (MQ=255)
                     gACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAgcgg                                                                                                         >  1:480684/1‑133 (MQ=255)
                      aCCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAgcggg                                                                                                        >  1:495697/1‑133 (MQ=255)
                        cGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCCCCACGTAGTCCTGTAgcgggcg                                                                                                      >  1:209692/1‑133 (MQ=255)
                          ttGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGa                                                                                                    <  1:451481/133‑1 (MQ=255)
                                         tCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGAAACAGCAGAt                                                                                     <  2:471342/133‑1 (MQ=255)
                                            tccAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCGAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGAc                                                                                  <  2:391523/131‑1 (MQ=255)
                                               aGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCaa                                                                               <  1:479695/133‑1 (MQ=255)
                                                   gTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAg                                                                           >  2:98133/1‑133 (MQ=255)
                                                      ccTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACTGGGAACAGCATATGAACAATAAgtgt                                                                        >  2:229522/1‑133 (MQ=255)
                                                         ggCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAg                                                                     <  2:56106/133‑1 (MQ=255)
                                                           cTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGa                                                                   >  1:383483/1‑133 (MQ=255)
                                                              tGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGCGCGGACCGAGAAGAGCAGAGGACCAATAAGGTTCAGGAcct                                                                >  2:246450/1‑131 (MQ=255)
                                                                        tatCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGAc                                                      >  2:286727/1‑133 (MQ=255)
                                                                          tCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACgg                                                    >  2:62335/1‑133 (MQ=255)
                                                                          tCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACgg                                                    >  1:287593/1‑133 (MQ=255)
                                                                               aTCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCg                                               >  1:272946/1‑133 (MQ=255)
                                                                               aTCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCg                                               >  1:16306/1‑133 (MQ=255)
                                                                                 cAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGc                                             <  2:340109/133‑1 (MQ=255)
                                                                                  aaaCATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGcc                                            <  2:363602/133‑1 (MQ=255)
                                                                                  aaaCATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGcc                                            <  2:495697/133‑1 (MQ=255)
                                                                                  aaaCATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGAAGTTGTACCAGACGGCACCGGcc                                            >  1:184031/1‑133 (MQ=255)
                                                                                    aCATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCcag                                          >  2:362603/1‑133 (MQ=255)
                                                                                        gcgcGACTTGAGCTCAAGATAACGCGCCAGATGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGTGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGa                                      >  2:434490/1‑133 (MQ=255)
                                                                                         cgcgACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGAt                                     >  1:89236/1‑133 (MQ=255)
                                                                                         cgcgACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGAt                                     >  1:411063/1‑133 (MQ=255)
                                                                                          gcgACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATa                                    >  2:499330/1‑133 (MQ=255)
                                                                                              cTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATc                                >  1:5561/1‑133 (MQ=255)
                                                                                                tGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCgg                              <  1:21705/133‑1 (MQ=255)
                                                                                                  aGCTCAAGATAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGgct                            <  2:89236/133‑1 (MQ=255)
                                                                                                           tAACGCGCCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATAcc                   <  1:98133/133‑1 (MQ=255)
                                                                                                               gcgcCAGTTGTTCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTg               <  2:446393/133‑1 (MQ=255)
                                                                                                                   cAGTTGTTCATAGCATCGCGCCAGTTTGTCCTGCAGCGGGCCGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCa           <  2:209692/133‑1 (MQ=255)
                                                                                                                          tCATAGCCGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTg    <  2:480684/133‑1 (MQ=255)
                                                                                                                          tCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTg    <  2:287593/133‑1 (MQ=255)
                                                                                                                          tCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTg    <  2:418103/133‑1 (MQ=255)
                                                                                                                          tCATAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTg    <  1:362603/133‑1 (MQ=255)
                                                                                                                            aTAGCAGCGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTGAt  <  1:286727/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                              
GATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGATACCAGTGCTCATACAGTGAT  >  NC_000913/1023238‑1023494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: