Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,096,479 C→A 16.9% L42F (TTG→TTT)  U712_10420 ← putative coenzyme A transferase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,096,4790CA16.9% 76.8 / 2.2 38L42F (TTG→TTT) U712_10420putative coenzyme A transferase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (22/9);  new base A (2/5);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 7.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGCTACTGGATTGGAACCGATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCATTTAATGT  >  CP006881/2096344‑2096611
                                                                                                                                       |                                                                                                                                    
tCCGGCTACTGGATTGGAACCGATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCaa                                                                                                                                    >  5:295674/1‑138 (MQ=255)
tCCGGCTACTGGATTGGAACCGATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCaa                                                                                                                                    >  4:158395/1‑138 (MQ=255)
              tGGAACCGATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTtcaatca                                                                                                                      <  1:166414/138‑1 (MQ=255)
                     gATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAag                                                                                                                >  4:299915/1‑137 (MQ=255)
                      atatGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAaggag                                                                                                             >  1:36529/1‑139 (MQ=255)
                      atatGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCAATGCAAATCACACTCAAAACCCTTACACCGCAGTCCAAAATCGCTTCTATCAAAGAaggag                                                                                                             >  3:213169/1‑139 (MQ=255)
                           cGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGAc                                                                                                        >  1:194307/1‑139 (MQ=255)
                                cgttCAGGGTTTTGACCCGTTGATTAACAAAACGCGGGCCGATTCCGATATCCGGGAACCGGGCGTCATTGCAAATCACAGTAAAATCCGTAAAACCGCTGGCAAATAGCGCTTCAATCATTGAAGGAGGCGACCCGAc                                                                                                   <  3:260019/136‑1 (MQ=255)
                                cgagccggggTTTGACCCGTTGATTACCAAAAACCGGGCCGATTCCGATATCGGGCAACCCGGCGGCATTGCAAATCCCAGTCAAACCCGTACCACGGCTGTCAAAAATCGCTTCAATCATTCACGCCGGCGACCCGAc                                                                                                   <  4:181946/133‑1 (MQ=255)
                                                ccGTTTATTAACAATCAGCTGTCCGATTCAGATCTCCGGTAAACCGTCGTCATTTCAAATCACATTCAAATCCTTTACCCCGTTGTCAAATAGCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCa                                                                                   <  6:192965/139‑1 (MQ=255)
                                                ccGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCa                                                                                   <  3:158395/139‑1 (MQ=255)
                                                 cGTTGATTACCAATAAGCGGGCCGATTCAGATATCCGGGAAACCGGCGTCTTTGCAAATCACAGTCATATCCGATAACCCGCGGTCCAATATCGCGTCAATCAATGAAGGAGGCGCCCCGACTCCCCCAAACCCGCCa                                                                                   >  4:41061/1‑138 (MQ=255)
                                                 cGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCa                                                                                   >  7:46979/1‑138 (MQ=255)
                                                   ttGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCaaa                                                                                 >  7:12872/1‑138 (MQ=255)
                                                              tAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCCATATAGCTTCAAACAATCCAGGAGGCGACCCGACTACCCCAAAACCCCCAACCATCAGAAccc                                                                      >  1:131514/1‑137 (MQ=255)
                                                                     gCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATGAAAGTCAAAACCGTTACACCGCTGTCCAATATCGCTTCAATCATTGAACGAGGAGACCCGACTCCCCCAAAACCGCCAAACATCATAACCGaacaatcc                                                              >  3:170254/1‑139 (MQ=255)
                                                                     gCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCt                                                                                                                             <  2:123173/76‑1 (MQ=255)
                                                                     gCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCt                                                                                                                             >  1:123173/1‑76 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTtca                                                                                                                          <  2:82649/78‑1 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTtca                                                                                                                          >  1:82649/1‑78 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATTAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGatacaaacc                                                             >  8:192038/1‑135 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCa                                                                                   <  7:282882/117‑1 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCa                                                                                   >  8:282882/1‑117 (MQ=255)
                                                                      ccGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCAAAGTCAAATACGTTACACCGCAGTCCAATATATCTTCAATCAAAGAAGGAGGAGACCCGCCTACCCCAAAACCCCCAAACAAAAGAACCCAACCATccc                                                             >  6:43467/1‑139 (MQ=255)
                                                                            cctatatCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCGTTAACACCGGTGTCAAATATCCCTTCAAACAATCAAGGCGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACa                                                        <  3:299915/135‑1 (MQ=255)
                                                                               atatCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTg                                                    >  5:152320/1‑139 (MQ=255)
                                                                               atatCCGGCAACCCGGCGTCATTGAAAATCACAGTCAAATCCGTTACACCGCTGTCAAATCTCGCTACAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTg                                                    <  7:192038/139‑1 (MQ=255)
                                                                                  tcggggAAACCGGCTGCATTGAAAATCACAGTCAAATCCTTTACACCGCTGTCCAATATCGCTTCAAACAATTACGTAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCaa                                                 <  8:46979/136‑1 (MQ=255)
                                                                                      ggAAACCGGGCTCATCGCAAATCCAAGTCAAATCCTTCACGCCGCTTGCCAATCTCGCTGCATCCATCGACGGCGGCGACCCGACTCCCCCAAAACCGCCCAACCTCAGACCCGATCCATCCCGTACATCTGCAATTGc                                             <  4:330725/139‑1 (MQ=255)
                                                                                                               aGTCAAATCCGTTACACCGCTGTCAAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAACACCGCCAAACATCAGAACCGATCCATCCAGAACATCTGCAATTGCTGTGTCAATGCTTATTGATTTTTg                     >  1:234315/1‑138 (MQ=255)
                                                                                                                 tCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGCACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGaaa                  >  1:360829/1‑139 (MQ=255)
                                                                                                                      tCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGTAGGCTACCCGACTCCCCCAAAACAGCCAAACATCAGAACCGATCCATCCAGAAAATCTGCAATTTCTGTGTCAATGCTGATTGCTTTTTGAAATGgcg             >  3:281168/1‑139 (MQ=255)
                                                                                                                           tACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCAtt        >  3:10956/1‑139 (MQ=255)
                                                                                                                           tACACCGCTGTCAAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCAtt        >  7:115949/1‑139 (MQ=255)
                                                                                                                             cacCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCAttt       <  6:152320/138‑1 (MQ=255)
                                                                                                                             cacCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCAttt       <  2:36529/138‑1 (MQ=255)
                                                                                                                                 gCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCATTTAATg   >  5:76605/1‑138 (MQ=255)
                                                                                                                                  cTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCATTTAATGt  >  7:77232/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                    
TCCGGCTACTGGATTGGAACCGATATGCGAGGCGATCAGGGTTTTGACCCGTTGATTAACAATAAGCGGGCCGATTCCGATATCCGGGAAACCGGCGTCATTGCAAATCACAGTCAAATCCGTTACACCGCTGTCCAATATCGCTTCAATCAATGAAGGAGGCGACCCGACTCCCCCAAAACCGCCAAACATCAGAACCGATCCATCCCGAACATCTGCAATTGCTGTGTCAATGCTGATTGCTTTTTGAAATGGCGCCATTTAATGT  >  CP006881/2096344‑2096611

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