Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,408,175 C→T L508L (CTG→TTG)  rnr → RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,408,1750CT96.0% 64.6 / NA 25L508L (CTG→TTG) rnrRNase R
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/9);  minor base G (0/1);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAA  >  NC_000913/4407941‑4408416
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
gcgagCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAATGTATCGAACAGACCCAGCGCCCCGACGCGCAAACATTAATTGAAGAGTGCATGATTCTGGAGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGACAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGACGGGCGGTAACa                                                                                                                                                                                                                                    >  1:156374/3‑250 (MQ=255)
                          naaCGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCAAAAATTAATTGAAGAGTGCATGATTCTGGCGCATATATCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGACCCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAgc                                                                                                                                                                                                         >  1:64025/2‑251 (MQ=255)
                                   aCGCCNNATNNNACNNATCGAACAGACCCAGCGTAACGANGCGCACAAATTAATTGAAGCGTGCATGCTNCNGNCGAANNTCTNGGCGGAGCGTNTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCNGNGCACCGAAGCGATTACCTCTTTCCGGTCAGTNCTGGCGGTGCTTGGGCTTTAATTGNCGGGNnggnannnnnnnnn                                                                                                                                                                                                                             >  2:246918/1‑210 (MQ=255)
                                   aCGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCTAATATCTCGGCGGCGCGTTGGGTTTAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGcc                                                                                                                                                                                                                                <  1:246918/219‑1 (MQ=255)
                                                cGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATCTCTCGGCGGCGCGTTTCGTTGAGCAAGCGAAAGAACCGGCAATGTTCCGTATTCACGACAAGCCGAGAACCGAAGCGATTACATCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCg                                                                                                                                                                                   >  1:311408/1‑251 (MQ=255)
                                                cGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGa                                                                                                                                                                                                                             >  1:275343/1‑209 (MQ=255)
                                                               cAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTgg                                                                                                                                                                                                   >  1:87546/1‑220 (MQ=255)
                                                                            cgcACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGACAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGCGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAACCCAtg                                                                                                                                                       >  1:229904/1‑251 (MQ=255)
                                                                             gcACAAATTAATTGAAGAGTGCATGATTCTAGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTACTCACGACAAGCCGAGCACCGAAGCGATTACCTCGTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTTACTACGCGGAGCTGCTGGAGTCGGTTGCCTATCGTCCTGATGCAGAAATGCTGCAAACCAtgc                                                                                                                                                      <  1:294697/251‑1 (MQ=255)
                                                                                      aattGAAGAGTGCATGAGTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCTAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATGACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAAGTGCCGGGCTGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCTGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTgcgc                                                                                                                                              <  1:218593/250‑1 (MQ=255)
                                                                                            agagTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTAGTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGGTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTAcgcn                                                                                                                                                                                                             <  1:64444/181‑2 (MQ=255)
                                                                                                 gCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGGACTGTTCCGTATTCACGACAAGCCGAGAACAGAAGCGATTACCTCTTTCCGTTCCGTGCTGGCGGAGCTGGGGCTGGAATTGcc                                                                                                                                                                                                                                               >  1:315900/1‑142 (MQ=255)
                                                                                                                               gTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTTCTGGCGGAGCTGTGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCTGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCAtgc                                                                                                                                                      <  1:224829/201‑1 (MQ=255)
                                                                                                                                         naaaGCGAAAGAACCGGCACTGTTCCGTATTCACGACACGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGAGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAAAAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCg                                                                                                                                                                                   >  1:64060/2‑162 (MQ=255)
                                                                                                                                                        ggCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCTATTACCTCTTTCCGTTCAGTGCTGGCGTAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAAc                                                                                                                                                                                                                           <  1:220008/107‑1 (MQ=255)
                                                                                                                                                                 ccGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTCCGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGc                                                                                    >  1:1589/1‑233 (MQ=255)
                                                                                                                                                                                            aGCGATTACCGCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCGGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAGACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTTGCCTGGCATTGCAGTCCTATGCGCGCTTTTCTTCGCCGATTCGTCGTTAGCCAGAcc                                       <  1:270489/251‑1 (MQ=255)
                                                                                                                                                                                              cGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTGCCGGGCGGTACCACGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTg                                     >  1:325823/1‑251 (MQ=255)
                                                                                                                                                                                                       ctTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAATTACAGGGCGGTAACAAGCCGGCACCGCGTGACTACGCGGAGCTGCTGGAGACGGTTGCCGAACGTCCTGATGCAGCAATGCTGCAAACCATGCTGATGCGCTCGATGAAACAGGCGATTTACGATCCCGAAAACCGTGGTCACTTTGGcc                                                                                             >  1:290608/1‑186 (MQ=255)
                                                                                                                                                                                                                         cGGAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTACAAACCATGCAGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGcn                                                           <  1:63948/202‑2 (MQ=255)
                                                                                                                                                                                                                           gAGCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGCACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATACTGCTGCGCTCGATGAAACAGGCGATTTACGATACAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATCTctgct        >  1:296534/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTGGGGCTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGc                                                                                         <  1:220399/168‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCTGNAATTGNCNNCCGGTAACAAGCCGGAACCNCGTGACTACGCGGAGNTGCTNGAGTCngg                                                                                                                                                                                            >  2:140099/1‑63 (MQ=255)
                                                                                                                                                                                                                                    cTGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCTCTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGaa  >  1:328658/1‑248 (MQ=255)
                                                                                                                                                                                                                                     tGGAATTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAAcn                                                                                                                                                          <  1:78293/95‑2 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
GAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCACTGTTCCGTATTCACGACAAGCCGAGCACCGAAGCGATTACCTCTTTCCGTTCAGTGCTGGCGGAGCTGGGGCTGGAACTGCCGGGCGGTAACAAGCCGGAACCGCGTGACTACGCGGAGCTGCTGGAGTCGGTTGCCGATCGTCCTGATGCAGAAATGCTGCAAACCATGCTGCTGCGCTCGATGAAACAGGCGATTTACGATCCAGAAAACCGTGGTCACTTTGGCCTGGCATTGCAGTCCTATGCGCACTTTACTTCGCCGATTCGTCGTTATCCAGACCTGACGCTGCACCGCGCCATTAAATATCTGCTGGCGAA  >  NC_000913/4407941‑4408416

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