Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,044,722 C→G 13.3% G344R (GGG→CGG) ‡ PP_3566 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,044,7220CG13.3% 129.1 / 9.9 45G344R (GGG→CGG) ‡PP_3566MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base G (3/3);  total (21/24)
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

ACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCC  >  NC_002947/4044578‑4044861
                                                                                                                                                |                                                                                                                                           
acacAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGc                                                                                                                                        >  1:1131659/1‑150 (MQ=255)
  acaACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTg                                                                                                                                      >  1:1062474/1‑150 (MQ=255)
  acaACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTg                                                                                                                                      >  1:423402/1‑150 (MQ=255)
   caACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTgg                                                                                                                                     <  1:1058324/150‑1 (MQ=255)
        aGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCGGGCGGCTgaatgcg                                                                                                                                <  1:628050/150‑7 (MQ=255)
        aGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCGGGCGGCTgaatgcg                                                                                                                                <  1:320523/150‑7 (MQ=255)
            ccACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCGGGCGGCTgaat                                                                                                                                   >  2:1223322/1‑140 (MQ=255)
            ccACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCGGGCGGCTgaat                                                                                                                                   <  1:1223322/143‑4 (MQ=255)
            ccACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCa                                                                                                                            >  2:938079/1‑150 (MQ=255)
                 cccGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCg                                                                                                                       <  2:423402/150‑1 (MQ=255)
                  ccGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGt                                                                                                                      <  2:507716/150‑1 (MQ=255)
                   cGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTa                                                                                                                     <  2:1062474/150‑1 (MQ=255)
                           ccGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCAc                                                                                                             >  1:1239119/1‑150 (MQ=255)
                               ccccACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTg                                                                                                         >  2:341631/1‑150 (MQ=255)
                                             cGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCGGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGt                                                                                           <  1:1362159/150‑1 (MQ=255)
                                               ccTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGNGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTc                                                                                         <  1:894412/150‑1 (MQ=255)
                                                  tGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCGTTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATc                                                                                      <  1:385019/150‑1 (MQ=255)
                                                     cagccaAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGACGGGTCCAGTGCGTACACtt                                                                                                                >  2:1017540/1‑121 (MQ=255)
                                                     cagccaAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGACGGGTCCAGTGCGTACACtt                                                                                                                <  1:1017540/121‑1 (MQ=255)
                                                          aaaCTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCCTCTGCACCGc                                                                              >  1:463042/1‑150 (MQ=255)
                                                             cTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGGGTCAGCCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCaaa                                                                           <  1:1084953/150‑1 (MQ=255)
                                                                 cTGAACAGCGCACTGGCCACCGGCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCaaa                                                                           <  1:714035/146‑1 (MQ=255)
                                                                 cTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCaaa                                                                           >  2:714035/1‑146 (MQ=255)
                                                                  tGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCa                                                                      >  2:250261/1‑150 (MQ=255)
                                                                  tGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAACTccc                                                                      >  2:1148399/1‑149 (MQ=255)
                                                                                         gACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAgg                                               <  2:215002/150‑1 (MQ=255)
                                                                                          aCCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGcc                                                                                                      <  1:190284/94‑1 (MQ=255)
                                                                                          aCCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGcc                                                                                                      >  2:190284/1‑94 (MQ=255)
                                                                                          aCCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCACATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGc                                              <  1:232350/150‑1 (MQ=255)
                                                                                           ccccAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAg                                                      <  1:1378503/141‑1 (MQ=255)
                                                                                           ccccAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAg                                                      >  2:1378503/1‑141 (MQ=255)
                                                                                              cAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCAc                                                                                 >  1:111604/1‑111 (MQ=255)
                                                                                              cAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCAc                                                                                 <  2:111604/111‑1 (MQ=255)
                                                                                              cAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAAt                                          <  1:341631/150‑1 (MQ=255)
                                                                                                    cccGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCACCATGCCGGTAAGGCCCAGGCCAATCACGCt                                    <  1:399541/150‑1 (MQ=255)
                                                                                                             gAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCg                                                           >  1:690001/1‑118 (MQ=255)
                                                                                                             gAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCg                                                           <  2:690001/118‑1 (MQ=255)
                                                                                                                           tGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCt                                    >  2:1291423/1‑127 (MQ=255)
                                                                                                                           tGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCt                                    <  1:1291423/127‑1 (MQ=255)
                                                                                                                            gTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAg            >  1:1325321/1‑150 (MQ=255)
                                                                                                                                  cGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCaa                                           <  1:43051/113‑1 (MQ=255)
                                                                                                                                  cGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCaa                                           >  2:43051/1‑113 (MQ=255)
                                                                                                                                    cATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCCCGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAgg    <  1:1353835/150‑1 (MQ=255)
                                                                                                                                     ggccAGCCGCGGGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGgc   >  2:935234/4‑150 (MQ=255)
                                                                                                                                      gccAGCCGCGGGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGgcc  >  1:47709/3‑150 (MQ=255)
                                                                                                                                                cgcgGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCt                                    >  2:1101518/1‑106 (MQ=255)
                                                                                                                                                cgcgGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCt                                    <  1:1101518/106‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
ACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCACCGCCGACCCCAGCGCCCCGCCGAGGAAGATGCTGGTCATGTACACCGCATTCAGCCGCCCGCGGCTGGCCGGGTCCAGTGCGTACACTTCACGCTGGCCGATGACCATGTTCATCTGCACCGCAAAGTCCAGCAACACGCCGGTAAGGCCCAGGCCAATCACGCTGTAACCCGGCACTGTCAGCCCCAGCAACAAGGCC  >  NC_002947/4044578‑4044861

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