Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,174,296 T→C Q210Q (CAA→CAG)  U712_10875 ← CCA‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,174,2960TC97.4% 113.8 / NA 38Q210Q (CAA→CAG) U712_10875CCA‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/18);  minor base A (1/0);  total (20/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 = 1.00e+00

AAGCGCAGCCCAAAGCTCTTCACGGGATGTTAAAGAGTAAAACGGAAACTCACTCGCCGCAAGCAAGTTCTCTCGTTTATGATAAAAACCAGGGAGTTCCTCGTACAGCCCCGTTTGAATGAGCGTTTGAAGCGCTTGGCGGGAGGCTCTTCCCTGCAGCATTTTCTCAAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCT  >  CP006881/2174161‑2174430
                                                                                                                                       |                                                                                                                                      
aaGAGCAGCCCAAAGCTCCTCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTgg                                                                                                                                     >  7:102138/1‑139 (MQ=255)
  gagcAGCCCAAAGCTCCTCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCg                                                                                                                                   <  7:11422/137‑1 (MQ=255)
  gagcAGCCCAAAGCTCCTCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCg                                                                                                                                   <  6:100129/137‑1 (MQ=255)
  gagcAGCCCAAAGCTCCTCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCg                                                                                                                                   >  5:38443/3‑139 (MQ=255)
    gcAGCCCAAAGCTCCTCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCgag                                                                                                                                 <  7:19512/139‑3 (MQ=255)
                  ctCACGGGATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGTTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTATGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTg                                                                                                                    >  6:275786/2‑138 (MQ=255)
                        ggATGTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTgcagcaat                                                                                                             <  4:241331/139‑3 (MQ=255)
                            gTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTtc                                                                                                          >  5:401310/1‑138 (MQ=255)
                               aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                       <  6:318709/138‑1 (MQ=255)
                               aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                       <  8:88562/138‑1 (MQ=255)
                               aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                       >  7:232617/1‑138 (MQ=255)
                                aaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCaaa                                                                                                     <  5:114376/139‑1 (MQ=255)
                                     aagaatggAATCTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                >  7:126220/7‑139 (MQ=255)
                                     aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAACCTCTa                                                                                                >  5:97589/7‑139 (MQ=255)
                                     aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                >  7:350820/7‑139 (MQ=255)
                                     aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                >  1:227551/7‑139 (MQ=255)
                                       gaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                >  7:383888/5‑137 (MQ=255)
                                          tggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                       <  3:318643/126‑1 (MQ=255)
                                          tggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCATGGCGAGAAGCTCTTACCTGCAGCAATTTCTCa                                                                                                       >  4:318643/2‑127 (MQ=255)
                                                  cactcGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTc                                                                                   >  2:432906/1‑139 (MQ=255)
                                                  cactcGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTc                                                                                   >  3:293236/1‑139 (MQ=255)
                                                            aatcaggTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACa                                                                         <  7:427466/133‑1 (MQ=255)
                                                            aatcaggTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACa                                                                         <  2:66109/133‑1 (MQ=255)
                                                                     ttctcGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGAc                                                                                <  6:337077/122‑1 (MQ=255)
                                                                     ttctcGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGAc                                                                                >  5:337077/2‑123 (MQ=255)
                                                                                         cAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCgcaat                                            <  8:102138/139‑5 (MQ=255)
                                                                                                               cGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGa                       >  6:118927/1‑138 (MQ=255)
                                                                                                                   gacaggAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACGCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCCGTGAAAGCg                  <  5:436969/134‑1 (MQ=255)
                                                                                                                     aaTGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGt                 <  5:118927/138‑1 (MQ=255)
                                                                                                                           tgTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCgca                                              >  7:368223/2‑101 (MQ=255)
                                                                                                                           tgTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCgca                                              <  8:368223/102‑3 (MQ=255)
                                                                                                                           tgTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGcc           >  3:36015/2‑138 (MQ=255)
                                                                                                                             tttGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAg        <  7:33632/139‑1 (MQ=255)
                                                                                                                               tGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCt      >  7:249340/1‑139 (MQ=255)
                                                                                                                               tGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTGTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCt      >  2:48576/1‑139 (MQ=255)
                                                                                                                                gacgcgcCTGGCGAGAACCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGATCGACTGAAACAAGGGAGCGGAGTGATTTTTCTTTCGCAATGGCTGCTTCTGGTTCCGCTGAAAGCGTAACGCCAAGCt      <  8:430035/135‑1 (MQ=255)
                                                                                                                                gAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCt      <  7:34147/138‑1 (MQ=255)
                                                                                                                                   gcgcCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTTTTTACGGTGAAAGCGTAAAGCCAAGcctgtc  >  6:436969/1‑134 (MQ=255)
                                                                                                                                       |                                                                                                                                      
AAGCGCAGCCCAAAGCTCTTCACGGGATGTTAAAGAGTAAAACGGAAACTCACTCGCCGCAAGCAAGTTCTCTCGTTTATGATAAAAACCAGGGAGTTCCTCGTACAGCCCCGTTTGAATGAGCGTTTGAAGCGCTTGGCGGGAGGCTCTTCCCTGCAGCATTTTCTCAAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCT  >  CP006881/2174161‑2174430

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