Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,174,322 T→A M202L (ATG→TTG)  U712_10875 ← CCA‑adding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,174,3220TA97.8% 138.8 / NA 46M202L (ATG→TTG) U712_10875CCA‑adding enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (21/24);  minor base C (0/1);  total (21/25)
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

GTTAAAGAGTAAAACGGAAACTCACTCGCCGCAAGCAAGTTCTCTCGTTTATGATAAAAACCAGGGAGTTCCTCGTACAGCCCCGTTTGAATGAGCGTTTGAAGCGCTTGGCGGGAGGCTCTTCCCTGCAGCATTTTCTCAAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACAGCTCT  >  CP006881/2174189‑2174448
                                                                                                                                     |                                                                                                                              
gTTAAAGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTtc                                                                                                                            >  5:401310/1‑138 (MQ=255)
   aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                                         <  8:88562/138‑1 (MQ=255)
   aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                                         >  7:232617/1‑138 (MQ=255)
   aaaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCa                                                                                                                         <  6:318709/138‑1 (MQ=255)
    aaGAGAAGAATGGAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCaaa                                                                                                                       <  5:114376/139‑1 (MQ=255)
         aagaatggAATCTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                                  >  7:126220/7‑139 (MQ=255)
         aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAACCTCTa                                                                                                                  >  5:97589/7‑139 (MQ=255)
         aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                                  >  1:227551/7‑139 (MQ=255)
         aagaatggAAACTCACTCGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTa                                                                                                                  >  7:350820/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                                                                                                     >  3:293236/1‑139 (MQ=255)
                      cactcGTTGAAATCAGGTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTc                                                                                                     >  2:432906/1‑139 (MQ=255)
                                aatcaggTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACa                                                                                           <  2:66109/133‑1 (MQ=255)
                                aatcaggTTTTCTCGTTTATGATAAAAACCAGGAAGCTCTTCATACAATCTCGTTTGAATGAGTGTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACa                                                                                           <  7:427466/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                                                                <  8:368223/102‑3 (MQ=255)
                                                                                               tgTTTGAAGCGCCTGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCgca                                                                >  7:368223/2‑101 (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)
                                                                                                         gcctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCTGTGAAAGCGTAAAGACAAGCAGGCTc                   >  8:206901/4‑138 (MQ=255)
                                                                                                         gcctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCa                  >  1:60865/4‑139 (MQ=255)
                                                                                                          cctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAg                          <  4:41874/128‑1 (MQ=255)
                                                                                                          cctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAg                          >  3:41874/3‑130 (MQ=255)
                                                                                                           ctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATa                <  6:153100/138‑1 (MQ=255)
                                                                                                           ctGGCGAGAAGCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATa                <  8:126220/138‑1 (MQ=255)
                                                                                                                  gaatctctTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTa         <  5:228571/135‑1 (MQ=255)
                                                                                                                  gaagCTCTTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTa         <  8:232617/136‑1 (MQ=255)
                                                                                                                        ctTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCt                        >  3:322166/1‑118 (MQ=255)
                                                                                                                        ctTCCCTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCt                        <  4:322166/118‑1 (MQ=255)
                                                                                                                           cccTGCAGCACTTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGctc   <  4:440333/136‑1 (MQ=255)
                                                                                                                           cccTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGctct  <  1:355869/137‑1 (MQ=255)
                                                                                                                           cccTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGctct  <  1:425245/137‑1 (MQ=255)
                                                                                                                           cccTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGctct  <  1:187081/137‑1 (MQ=255)
                                                                                                                           cccTGCAGCAATTTCTCAAACTCTATAGTTTTTCGTTCGACTGAAACATGGGAGAGGAGTGATTTTTCTTTCGCAATGGCTTCTTCTGTTTCCGGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACGGctct  <  8:249340/137‑1 (MQ=255)
                                                                                                                                     |                                                                                                                              
GTTAAAGAGTAAAACGGAAACTCACTCGCCGCAAGCAAGTTCTCTCGTTTATGATAAAAACCAGGGAGTTCCTCGTACAGCCCCGTTTGAATGAGCGTTTGAAGCGCTTGGCGGGAGGCTCTTCCCTGCAGCATTTTCTCAAACTCTATTGTTTTTCGTTCAACTGAAACATGGGATAGGAGTGATTTTTCTTTCGTGATGGCTTCTTTTGTTTCTAGTGAAAGCGTAAAGCCAAGCTGGCTCATAAAGCGTACAGCTCT  >  CP006881/2174189‑2174448

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