Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,298,393 A→G 16.8% intergenic (+521/+194) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,298,3930AG16.8% 50.1 / 2.2 27intergenic (+521/+194)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (13/9);  new base G (3/2);  total (16/11)
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

ACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTT  >  CP009273/2298161‑2298643
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTCTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCCGTTACGCCGCATCCGACACCCAACGCCCAAGCCGGTTGCCTCATGAGACGCACGCGCGCATTATCAGGCCt                                                                                                                                                                                                                  >  1:203988/1‑275 (MQ=11)
            gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGACGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTCGTCAGGCCTACCCCGCtctg                                                                                                                                                                                                       >  2:684155/1‑274 (MQ=255)
                              ggCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAACGCTGTGAAGTGCTCCACCCCGCAgg                                                                                                                                                                                    >  1:456999/1‑275 (MQ=38)
                                                   acCGCTGTGAAGTGCTCCATCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAATGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGTTGGATAAGGCGGTTAcgcc                                                                                                                                                                >  1:591236/1‑274 (MQ=255)
                                                                    cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAGATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGATGTGAAGTGCTCCACCCCGTAGGTTGTATCAGGCGGTTACGCCCCAGCCAACATCTACCg                                                                                                                                               >  1:640229/1‑274 (MQ=255)
                                                                      ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCACAAGCCGGTTGCCTGATGCGACGCTGGCGAGTCTTATAAGGCATACACCGCTGTGCAGTGCTCCACCCCGTAGTTTGGATAAGGCGATTACGCCGCATCCAACAACTAACGcc                                                                                                                                             >  1:206377/1‑274 (MQ=255)
                                                                      ccccGTAGGTTGGATAAGGAGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGc                                                                                                                                                                                                                                           >  1:590549/1‑180 (MQ=255)
                                                                          gTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                                                                                                                                                       <  2:590536/180‑1 (MQ=39)
                                                                                   aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCATGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATCGCTGTGAAGTGCTCCACCCCGTAGGTTGTATAAGGCGGTTACGCCGCATCCCAc                                                                                                                                                       >  1:707778/1‑251 (MQ=18)
                                                                                             tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGcc                                                                                                                                                                                                                   >  1:594404/1‑181 (MQ=32)
                                                                                             tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTACACCGCGGTGAAGTACTCCACCCCGTCGGCTGGATAAGGCGGTTACGCCGCATCCACCATATAACGCACGAGCCGGGTGCCCGATGctgta                                                                                                                      >  1:225842/1‑270 (MQ=11)
                                                                                                 ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                               <  2:594391/181‑1 (MQ=31)
                                                                                                               cAACGTCCGGGCCGGTTGCCTGATGCGACGCTGGCGCTTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCa                                                                                                   <  1:594816/275‑1 (MQ=25)
                                                                                                                accggccGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACNCCCAAGCCGGCGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCAt                                                                                                  <  1:684169/270‑1 (MQ=19)
                                                                                                                            agagaaagagcgcggcgCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCTACCCCGTAGGTCGGAGAAGGCCGCTACGCCGCATCCGACATCTAACGGCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGTGCTCCACCCCTTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCACGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGCATACGCCACATCCg                                                                                      <  1:308391/264‑1 (MQ=14)
                                                                                                                                                                            gAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAACAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCg                                                                                                                                           <  1:457854/174‑1 (MQ=25)
                                                                                                                                                                                              ggttggtgaaGGCGGTCGCGGCGGTTGCGACAGCGGACGCCCGGGCCGGGTTCCGGCTGTGGTGCTGGTGAGTGTAATCATGGCTACACCGCTGTGAAGGGCACCGCCCCGTAGGTTGGAGAAAGCGGTTACGCCGCATCCAACATCTAACTCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGtt                     <  2:206375/266‑1 (MQ=18)
                                                                                                                                                                                                        ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCa                                                                                                   <  1:646272/186‑1 (MQ=37)
                                                                                                                                                                                                         gCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACg                                                                                                                      >  2:72828/1‑166 (MQ=35)
                                                                                                                                                                                                            gTTACGCCGCATCCGAGATCTAACGCCCGAGCCGGTTGCCTGATACGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATNCGCCACATCCGGCATa                                                                                 >  1:531862/1‑200 (MQ=255)
                                                                                                                                                                                                             ttACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                  <  1:72829/166‑1 (MQ=35)
                                                                                                                                                                                                              caCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCTCATCCGCCATACCATCCCGGATGCTGCGCATCATGACAACGAAACAACCACCAGCTCTTTTTGAACGggtgtcggtgtcttatgcg      >  2:90148/2‑266 (MQ=14)
                                                                                                                                                                                                                  ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGttt  <  1:90149/273‑1 (MQ=255)
                                                                                                                                                                                                                      aTCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGa                                                                                                                                          >  1:55759/1‑133 (MQ=37)
                                                                                                                                                                                                                          gACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCg                                                                                                                                      <  2:55758/133‑1 (MQ=37)
                                                                                                                                                                                                                                  aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGNCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGt    >  1:549233/1‑255 (MQ=255)
                                                                                                                                                                                                                                      ccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAATGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTATTATCAACATACCCcaca                                                                                          >  1:655133/1‑165 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
ACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTT  >  CP009273/2298161‑2298643

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