Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,441,060 T→C 14.9% K206E (AAG→GAG)  rpsD ← 30S ribosomal subunit protein S4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,441,0600TC14.9% 54.5 / 4.5 27K206E (AAG→GAG) rpsD30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base C (2/2);  total (13/14)
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.01e-02

ACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGA  >  NC_000913/3440931‑3441202
                                                                                                                                 |                                                                                                                                              
aCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtta                                                                                                                                   <  2:32007/143‑3 (MQ=255)
         aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAgg                                                                                                                          >  2:371575/1‑143 (MQ=255)
          gtcgtaaGCGTCACCTTGGCGTGCGTCGAACTCACTTGTTCGATATCAACCAGGCGCGGTTTTAGAAACGCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt                                                                                                                         <  2:340804/138‑1 (MQ=255)
                         ttGGGGTGCGTCGAACTCGCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCg                                                                                                          <  2:151063/143‑1 (MQ=255)
                            gcgtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcag                                                                                                       <  1:51442/143‑1 (MQ=255)
                                       aCTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGTTCAGGTCTTCTTTAATGTCCGGAGACAGATCAGaa                                                                                            >  2:153251/1‑143 (MQ=255)
                                       aCTCACTTGCGCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGaa                                                                                            <  2:450235/143‑1 (MQ=255)
                                         tCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACg                                                                                          <  1:168720/143‑1 (MQ=255)
                                           aCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATGGTGTCCTCTCTTTGGTACTAATCTTTACTTGGCGCCCAGTTCGAAGAGCAGCTGCTCGTTAATGTACGCCGACAGCACAGAACGCt                                                                                        >  2:213309/1‑143 (MQ=255)
                                             ttGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTcc                                                                                      >  1:376023/1‑143 (MQ=255)
                                                 tCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCGGGCt                                                                                  >  1:218448/1‑143 (MQ=255)
                                                 tCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtt                                                                                                                                    >  1:522695/1‑92 (MQ=255)
                                                 tCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtt                                                                                                                                    <  2:522695/93‑2 (MQ=255)
                                                                   ggTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTCAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCtt                                                                >  2:285326/1‑143 (MQ=255)
                                                                            aaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGc                                                       <  1:27049/143‑1 (MQ=255)
                                                                              ctctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCa                                                     <  2:376023/143‑1 (MQ=255)
                                                                                  gTCACAGCACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCGCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCAt                                                 <  2:196527/143‑1 (MQ=255)
                                                                                            gcTGCATTGTGTCGGCTCTTTGGTAGTAGGCTTTACTTGTGGTAGAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCGCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccag                                       <  2:235298/142‑1 (MQ=255)
                                                                                            ccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCTAGCATCAACCTccag                                       >  2:122279/1‑143 (MQ=255)
                                                                                            ccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccag                                       >  1:490941/1‑143 (MQ=255)
                                                                                             cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagc                                      >  1:497060/1‑143 (MQ=255)
                                                                                              tGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagcc                                     >  1:498466/1‑143 (MQ=255)
                                                                                              tGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCACCTTGCCAGCACCAACTTccagcc                                     >  2:290016/1‑143 (MQ=255)
                                                                                                       tCCTCTCTTTGGTGCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGGTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCt                            <  2:107077/143‑1 (MQ=255)
                                                                                                              tttGGTACTAAGCTTTACGTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACg                     <  1:371575/143‑1 (MQ=255)
                                                                                                                       aaGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTGCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagc            <  2:349016/143‑1 (MQ=255)
                                                                                                                        aGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGCCTCTGCACGCTGCTcagcc           >  2:270693/1‑143 (MQ=255)
                                                                                                                                 tGGAGTAAAGCTCGACGATCAAGAGTTCGGTAATGTCCGCAGACAGATCAGAACGCTCCGACTTACGATTAAACGTACCTTCCCTCTTCCCACCATCAACTTCCCGCCAGGCTGCGTTTCTACACTCCACAGCAAGCTccccg  >  2:69437/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                              
ACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGA  >  NC_000913/3440931‑3441202

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