Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,473,614 C→T 100% W156* (TGG→TGA rpsG ← 30S ribosomal subunit protein S7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,473,6140CT97.9% 117.7 / ‑5.7 40W156* (TGG→TGArpsG30S ribosomal subunit protein S7
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/22);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGT  >  NC_000913/3473475‑3473752
                                                                                                                                           |                                                                                                                                          
cGATGTCACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTCCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCaa                                                                                                                                         <  1:176779/143‑1 (MQ=255)
               gTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATACTCTACGGTAGTGTGCGAAc                                                                                                                          >  1:71385/1‑143 (MQ=255)
                  cGATGGGGGTTGTACGAGCCATTTGTTTCCGCGTTTATCTTTTAGGCGTGCAATTTAAGTAGCCCAAAGCGGGCTGCTGACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGcc                                                                                                                       <  1:405366/143‑1 (MQ=255)
                    aTGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCtt                                                                                                                     <  1:578789/143‑1 (MQ=255)
                               aCGCGCCAGTTTTCTGCTCTTGTATCTCTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAATGGATAATCAACGGTAGTGTGCTAACGCCTTGTTGGCTTCgg                                                                                                          <  2:36053/143‑1 (MQ=255)
                               aCGAGCCATTTGTTTCCTCGGTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCgg                                                                                                          <  2:6634/143‑1 (MQ=255)
                                cGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATATAAGTAGCCTAAAGCGGACTGCTTACTGGAAGCGCCCGCCTGGTGACTAAACCTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCCTGTTGGCTTCGGc                                                                                                         >  2:429771/1‑143 (MQ=255)
                                 gAGCCATTTGTTTCCTCTTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGcc                                                                                                        <  1:131672/143‑1 (MQ=255)
                                    ccATTTGTTTCCTCGTTTATCTTTTAGGCGGTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGGCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGCGTGCGAACGCCTTGTTGGCTTCGGCCATa                                                                                                     <  1:292766/143‑1 (MQ=255)
                                    ccATTTGTTTCCTCGTTTAACTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACAGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGcacta                                                                                                     >  1:138768/1‑139 (MQ=255)
                                       tttGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAGGCGGGCTGCTTACTGGAAGCGACCGCCTGATCACTAAAACTCCGAACGGATCATCCACGGCAGGGTGCGAACGCGCTGTGGGCTTCGGCCATACgg                                                                                                  >  2:200321/1‑143 (MQ=255)
                                        ttGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCGCGCCTGGTGACGAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGt                                                                                                 <  1:363388/143‑1 (MQ=255)
                                           tttCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGaa                                                                                              <  1:238413/143‑1 (MQ=255)
                                            ttCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCGGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAAc                                                                                             <  1:611419/143‑1 (MQ=255)
                                              ccTCGTTTATCTTTTAGGCGTTCAAGTTAAGGAGCCCAAAGCGGGCTGCTGACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGAATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGt                                                                                           <  2:367545/143‑1 (MQ=255)
                                                tCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGGACGTCt                                                                                         >  1:253089/1‑143 (MQ=255)
                                                  tgttATCTTTTAGGCGTTCAATTTAAGTAGCCAAAAGCGGGATGAGTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTc                                                                                       <  1:347348/141‑1 (MQ=255)
                                                      aTCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGGAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGt                                                                                   <  1:513092/143‑1 (MQ=255)
                                                      aTCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGt                                                                                   >  1:156719/1‑143 (MQ=255)
                                                        cTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTCCTGGACGCGCCCGCGTGGTGACTAAAAGTCCGAAGGGATCGCCAACGGGAGTGTGCGCACGCCTTGTTGGCTTCGGCCATACGGTGAGGGTCTTCGCGtct                                                                                 >  1:26608/1‑141 (MQ=255)
                                                                                cccAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGCAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGtgtcc                                                         >  2:438010/1‑139 (MQ=255)
                                                                                  cAAATCGGGCTGCTTACTGGAAGCGCCCGCCTGGCGACTAAAACTCCGAAGGGAGAATCACCGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATTATGTTAACGTCTCCACGTCATCCAACTGCTGTTCCTTTATTCTCTg                                                       >  2:426367/1‑143 (MQ=255)
                                                                                  cAAAGCGGGATGCTTACGGGAAGCGGACGATTGGTCCTTGTATCGCCGAAGGGATAATCAACGGTAGGGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTCTCTGAACTGCAGTGCCTTTGTTTTCTg                                                       >  2:470361/1‑143 (MQ=255)
                                                                                       cGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTgcagcc                                                  >  2:339478/1‑142 (MQ=255)
                                                                                        gggCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCAt                                                 <  1:334607/143‑1 (MQ=255)
                                                                                          gctgctTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGGCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCa                                               >  1:24607/1‑143 (MQ=255)
                                                                                              ctTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGaaa                                           <  1:498111/143‑1 (MQ=255)
                                                                                              ctTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGaaa                                           >  1:20771/1‑143 (MQ=255)
                                                                                                 aCTGGACGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgtt                                        <  1:27018/143‑1 (MQ=255)
                                                                                                  cTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttc                                       <  1:644658/143‑1 (MQ=255)
                                                                                                   tGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTCACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttcg                                      >  1:508279/1‑143 (MQ=255)
                                                                                                          gctcgccTGCTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCa                               <  2:98155/140‑1 (MQ=255)
                                                                                                                  ggTGACTAAAACTCCGAAGGGATAATCGACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCaga                       <  1:251390/143‑1 (MQ=255)
                                                                                                                   gTGACTAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCGTCGGAAAGTTCGTTCGCCAGGCGCagag                      >  1:279882/1‑143 (MQ=255)
                                                                                                                        tAAAACTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGCACGTATTCACGTTTTTTATCTGCATTAGCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATg                 <  1:151901/143‑1 (MQ=255)
                                                                                                                           aaCTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGAt              <  2:226899/143‑1 (MQ=255)
                                                                                                                           aaCTCCGAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCACCAGAAAGTTCGCTCTCCAGGCGCAGAGCCATGGAt              >  2:416239/1‑143 (MQ=255)
                                                                                                                              ttcGAGGGGATAGTCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTa           <  2:348376/141‑1 (MQ=255)
                                                                                                                                 gAAGGGATAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCa        >  2:280201/1‑143 (MQ=255)
                                                                                                                                       aTAATCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGACAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGt  >  2:340740/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                          
CGATGTTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCACCGCGT  >  NC_000913/3473475‑3473752

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