Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,473,615 C→T 50.2% W156* (TGG→TAG)  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,6150CT50.2% ‑6.5 / 47.6 42W156* (TGG→TAG) rpsG30S ribosomal subunit protein S7
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base T (11/10);  total (21/21)
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.80e-01

TTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCAC  >  NC_000913/3473480‑3473747
                                                                                                                                       |                                                                                                                                    
ttACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTCAAACTCCGAAGGGATAACTAACGGTa                                                                                                                               >  1:761274/1‑143 (MQ=255)
  aCGGTAGCGTGCGATGGGTGTTGTACGGGCCATGTGTGTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAgt                                                                                                                             >  1:186161/1‑143 (MQ=255)
   cGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTCgtg                                                                                                                            >  2:544397/1‑143 (MQ=255)
   cGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAgtg                                                                                                                            >  1:433206/1‑143 (MQ=255)
       aGCGTGCGATGGGGGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACGGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCg                                                                                                                        <  1:298853/143‑1 (MQ=255)
             cGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGcc                                                                                                                  <  1:450330/143‑1 (MQ=255)
             cGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGcc                                                                                                                  <  2:744505/143‑1 (MQ=255)
             cGATGGGTGTGGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGcc                                                                                                                  <  2:745332/143‑1 (MQ=255)
                  ggTGTTGTACCAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGCAGCCCAAAGCGGGCTGCTTACTGGAAGCGCGCGGCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCttgtt                                                                                                             <  2:470107/143‑1 (MQ=255)
                         tACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCg                                                                                                      >  1:770583/1‑143 (MQ=255)
                               ccATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATCTAAGTCGCCCAAAGCGGGCTGCTTACTGGACGCGCCCGCCTGGTGACTAACACTCCGAAGGGCTAACCAACGGTAGTGTCCGAACGCCCTGTCGCCTACGGCCacc                                                                                                >  2:212994/1‑141 (MQ=255)
                                      tttCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGCATAACTAACGCTCGTGTGCGAACGCGTTGGTGGCTTCGGCAAGAAGGTGaa                                                                                         >  2:82979/1‑143 (MQ=255)
                                         ccTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATAAGGTGAACGt                                                                                      >  1:81144/1‑143 (MQ=255)
                                         ccTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTAGGACACACGGAGAACgc                                                                                      >  2:715523/1‑142 (MQ=255)
                                           tCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGATTCAAACGCCGAAGGGTTAAACAACATCAGTGCGCGAACGCCTTGTTCGCTTCGGCCATACGGTGAACGTCt                                                                                    >  2:349031/1‑143 (MQ=255)
                                            cgcttctCTTTTAGGCGTGCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCtt                                                                                   <  2:294330/137‑1 (MQ=255)
                                                 aTCTTTTAGGCGTTCAGTTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGGGACTAAAACTCCGAAGGGATAACTCACGGTAGTGTGCGAACGCCTTGGTGGCTTCGGCCATACGGTGAACGCCTTCACGt                                                                              >  2:116003/1‑143 (MQ=255)
                                                 aTCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGt                                                                              >  2:283374/1‑143 (MQ=255)
                                                    ttttAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTc                                                                           >  1:731036/1‑143 (MQ=255)
                                                     tttAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCt                                                                          <  1:112835/143‑1 (MQ=255)
                                                           cGTTCAATCTACGTAGCCCAACGCGGGCTGCTTACTGGACGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTg                                                                    <  1:596530/143‑1 (MQ=255)
                                                            gTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGGCATACGGTGAACGTCTTCACGGTTCTTAACTGc                                                                   >  1:242369/1‑143 (MQ=255)
                                                              ttaaTTTAAGTTGCCAAGAGCGGGCTGTGTACTGGAAGCGCCCGCCTGGTGACTAAAACGCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAg                                                                 <  1:677876/141‑1 (MQ=255)
                                                                  tttAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTAcc                                                             >  1:494072/1‑143 (MQ=255)
                                                                        tAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGtt                                                       <  1:764765/143‑1 (MQ=255)
                                                                        tAGCCCAAAGCGGGCTGCTGACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGtt                                                       <  2:115503/143‑1 (MQ=255)
                                                                           gccAAAGCGGGCTGCTTACGGGACGCGCCCGCCTGGTGACTAACACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTc                                                    <  2:712019/142‑1 (MQ=255)
                                                                           cccAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGG‑TAACTAACGGTAGTGTTCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCt                                                   >  1:93500/1‑143 (MQ=255)
                                                                            ccAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCt                                                   <  2:302237/143‑1 (MQ=255)
                                                                             cAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGCACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTCCACGTTTCTTCACTGCCGTACCTCTGTTTTCTg                                                  >  2:520606/1‑143 (MQ=255)
                                                                                              tGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAgttcg                                 <  2:540320/143‑1 (MQ=255)
                                                                                                gAAGCGCCCGCCTGGCGAGTACTACTCCGAAGGGATCTCCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGTATCAGAAAgttcgtt                               <  2:665477/143‑1 (MQ=255)
                                                                                                   gcgcCCGCCTGGTGACTAAATCTCCGAGTGGATACCCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGc                            <  2:341230/143‑1 (MQ=255)
                                                                                                    cgccctccCGGTGACTAAAACTCCTAAGGGATAGCCAACGGGAGTGTGCGAACGCCTTGTTGGCTGCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGcc                           <  2:246575/143‑1 (MQ=255)
                                                                                                     gcccgccTGGTGACTAAAGCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCGTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCa                          <  2:562978/143‑1 (MQ=255)
                                                                                                        cgccTGGTGACCAAAACGCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCCGCAGCATCAGAAAGTTCGTTCGCCAGgc                       <  1:207971/143‑1 (MQ=255)
                                                                                                            tGGTGACTAAATCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCag                   <  2:643535/143‑1 (MQ=255)
                                                                                                                  cTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCGGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCAt             <  2:186161/143‑1 (MQ=255)
                                                                                                                  cTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCAt             >  1:14291/1‑143 (MQ=255)
                                                                                                                     aaaCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGa          >  1:341179/1‑143 (MQ=255)
                                                                                                                     aaaCTCCGAAGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGa          >  1:148932/1‑143 (MQ=255)
                                                                                                                             aaGGGATAACTAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCGTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCGAGGCCATGGATTTATCAc  >  1:125714/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                    
TTACGGTAGCGTGCGATGGGTGTTGTACGAGCCATTTGTTTCCTCGTTTATCTTTTAGGCGTTCAATTTAAGTAGCCCAAAGCGGGCTGCTTACTGGAAGCGCCCGCCTGGTGACTAAAACTCCGAAGGGATAACCAACGGTAGTGTGCGAACGCCTTGTTGGCTTCGGCCATACGGTGAACGTCTTCACGTTTCTTAACTGCAGTACCTTTGTTTTCTGCAGCATCAGAAAGTTCGTTCGCCAGGCGCAGAGCCATGGATTTATCAC  >  NC_000913/3473480‑3473747

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