Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,498,534 C→T 16.6% L121L (CTG→CTA PP_5563 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,498,5340CT16.6% 47.6 / 2.9 24L121L (CTG→CTAPP_5563membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base T (2/2);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTGGCAACGATCA  >  NC_002947/3498386‑3498676
                                                                                                                                                    |                                                                                                                                              
ttGCCTTGCACTGCCAATTTTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTTCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCa                                                                                                                                               <  2:25533/150‑1 (MQ=255)
           ggCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCAtt                                                                                                                                    >  2:93939/1‑150 (MQ=255)
            gCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTg                                                                                                                                   >  1:258838/1‑150 (MQ=255)
             ccAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCAtt                                                                                                                                    <  1:456769/148‑1 (MQ=255)
             ccAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCAtt                                                                                                                                    >  2:456769/1‑148 (MQ=255)
              cAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTc                                                                                                                                 <  1:532747/150‑1 (MQ=255)
               aaTTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCg                                                                                                                                <  2:258838/150‑1 (MQ=255)
                              aaNGATCCGGCTTGAAGCACACCCGTTNATCGTTGACCCANAACATCNGNCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGTCGAAAGCGTAGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCa                                                                                                                 >  1:498979/1‑150 (MQ=255)
                                tGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATc                                                                                                               <  2:463381/150‑1 (MQ=255)
                                                   cccGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACAACCAGTAGACCTTTCACTAAAGCGTGGCGCTAGGCAGCTCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCa                                                                                            <  1:91375/150‑1 (MQ=255)
                                                           tCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCa                                                                                                                                      <  1:263585/100‑1 (MQ=255)
                                                           tCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCa                                                                                                                                      >  2:263585/1‑100 (MQ=255)
                                                             gTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAg                                                                                  >  1:68078/1‑150 (MQ=255)
                                                                       aaCATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCCTTCTGGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCa                                                                        <  2:290269/150‑1 (MQ=255)
                                                                             tGCCAACCGGAGCCCGGGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGc                                                                  >  2:728149/1‑150 (MQ=255)
                                                                               ccAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCa                                                                    <  1:507722/146‑1 (MQ=255)
                                                                               ccAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCa                                                                    >  2:507722/1‑146 (MQ=255)
                                                                                 aaCCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTTAGCGCAATCATTGt                                                                                                                                  >  2:31591/1‑82 (MQ=255)
                                                                                 aaCCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTTAGCGCAATCATTGt                                                                                                                                  <  1:31591/82‑1 (MQ=255)
                                                                                                                       tagtCCTTTGCCGAAAGCGTGGCTGTGGTTAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTGGTGGGTCAAACTTTCGACCGGTCCATgcgg                       >  2:594524/1‑150 (MQ=255)
                                                                                                                                   gAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTg            >  1:497379/1‑150 (MQ=255)
                                                                                                                                   gAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTg            >  2:56489/1‑150 (MQ=255)
                                                                                                                                    aaaGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTgg           <  1:93939/150‑1 (MQ=255)
                                                                                                                                            gCTGTGGTTAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTGGCAACGATCa  <  1:594524/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCCGTTCATCGTTGACCCAGAACATCTGCCAACCGGAGCCCGAGATCCTGAGCTTCGGGCACGACCAGTAGTCCTTTGCCGAAAGCGTGGCTGTGGTCAGCGCAATCATTGTCGGAGCAACGATAACCATCAGCACAACTGTCAGTTGCAGCGAGCGTAGCTGGAACCCATGCAGCCAAGAGGTTTTTGGTGGGTCAAACTTTCGACCGGTCCATGCGGCGATCAGGCTGGCAACGATCA  >  NC_002947/3498386‑3498676

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