Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,274,879 T→A 14.3% H198Q (CAT→CAA eamA → O‑acetylserine/cysteine export transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,274,8790TA14.3% 87.1 / 8.3 35H198Q (CAT→CAAeamAO‑acetylserine/cysteine export transporter
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base A (3/2);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGT  >  NC_002947/3274734‑3275025
                                                                                                                                                 |                                                                                                                                                  
gTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCCCTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCt                                                                                                                                                >  2:751254/1‑150 (MQ=255)
       gtgGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCaa                                                                                                                                         <  2:485344/150‑1 (MQ=255)
         ggggCCTCGCCAACATCATTACCCGTCGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAata                                                                                                                                       <  2:281226/150‑1 (MQ=255)
            gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATAACGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCa                                                                                                                                    <  2:589579/150‑1 (MQ=255)
            gCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCGGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCa                                                                                                                                    <  1:751254/150‑1 (MQ=255)
             ccTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAg                                                                                                                                   >  2:557920/1‑150 (MQ=255)
                 ggCAACATCATTACCCGACGGTACGGCTCGGTCTACGTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTgg                                                                                                                               >  1:19331/1‑150 (MQ=255)
                  gCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGa                                                                                                                              <  2:19331/150‑1 (MQ=255)
                       atcatTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCg                                                                                                                         >  1:24246/1‑150 (MQ=255)
                            tACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCt                                                                                                                    >  2:742981/1‑150 (MQ=255)
                                   cGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTg                                                                                                             <  1:253091/150‑1 (MQ=255)
                                     gTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGGTGGCCCTGGc                                                                                                           >  1:176098/1‑150 (MQ=255)
                                           cTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTAttt                                                                                                     <  1:557920/150‑1 (MQ=255)
                                                    ccTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTtgt                                                                                            <  2:577574/150‑1 (MQ=255)
                                                    ccTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTtgt                                                                                            <  1:576136/150‑1 (MQ=255)
                                                     cTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCAat                                                                                                                                        >  2:246636/1‑105 (MQ=255)
                                                     cTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCAAGGCATGGCCAat                                                                                                                                        <  1:246636/105‑1 (MQ=255)
                                                                     tCTGGGGTGGGCTGATACCGCCACTTCNATTCCTGGCACTNTCCTGGNGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCt                                                                           >  1:545818/1‑150 (MQ=255)
                                                                      cTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCAATAGCGTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTa                                                                          <  1:58239/150‑1 (MQ=255)
                                                                               gCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGt                                                                                                                                  <  1:659442/85‑1 (MQ=255)
                                                                               gCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGt                                                                                                                                  >  2:659442/1‑85 (MQ=255)
                                                                                             ttCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGt                                                   <  2:24246/150‑1 (MQ=255)
                                                                                                cATTCCTGGCACTGTCCTGGTGCATGGACGGCCCAGTGTGCATTGCCCATGCCTCAGCCCATCTCTTAAAGAGTTTGGTGGTGGTGTGGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTcgc                                                <  2:816797/150‑1 (MQ=255)
                                                                                                     cTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCAc                                           >  2:360242/1‑150 (MQ=255)
                                                                                                                        tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtgtg                                                                 >  2:282160/1‑109 (MQ=255)
                                                                                                                        tGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCtgtg                                                                 <  1:282160/109‑1 (MQ=255)
                                                                                                                         ggAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGcc                       >  2:191968/1‑150 (MQ=255)
                                                                                                                                        cATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGt        >  2:720999/1‑150 (MQ=255)
                                                                                                                                         aTTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTa       >  2:579952/1‑150 (MQ=255)
                                                                                                                                         aTTGGCCAAGGCATGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTa       >  2:680893/1‑150 (MQ=255)
                                                                                                                                         aTTGGCCAAGGCATGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTa       >  2:693366/1‑150 (MQ=255)
                                                                                                                                             gCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGt                             >  2:6568/1‑124 (MQ=255)
                                                                                                                                             gCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGt                             <  1:6568/124‑1 (MQ=255)
                                                                                                                                             gCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCACTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCgg   >  2:766760/1‑150 (MQ=255)
                                                                                                                                              ccATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGt  >  1:2958/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTGCGATGTGGGGCCTGGGCAACATCATTACCCGACGGTTCGGCTCGGTCGACCTGGTCGCGTTGGTGATCTGGGGTGGGCTGATACCGCCACTTCCATTCCTGGCACTGTCCTGGTGGCTGGAAGGCCCCGAGCGCATTGGCCATGCCTTGGCCAATATCAGTTGGAGTTCGGTGCTGGCCCTGGCGTATTTGGCCTTTGTGGCCACCATGCTCGGCTACAGCCTGTGGAGCAAGTTGCTGTCGCGTCACCCGGCAGGCAAAGTGGCGCCGTTCTCGCTGCTGGTACCGGT  >  NC_002947/3274734‑3275025

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: