Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,304,981 +C intergenic (+588/+127) eco → / ← mqo serine protease inhibitor ecotin/malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,9771.C96.3% 172.6 / ‑0.6 54intergenic (+584/+131)eco/mqoserine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (25/27);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTA  >  NC_000913/2304875‑2305110
                                                                                                       |                                                                                                                                      
cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAACCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCa                                                                                                   >  1:1642918/1‑140 (MQ=2)
 ggATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAAc                                                                                                 >  2:1394547/1‑140 (MQ=32)
   aTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACa                                                                                                >  1:247694/1‑140 (MQ=17)
      aGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCt                                                                                             <  1:926176/140‑1 (MQ=18)
       ggCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGCGCTCCACCCCCGTAGGTTGGAGAAGGCGGCTACGCCGCCTCCAACATcca                                                                                            >  2:910141/1‑138 (MQ=11)
                      tCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGCTACATCGCATCCAACATCTAACGCCCACGCCGGtt                                                                             >  2:1360743/1‑140 (MQ=11)
                                 aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGa                                                                  >  2:101914/1‑140 (MQ=255)
                                   gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACg                                                                <  2:473971/140‑1 (MQ=255)
                                      cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTg                                                             >  2:422375/1‑140 (MQ=255)
                                      cAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTg                                                             >  2:1540672/1‑140 (MQ=255)
                                       aaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                            >  1:534097/1‑140 (MQ=255)
                                         gCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcg                                                          <  1:1285609/140‑1 (MQ=255)
                                         gCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcg                                                          <  1:20768/140‑1 (MQ=255)
                                          ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGgcgc                                                         <  1:422375/140‑1 (MQ=255)
                                           cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                        <  1:117075/140‑1 (MQ=255)
                                            ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACGTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                       <  2:23402/140‑1 (MQ=255)
                                            ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                       >  2:276483/1‑140 (MQ=255)
                                             gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                      >  1:742744/1‑140 (MQ=12)
                                              ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCt                                                     <  2:1546802/140‑1 (MQ=255)
                                               tGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt                                                    >  2:1285609/1‑140 (MQ=255)
                                                gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTa                                                   <  1:209686/140‑1 (MQ=255)
                                                gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  1:1131895/140‑1 (MQ=255)
                                                gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  1:1394547/140‑1 (MQ=37)
                                                 ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  2:1642918/140‑1 (MQ=255)
                                                 ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  2:264816/140‑1 (MQ=255)
                                                 ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  1:775670/140‑1 (MQ=255)
                                                 ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                  <  1:924151/140‑1 (MQ=255)
                                                  cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                 >  2:117075/1‑140 (MQ=255)
                                                    gATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                               <  2:742744/140‑1 (MQ=12)
                                                      tGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGCa                                             >  2:931618/1‑140 (MQ=255)
                                                           cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGc                                        <  1:1258165/140‑1 (MQ=255)
                                                           cGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGc                                        <  1:910141/140‑1 (MQ=255)
                                                            gCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGcc                                       <  1:1639037/140‑1 (MQ=255)
                                                                gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTTCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACAt                                   >  1:1137721/1‑140 (MQ=255)
                                                                gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACAt                                   <  2:534097/140‑1 (MQ=255)
                                                                gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACAt                                   <  2:247694/140‑1 (MQ=255)
                                                                     tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGc                              <  1:101914/140‑1 (MQ=255)
                                                                              gCCTACACCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCGCATACCATGc                     >  2:1390544/1‑140 (MQ=255)
                                                                                  acacCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCgg                               <  2:1449948/126‑1 (MQ=255)
                                                                                  acacCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCgg                               >  1:1449948/1‑126 (MQ=255)
                                                                                   cacCGCTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCCGAt                >  2:1606189/1‑140 (MQ=255)
                                                                                       gttgttAAGTGCGCCGCCCCGGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg            <  1:1360743/138‑1 (MQ=21)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATa                           <  1:112982/124‑1 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATa                           >  2:112982/1‑124 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtg               <  1:132158/136‑1 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtg               >  2:132158/1‑136 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc           >  1:162291/1‑140 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc           >  1:1176929/1‑140 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc           >  2:638012/1‑140 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc           >  1:1047664/1‑140 (MQ=255)
                                                                                        cTGTGAAGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCAGATGTGGc           >  2:906558/1‑140 (MQ=255)
                                                                                             aaGTGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATc      >  1:646781/1‑140 (MQ=255)
                                                                                                tGCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCAtt   >  1:616418/1‑140 (MQ=255)
                                                                                                 gCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCTACATCCGGCATACCATGCCGGATGTGGCGTATCATTa  <  2:646781/140‑1 (MQ=255)
                                                                                                 gCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTa  <  1:1540672/140‑1 (MQ=255)
                                                                                                 gCTCCACCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCAGATGTGGCGTATCATTa  <  1:906558/140‑1 (MQ=255)
                                                                                                       |                                                                                                                                      
CGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTA  >  NC_000913/2304875‑2305110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: