Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,439,996 G→A R2C (CGC→TGC)  rplQ ← 50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,9960GA97.7% 166.6 / NA 44R2C (CGC→TGC) rplQ50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/24);  minor base C (0/1);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTT  >  NC_000913/3439851‑3440141
                                                                                                                                                 |                                                                                                                                                 
ggCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATg                                                                                                                                               <  1:604847/150‑1 (MQ=255)
      aCTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTAt                                                                                                                                        >  1:198254/1‑151 (MQ=255)
                 cAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTa                                                                                                                             <  1:19306/151‑1 (MQ=255)
                  aGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTa                                                                                                                             >  2:607677/1‑150 (MQ=255)
                      ctTTCGCTTTAGGCAGAGTCTTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTaaaa                                                                                                                          <  1:566214/149‑1 (MQ=255)
                        ttCGCTTTAGGCAGAGTCTTCTTGATGATTTCATGACGAACCGGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAAcc                                                                                                                        <  1:27123/149‑1 (MQ=255)
                           gCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAAc                                                                                                                   <  2:214649/151‑1 (MQ=255)
                              ttAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCt                                                                                                                 <  1:147669/150‑1 (MQ=255)
                                aGGCAGAGTCGTCTGGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCtgt                                                                                                               <  1:164829/150‑1 (MQ=255)
                                 ggCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGa                                                                                                             <  2:593369/151‑1 (MQ=255)
                                    agagTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATcc                                                                                                          >  2:75276/1‑151 (MQ=255)
                                     gagCCGTCTTAATGATTTAATGACGAGCGAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCGCTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCg                                                                                                         <  1:56142/151‑1 (MQ=255)
                                               gatgatTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGt                                                                                                <  2:635597/150‑1 (MQ=255)
                                                 tgatTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTc                                                                                               >  2:447119/1‑149 (MQ=255)
                                                     ttCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGAt                                                                                         >  2:319959/1‑151 (MQ=255)
                                                         tGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCtt                                                                                     >  1:138582/1‑151 (MQ=255)
                                                              aaCCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCg                                                                                 <  1:156619/150‑1 (MQ=255)
                                                                 cAGTGACCCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGTCCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGt                                                                               <  2:89680/149‑1 (MQ=255)
                                                                 cAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTgg                                                                             >  2:573677/1‑151 (MQ=255)
                                                                          tGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTc                                                                    >  1:407072/1‑151 (MQ=255)
                                                                            ccATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCa                                                                  <  2:387519/151‑1 (MQ=255)
                                                                                 ttGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGgcgc                                                             >  1:570750/1‑151 (MQ=255)
                                                                                    cGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATg                                                          <  1:615372/151‑1 (MQ=255)
                                                                                    cGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATg                                                          <  2:477819/151‑1 (MQ=255)
                                                                                     ggAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGc                                                         >  2:45480/1‑151 (MQ=255)
                                                                                           tAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCaga                                                    <  1:534669/150‑1 (MQ=255)
                                                                                            aGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAgaga                                                  <  2:309528/151‑1 (MQ=255)
                                                                                             gCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGAc                                                 >  2:380037/1‑151 (MQ=255)
                                                                                              ccTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCCCATGACCTTATCCTTCTCAGTAAACCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACa                                                <  1:93941/151‑1 (MQ=255)
                                                                                                  gCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTcc                                            <  1:50535/151‑1 (MQ=255)
                                                                                                    gATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCAc                                          >  1:371385/1‑151 (MQ=255)
                                                                                                     aTGGCTGCTGTTGCGGTTTAGTTGACGACCACTCTTGCTATGGCACATGACCTTATCCTTCTCAGTAAATCCTTAACCTGTGAGCCGGTTACTCGTCAGGGATGCATGCCGGTGGCCAGTTTTCCAGGCGCATTCCCAGAGACATTCCACg                                         <  1:124495/151‑1 (MQ=255)
                                                                                                       ggCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACggg                                       <  1:237724/151‑1 (MQ=255)
                                                                                                              gTTGCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCa                                 <  2:37326/150‑1 (MQ=255)
                                                                                                                 gCGGTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCAc                             <  2:122646/151‑1 (MQ=255)
                                                                                                                   ggTTCAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGt                           >  2:520492/1‑151 (MQ=255)
                                                                                                                       cAGTTGACGACCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCttt                       >  2:437043/1‑151 (MQ=255)
                                                                                                                               gacCACTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtctc                 >  1:561794/1‑149 (MQ=255)
                                                                                                                               gacCACTCGTACGATGGCACATGACCTTATCCTTCTCCGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGACGGTGGCCAGTTTTCCAGGCGCATGCCCACAGACAGTCCACGGGAAGCCCGCACGTCTTTAAtctc                 >  1:22114/1‑149 (MQ=255)
                                                                                                                                   aCTCTTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTa             >  1:675880/1‑149 (MQ=255)
                                                                                                                                      ctTACGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga          <  2:574399/149‑1 (MQ=255)
                                                                                                                                         acgacgGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAttt     >  1:610942/1‑151 (MQ=255)
                                                                                                                                         aCGATGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAttt     >  1:315772/1‑151 (MQ=255)
                                                                                                                                            aTGGCACATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt  <  1:226656/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCTCAACTACGCGGCGCAGCTCTTTCGCTTTAGGCAGAGTCGTCTTGATGATTTCATGACGAACCAGTGAACCTGCCATATTGCGGAACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTT  >  NC_000913/3439851‑3440141

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑