Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,497,877 A→T 8.0% Y53F (TAC→TTC)  qmcA → Protein QmcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,497,8770AT 8.0% 155.8 / 4.7 50Y53F (TAC→TTC) qmcAProtein QmcA
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base T (2/2);  total (25/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 = 6.56e-01

ATCGCCCTGTTCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGCAG  >  ACB143/5497747‑5498010
                                                                                                                                  |                                                                                                                                     
aTCGCCCTGTTCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGTACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATgg                                                                                                                                  <  5:211006/136‑1 (MQ=255)
aTCGCCCTGTTCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATgg                                                                                                                                  <  5:152982/136‑1 (MQ=255)
    cccTGTTCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGt                                                                                                                              <  3:33223/136‑1 (MQ=255)
         ttCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCg                                                                                                                         >  3:179116/1‑136 (MQ=255)
                 gATCACCGGGGTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGc                                                                                                                 <  8:57260/136‑1 (MQ=255)
                   tCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCt                                                                                                               <  7:17895/136‑1 (MQ=255)
                       cGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCCTCCCCTACATGGACGTGGTCGCCTACCGCCTGCCg                                                                                                           <  3:233052/136‑1 (MQ=255)
                          gtTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGAc                                                                                                         >  7:142334/1‑135 (MQ=255)
                          gtTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGACTGCCGAc                                                                                                         >  6:29115/1‑135 (MQ=255)
                              aaGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAg                                                                                                     >  4:85302/1‑135 (MQ=255)
                               aGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGa                                                                                                   <  5:228200/136‑1 (MQ=255)
                                   ggTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGAcat                                                                                                >  5:141054/1‑135 (MQ=255)
                                     tGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGAcatcat                                                                                             <  6:115813/136‑1 (MQ=255)
                                     tGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGAcatcat                                                                                             <  2:164710/136‑1 (MQ=255)
                                        gcATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGCACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATcctcct                                                                                          >  4:137734/1‑136 (MQ=255)
                                          aTCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTc                                                                                         >  1:223649/1‑135 (MQ=255)
                                            cGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGAt                                                                                      <  8:147566/136‑1 (MQ=255)
                                               gCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGAtgtg                                                                                   >  4:34212/1‑136 (MQ=255)
                                                         gaggagTGGATCGTCGAACGCCTGGGTCGCTACCACAGCCCCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTAACGCCTGCCGACCAACGACATCAGACTCGATGTGCAGCAGCAgg                                                                         >  6:28040/1‑136 (MQ=255)
                                                            gagTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGaaa                                                                      <  6:141054/136‑1 (MQ=255)
                                                             agTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAat                                                                     <  3:42659/136‑1 (MQ=255)
                                                                    cgtcgAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCAc                                                               >  7:132175/1‑135 (MQ=255)
                                                                    cgtcgAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCAGACATGGACGGGGTCGCCTACCGCCGGCCGACCATGGGCATCATCATCGATTTGCCGCCGCAGGAAATCATcccc                                                              >  5:59682/1‑136 (MQ=255)
                                                                         aaCGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGa                                                         <  3:40886/136‑1 (MQ=255)
                                                                         aaCGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGa                                                         <  7:47007/136‑1 (MQ=255)
                                                                              cTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAAcg                                                    >  6:133169/1‑136 (MQ=255)
                                                                                   tCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGTTCGCCTACCGCCTGCCGACCAAGGa                                                                                                   >  2:35477/1‑84 (MQ=255)
                                                                                   tCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGa                                                                                                   <  1:35477/84‑1 (MQ=255)
                                                                                   tCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTa                                               <  3:140296/136‑1 (MQ=255)
                                                                                                aCCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACg                                  >  5:287022/1‑136 (MQ=255)
                                                                                                  ccTCAAGCCGGGCCTGAACATTGTCATCCCCTTCATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGcc                                >  1:135808/1‑136 (MQ=255)
                                                                                                  ccTCAAGCCGGGCCTGAACATTGTCATCCCCTTCATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGcc                                >  1:139917/1‑136 (MQ=255)
                                                                                                  ccTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGcc                                >  6:157241/1‑136 (MQ=255)
                                                                                                    tCAAGCCGGACCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACa                                                       <  8:154813/111‑1 (MQ=255)
                                                                                                    tCAAGCCGGACCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACa                                                       >  7:154813/1‑111 (MQ=255)
                                                                                                      aaGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCtg                             <  3:276591/135‑1 (MQ=255)
                                                                                                          cGGGCCTGCACATTGTCATCCCCTACATGGACGTGGTAGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCATCAGGAAAACATCACCACGGACCACGCGGTAATAGTCGCAAACGCCCTGTGttt                        >  7:273094/1‑136 (MQ=255)
                                                                                                          cGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGttt                        >  2:221891/1‑136 (MQ=255)
                                                                                                               cTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGcc                    >  4:257046/1‑135 (MQ=255)
                                                                                                                 gAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCGGCCGACCAAGGACATCATCCTCGAAGTGCAGCAGCAGGAAATCAACACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTTTGTTTGGCCAAg                 >  6:17578/1‑136 (MQ=255)
                                                                                                                        gTCATCCCCTTCATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTc           <  2:139917/135‑1 (MQ=255)
                                                                                                                        gTCATCCCCTTCATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTc           <  2:135808/135‑1 (MQ=255)
                                                                                                                        gTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTc           <  4:28518/135‑1 (MQ=255)
                                                                                                                          cATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGAc        >  8:192285/1‑136 (MQ=255)
                                                                                                                          cATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGCGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGAc        >  1:47133/1‑136 (MQ=255)
                                                                                                                              cccTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGc    <  6:64500/136‑1 (MQ=255)
                                                                                                                              accTACATGGACGTGGTCGCCTACCGCCTGCCGACAAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGc    <  1:110143/135‑1 (MQ=255)
                                                                                                                               ccTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCg     >  3:201800/1‑134 (MQ=255)
                                                                                                                                cTACCTGGACGTGGTCGCCTACCGCCTGACGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACTACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGCAg  <  5:70406/136‑1 (MQ=255)
                                                                                                                                cTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGCAg  <  3:257046/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
ATCGCCCTGTTCGTCCTGATCACCGTGTTCAAGGGGGTGCGCATCGTGCCCCAGGGCGAGGAGTGGATCGTCGAACGCCTGGGTCGCTACCACAGCACCCTCAAGCCGGGCCTGAACATTGTCATCCCCTACATGGACGTGGTCGCCTACCGCCTGCCGACCAAGGACATCATCCTCGATGTGCAGCAGCAGGAAATCATCACCAAGGACAACGCGGTAATCGTCGCCAACGCCCTGTGTTTCGCCAAGGTAGTCGACCCGCAG  >  ACB143/5497747‑5498010

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: ‑