Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,326,925 (C)6→7 coding (817/1185 nt) atoB → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,9191.C95.9% 169.5 / 1.2 49coding (811/1185 nt)atoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (21/26);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACA  >  NC_000913/2326820‑2327013
                                                                                                     |                                                                                                
aaaGCAGGAACAGACACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTa                                                                                                   >  1:1378255/1‑100 (MQ=255)
  aGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cc                                                                                               >  1:500549/1‑100 (MQ=255)
   gCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑ccc                                                                                              >  2:2217612/1‑100 (MQ=255)
      ggAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTagatcgg                                                                                             >  1:2101085/1‑94 (MQ=255)
      ggAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑cccccc                                                                                          >  1:1472577/1‑100 (MQ=255)
       gAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑ccccccc                                                                                         <  2:2173455/100‑2 (MQ=255)
       gAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑ccccccc                                                                                         <  1:487046/100‑2 (MQ=255)
         aCAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑ccccccctg                                                                                       <  2:1476092/100‑4 (MQ=255)
            gTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑ccccccctggc                                                                                     <  1:467600/99‑6 (MQ=255)
                ccGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCa                                                                                 >  2:288758/1‑100 (MQ=255)
                ccGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCa                                                                                 <  2:2150150/100‑1 (MQ=255)
                 cGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCAt                                                                                <  2:828999/100‑1 (MQ=255)
                      ggAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTaaaa                                                                           <  2:1827895/100‑1 (MQ=255)
                       gAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAg                                                                          >  1:1507392/1‑100 (MQ=255)
                        aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGt                                                                         >  2:1095252/1‑100 (MQ=255)
                          cgcgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTa                                                                       <  1:116200/100‑1 (MQ=255)
                            cgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATg                                                                     <  2:1852529/100‑1 (MQ=255)
                            cgTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATg                                                                     <  1:928874/100‑1 (MQ=255)
                              tCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGcc                                                                   <  1:840450/100‑1 (MQ=255)
                               cTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCa                                                                  >  2:760389/1‑100 (MQ=255)
                                 ggTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGc                                                                <  1:1175965/100‑1 (MQ=255)
                                  gTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCg                                                               <  2:1637123/100‑1 (MQ=255)
                                   tATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCgg                                                              <  2:703423/100‑1 (MQ=255)
                                            cgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGcc                                                     >  2:2194343/1‑100 (MQ=255)
                                            cgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGcc                                                     >  1:339386/1‑100 (MQ=255)
                                               tgctgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGccccc                                                  <  1:1302546/100‑1 (MQ=255)
                                                  tgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGc                                               >  1:162570/1‑100 (MQ=255)
                                                    ccGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAt                                             <  1:127849/100‑1 (MQ=255)
                                                     cGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCAtt                                            <  2:1650461/100‑1 (MQ=255)
                                                       ctctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGa                                          >  2:762595/1‑100 (MQ=255)
                                                          tGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATgg                                       >  1:1745206/1‑100 (MQ=255)
                                                           ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                      <  1:100525/100‑1 (MQ=255)
                                                           ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                      >  1:1417005/1‑100 (MQ=255)
                                                           ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                      <  2:782802/100‑1 (MQ=255)
                                                           ggTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCACCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATggg                                      <  2:917547/100‑1 (MQ=255)
                                                              gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTAt                                   >  1:1961604/1‑100 (MQ=255)
                                                              gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTAt                                   >  2:898250/1‑100 (MQ=255)
                                                                 tATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATggg                                <  2:1378257/100‑1 (MQ=255)
                                                                    ggAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGcc                             <  2:1726932/100‑1 (MQ=255)
                                                                    ggAAGAATCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGcc                             <  2:415622/100‑1 (MQ=255)
                                                                     gaagaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCa                            >  1:1587621/1‑100 (MQ=255)
                                                                     gaagaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCa                            >  1:1654434/1‑100 (MQ=255)
                                                                      aagaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAg                           <  1:713575/100‑1 (MQ=255)
                                                                       agaaTCTGCGGCGCTGGCAGCAGGCCTTACCCCCCCCTGGCTCGCATTAAAAGTTCTGCCAGCGGTGGCGTGCCCCCCGCATTGAGGGGTATGGGGCCAg                           >  1:1887001/1‑100 (MQ=255)
                                                                          aTCTGTGGCGCTGGCAGCAGGCCTTGC‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAcc                       <  2:2045723/100‑1 (MQ=255)
                                                                                        cagcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg         <  2:1100881/100‑1 (MQ=255)
                                                                                        cagcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg         >  1:1827389/1‑100 (MQ=255)
                                                                                          gcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg       >  2:1969648/1‑100 (MQ=255)
                                                                                          gcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg       >  2:702068/1‑100 (MQ=255)
                                                                                          gcagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCg       >  1:129930/1‑100 (MQ=255)
                                                                                           cagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGt      <  2:1012773/100‑1 (MQ=255)
                                                                                           cagGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGt      <  2:801231/100‑1 (MQ=255)
                                                                                            agGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGtt     <  1:194404/100‑1 (MQ=255)
                                                                                            agGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGtt     <  1:1791585/100‑1 (MQ=255)
                                                                                            agGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGtt     <  1:2217598/100‑1 (MQ=255)
                                                                                            agGCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGAAAAGCAAAAAGCGtt     <  1:940365/100‑1 (MQ=255)
                                                                                             ggCCTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTa    >  2:1333402/1‑100 (MQ=255)
                                                                                               ccTTACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACa  >  2:683543/1‑100 (MQ=255)
                                                                                                     |                                                                                                
AAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTTA‑‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACA  >  NC_000913/2326820‑2327013

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: