Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,641 T→A 12.5% Q351L (CAG→CTG)  PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,6410TA12.5% 87.3 / 4.3 32Q351L (CAG→CTG) PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (12/16);  new base A (2/2);  total (14/18)
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.23e-01

AACGTCCCTGGCTGGGCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTACATCGTTCATGGCA  >  NC_002947/3189494‑3189787
                                                                                                                                                    |                                                                                                                                                  
aaCGTCCCTGGCTGGGCGAACGGTTGGCGGGCGGCATTGTCGGGGATCATGTCGATCCCGCCCATCAGCCCGCCGCTTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGtgctgc                                                                                                                                                  <  1:192273/150‑1 (MQ=255)
    tCCCTGGCTGGGCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGATCGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACAt                                                                                                                                              >  2:616364/1‑150 (MQ=255)
       ctggctggGCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCATGTGCAGCACCTGGc                                                                                                                                           >  2:721768/1‑150 (MQ=255)
             gggCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCATGTGCAGCACCTGGcccccgt                                                                                                                                     <  1:721768/150‑6 (MQ=255)
             gggCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCATGTGCAGCACCTGGcc                                                                                                                                          <  2:769245/145‑1 (MQ=255)
             gggCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCATGTGCAGCACCTGGcc                                                                                                                                          >  1:769245/1‑145 (MQ=255)
               gCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGa                                                                                                                                   <  1:574085/150‑1 (MQ=255)
                    cGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATcgcg                                                                                                                              >  2:591439/1‑150 (MQ=255)
                                   cTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGcc                                                                                                                                          >  1:474724/1‑123 (MQ=255)
                                   cTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGcc                                                                                                                                          <  2:474724/123‑1 (MQ=255)
                                            aaCCAGTTCGATCCCGCCCATCAGCCCGACGCCCGAGAAGCATAAAGACGACAACTGGGGTTCTAACGCGCGTGGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTc                                                                                                      <  1:616364/150‑1 (MQ=255)
                                                aGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAgc                                                                                                  >  1:687832/1‑150 (MQ=255)
                                                                   gCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACgg                                                                               <  2:413814/150‑1 (MQ=255)
                                                                     ccGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACgggg                                                                             >  1:87191/1‑150 (MQ=255)
                                                                      cGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGt                                                                            <  2:670638/150‑1 (MQ=255)
                                                                                 aCGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAg                                                                 <  2:87191/150‑1 (MQ=255)
                                                                                          gACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGGGGCCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGcc                                                         >  1:816783/1‑150 (MQ=255)
                                                                                              aaCGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCNGCGCAANCGCCGTGGCGCCNGGTTGGCCCGTGCAGGTCAGGCCGTgg                                                     <  1:451202/150‑1 (MQ=255)
                                                                                              aaCGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTgg                                                     <  1:591439/150‑1 (MQ=255)
                                                                                                     ggtCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTCTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATg                                              >  1:649000/1‑150 (MQ=255)
                                                                                                       tCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGcc                                            >  2:367788/1‑150 (MQ=255)
                                                                                                         ttCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGcccc                                          >  1:76627/1‑150 (MQ=255)
                                                                                                                  gcgcAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGctg                                 >  2:716009/1‑150 (MQ=255)
                                                                                                                    gcAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGctggc                               <  2:795729/150‑1 (MQ=255)
                                                                                                                        gTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGctggctctg                           <  1:716009/150‑1 (MQ=255)
                                                                                                                         tGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCTACCAGATCGCGCTCTGCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGctggctctgg                          <  2:363196/150‑1 (MQ=255)
                                                                                                                         tGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGctggctctgg                          <  1:419389/150‑1 (MQ=255)
                                                                                                                                   acacGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTa                >  2:130705/1‑150 (MQ=255)
                                                                                                                                            ccAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTg                                                                           <  2:421722/82‑1 (MQ=255)
                                                                                                                                            ccAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTg                                                                           >  1:421722/1‑82 (MQ=255)
                                                                                                                                            ccAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTACATCGTTCa       <  1:686656/150‑1 (MQ=255)
                                                                                                                                              aGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTACATCGTTCATg     <  2:76627/150‑1 (MQ=255)
                                                                                                                                              aGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTACATCGTTCATg     <  1:130705/150‑1 (MQ=255)
                                                                                                                                                 tgctgcACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCCGGCGCTGGCTCTGGCTCACCAGTACATCGTTCATGGCa  >  1:672790/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACGTCCCTGGCTGGGCGAACGGTTGGCGGGTGGCCTTGTCGGCAACCAGTTCGATCCCGCCCATCAGCCCGACGCCTCTGACGTTGCCGACCAACGGGTGGTCTTCGAACGCGCGCAGGTGGCGCTGGAACACGGGGGCCAGGTGCTGCACATGGCCGACCAGATCGCGCTCTTCGATGATGCGGATGTTTTCCAGCGCAACCGCCGTGGCGACGGGGTGGCCCGTGCAGGTCAGGCCGTGGGCGAATGCCCCCAGGCGCTGGCTCTGGCTCACCAGTACATCGTTCATGGCA  >  NC_002947/3189494‑3189787

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