Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,685,100 A→C 100% pseudogene (1022/1645 nt) ilvG ← ECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,685,1000AC98.7% 255.4 / ‑5.0 73pseudogene (1022/1645 nt)ilvGECK3760:JW3740+JW3741:b4488; large subunit of acetolactate synthase II
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (33/39);  total (33/40)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGATCCGACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAGTGTTGCTGCCAGTCATTG‑‑ATTTAACGGCTGCTGTAATGCT  >  AP009048/3685033‑3685165
                                                                   |                                                                   
cGATCTGACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCAt                                                                   <  1:3564706/70‑1 (MQ=255)
cGATCTGACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATg                                                                  >  1:2930544/1‑71 (MQ=255)
cGATCTGACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATg                                                                  >  1:5188574/1‑71 (MQ=255)
   tctgACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGtt                                                                <  1:1083127/67‑1 (MQ=255)
   tctgACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGtt                                                                >  1:3244739/4‑70 (MQ=255)
   tctgACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGtt                                                                <  1:1663358/67‑1 (MQ=255)
     tgACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcat                                                             >  1:4478509/2‑71 (MQ=255)
       aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGt                                                                 <  1:3005196/65‑1 (MQ=255)
       aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcatc                                                            <  1:4824178/70‑1 (MQ=255)
       aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtca                                                              <  1:4456901/68‑1 (MQ=255)
       aCAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtc                                                               <  1:1742041/67‑1 (MQ=255)
        cAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCAc                                                          <  1:3741825/71‑1 (MQ=255)
         aGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACg                                                         >  1:673425/1‑71 (MQ=255)
          gttgttTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGc                                                        <  1:4712479/71‑1 (MQ=255)
              ttttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCt                                                    <  1:2439088/71‑1 (MQ=255)
               tttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTg                                                   >  1:2093256/1‑71 (MQ=255)
                ttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGt                                                                 <  1:103848/56‑1 (MQ=255)
                ttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgc                                                  >  1:3025038/1‑71 (MQ=255)
                ttAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgc                                                  >  1:3018975/1‑71 (MQ=255)
                  aacaacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcatca                                                           <  1:2745393/60‑1 (MQ=255)
                   acaacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcg                                               <  1:2283083/71‑1 (MQ=255)
                   acaacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcg                                               <  1:4563555/71‑1 (MQ=255)
                     aacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcg                                               >  1:4663244/1‑69 (MQ=255)
                     aacaacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcg                                                 <  1:3779452/67‑1 (MQ=255)
                       caacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgt                                           <  1:147034/71‑1 (MQ=255)
                        aacGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCAt                                                                   >  1:4651101/1‑46 (MQ=255)
                         acGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgtgt                                         <  1:3207033/71‑1 (MQ=255)
                         acGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgtgt                                         <  1:2355092/71‑1 (MQ=255)
                          cGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTtcatca                                                           >  1:2907474/1‑52 (MQ=255)
                          cGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgtgta                                        <  1:4719824/71‑2 (MQ=255)
                            gcgcGTAGATAGCGTCACCGGGATTGTCGTAACGCCATGCATGTTCATCACGCAGCTGCGCGCAg                                            <  1:5012747/65‑1 (MQ=255)
                              gcgTAGATATCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTagctg                                    >  1:4367778/1‑71 (MQ=255)
                              gcgTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcgc                                              <  1:4659979/61‑1 (MQ=255)
                              gcgTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTagctg                                    >  1:694240/1‑71 (MQ=255)
                                gTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTagctgc                                   >  1:4663149/1‑70 (MQ=255)
                                 tAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcg                                               <  1:3219331/57‑1 (MQ=255)
                                 tAGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTagctgc                                   >  1:3723704/1‑69 (MQ=255)
                                  aGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAg                                >  1:1363533/1‑71 (MQ=255)
                                  aGATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAg                                >  1:2098383/1‑71 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCa                                             >  1:1935695/1‑57 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               >  1:3717561/1‑71 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               <  1:3103130/71‑1 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               <  1:1445725/71‑1 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               <  1:927534/71‑1 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               <  1:4553081/71‑1 (MQ=255)
                                   gATAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               >  1:1042082/1‑71 (MQ=255)
                                    aTAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTagctgc                                   <  1:3043912/66‑1 (MQ=39)
                                    aTAGCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGt                               >  1:4836522/1‑70 (MQ=255)
                                    aTAGCGGCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCTGTGTAGCTGCCAGTc                              <  1:4043269/71‑1 (MQ=255)
                                     tAGCGTCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAGTGTTGCTGCCAGTCa                             <  1:440402/71‑1 (MQ=255)
                                       gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCATCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:4656863/71‑1 (MQ=255)
                                       gCGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           <  1:2116878/71‑1 (MQ=255)
                                        cGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCATTCa                             >  1:95048/1‑68 (MQ=39)
                                        cGTCACCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAtt                           >  1:2871414/1‑70 (MQ=255)
                                            aCCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcgc                                              <  1:5159559/47‑1 (MQ=255)
                                            aCCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcgc                                              <  1:1937699/47‑1 (MQ=255)
                                            aCCGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTG‑‑Atat                    >  1:1800685/1‑69 (MQ=255)
                                             ccGGGATGGTCGTAACGCCAGGCATTTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTG‑‑Atat                    >  1:1545280/1‑68 (MQ=255)
                                              cGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCAt                            <  1:518618/63‑1 (MQ=39)
                                              cGGGATGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTG‑‑Atatt                   >  1:3494401/1‑67 (MQ=255)
                                                  aTGGTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAg                                            >  1:1486134/1‑43 (MQ=255)
                                                     gTCGTAACGCCAGGCATGTTCATCACGCAGCTgcgcg                                               <  1:1159586/37‑1 (MQ=255)
                                                     gTCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgc              <  1:5023601/70‑1 (MQ=255)
                                                      tCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctg            >  1:4744227/1‑71 (MQ=255)
                                                      tCGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctg            <  1:689229/71‑1 (MQ=255)
                                                       cGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTg                          >  1:3813373/1‑56 (MQ=38)
                                                       cGTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGgctgc           >  1:3545106/1‑71 (MQ=255)
                                                        gTAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTg                          >  1:4249660/1‑55 (MQ=38)
                                                         tAACGCCAGGCATGTTCATCACGCAGCTGCGCGCAgt                                           <  1:2623550/37‑1 (MQ=255)
                                                           aCGCCAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTa       <  1:4638299/71‑1 (MQ=255)
                                                               cAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGc   <  1:2156418/71‑1 (MQ=255)
                                                               cAGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGc   >  1:1389621/1‑71 (MQ=255)
                                                                aGGCATGTTCATCACGCAGCTGCGCGCAGTGTAGCTGCCAGTCATTGATATTTAACGGCTGCTGTAATGCt  >  1:2117376/1‑71 (MQ=255)
                                                                   |                                                                   
CGATCCGACAGTTGTTTTAACAACAACGGCGCGTAGATAGCGTCACCGGGATGGTCGTAACGCCAGGAATGTTCATCACGCAGCTGCGCGCAGTGTTGCTGCCAGTCATTG‑‑ATTTAACGGCTGCTGTAATGCT  >  AP009048/3685033‑3685165

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