Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,166,935 2 bp→CG coding (312‑313/483 nt) b3023 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,9350GC99.3% 371.6 / NA 142L105V (CTT→GTT) b3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (80/61);  total (80/62)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,166,9360CG99.3% 363.6 / ‑5.0 141Q104H (CAG→CACb3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (79/61);  minor base A (1/0);  total (80/61)
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.71e-01

Only 100 of 142 total aligned reads displayed.

GCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATC  >  NC_000913/3166867‑3167005
                                                                    ||                                                                     
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 <  1:3488419/76‑1 (MQ=255)
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 >  2:3841171/1‑76 (MQ=255)
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 <  1:2266431/76‑1 (MQ=255)
 ccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCggg                                                                >  2:1852859/1‑76 (MQ=255)
 ccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCggg                                                                <  1:916578/76‑1 (MQ=255)
  cgattGCAGGCAGCCAGGGGCGAATGATGCCCCAGACAGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                               >  1:1207712/4‑76 (MQ=255)
   gcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGcc                                                              <  2:689571/76‑1 (MQ=255)
    cTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGGGGCGGGccc                                                             >  2:704459/1‑75 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            >  1:3130910/1‑76 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            >  2:4451753/1‑76 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            <  1:450067/76‑1 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            >  2:4683987/1‑76 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            >  1:3281839/1‑76 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCCCCCCAACGTGGCGGGCCAc                                                            >  2:1298664/1‑76 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACC‑‑GTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                          >  1:2690404/1‑76 (MQ=255)
      tGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                           >  2:4536614/1‑76 (MQ=255)
      ggcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                           <  2:1834061/75‑1 (MQ=255)
       gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                          >  1:3745027/1‑76 (MQ=255)
       gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                          >  2:2073348/1‑76 (MQ=255)
       gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                          >  1:2869387/1‑76 (MQ=255)
       gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCCCCCACGTGGCGGCCCACgc                                                          >  2:596836/1‑75 (MQ=255)
        caggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGc                                                         >  2:553707/1‑76 (MQ=255)
         aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                        >  1:3502718/1‑76 (MQ=255)
         aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                        <  1:1501865/76‑1 (MQ=255)
         aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                        >  1:4050046/1‑76 (MQ=255)
          ggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCAt                                                       <  2:4050046/76‑1 (MQ=255)
           gcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATa                                                      >  1:544021/1‑76 (MQ=255)
            cagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATac                                                     >  2:3388604/1‑75 (MQ=255)
              gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACg                                                   <  1:4031321/76‑1 (MQ=255)
              gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACg                                                   <  2:1083499/76‑1 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGt                                                  >  1:2796646/1‑76 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGGGGCGGGCCACGGCATAACgc                                                  >  2:226182/1‑75 (MQ=255)
                gaGTGGCGAATGGTGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAGCGTc                                                 <  2:546996/75‑1 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATTTCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                 <  2:111828/76‑1 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                 >  1:482312/1‑76 (MQ=255)
                 aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTcc                                                <  1:2955725/76‑1 (MQ=255)
                   tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTccca                                              >  1:1108332/1‑74 (MQ=255)
                   tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                              >  1:3039153/1‑76 (MQ=255)
                   gggCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                              <  1:4591234/75‑1 (MQ=255)
                    ggCGAATGATGCCGCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcc                                             <  1:3953188/76‑1 (MQ=255)
                    ggCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcc                                             <  1:456068/76‑1 (MQ=255)
                     gCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCg                                            >  2:214128/1‑76 (MQ=255)
                      cGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCgg                                           >  1:2127429/1‑76 (MQ=255)
                      cGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCgg                                           >  1:2549203/1‑76 (MQ=255)
                       gAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCgga                                          >  2:833076/1‑75 (MQ=255)
                       gAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCgga                                          >  2:3722763/1‑75 (MQ=255)
                        aaTGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTa                                         >  1:1760668/1‑76 (MQ=255)
                         atgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTaa                                        >  2:1823351/1‑76 (MQ=255)
                         atgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTaa                                        >  1:2432400/1‑76 (MQ=255)
                          tgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                       <  1:2752958/76‑1 (MQ=255)
                          tgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                       <  1:2900881/76‑1 (MQ=255)
                            atgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGtt                                     >  1:4346369/1‑76 (MQ=255)
                                cccATCCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAca                                 <  2:4447715/76‑2 (MQ=255)
                                cccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAcc                                 <  1:929461/76‑1 (MQ=255)
                                 ccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                >  1:1449464/1‑76 (MQ=255)
                                  cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGCTCACCAt                               >  2:4727103/1‑76 (MQ=255)
                                  cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGCCCACCGGTAAGTACACCAt                               >  2:4102112/1‑76 (MQ=255)
                                   aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAtt                              <  2:3850906/76‑1 (MQ=255)
                                    tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                             <  1:581284/76‑1 (MQ=255)
                                    tACCGTGTGTGAAATATCGTCCAGCTCGCCAACATGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                             >  2:2668010/1‑76 (MQ=255)
                                     aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGc                            >  2:1381593/1‑76 (MQ=255)
                                     aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAACATTGc                            >  2:2262817/1‑76 (MQ=255)
                                     aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGGTCAACAtttc                            >  1:1853342/1‑74 (MQ=255)
                                      ccGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCt                           <  2:1923963/76‑1 (MQ=255)
                                      ccGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCt                           >  2:4351440/1‑76 (MQ=255)
                                      ccGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCt                           <  1:1850565/76‑1 (MQ=255)
                                      ccGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGCAAGTTCCCCATTGcc                           >  2:3363437/1‑75 (MQ=255)
                                         tgtgggAAATCTGGTCCAGCTCGCGAACGTGGCGGTCCGCGGCGTAACGTCCACCGGTAAGCTCACCATTGCTccc                        <  2:3560860/76‑3 (MQ=255)
                                          gtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTcaca                       >  2:1475738/1‑76 (MQ=255)
                                          gggtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTcaca                       <  2:2103053/74‑1 (MQ=255)
                                          gcgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCAGACAGGCCACCGCTAAGTTCACCATTGCTcaca                       >  2:1951939/3‑76 (MQ=255)
                                           tgtgGAATATCTGCCAGCTCTCCAACGTGGCGTGCCACGGCATAACGTCCACCGGTAAGTTCCCCATTGCTCacac                      <  2:1448396/76‑1 (MQ=255)
                                           tgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCacac                      <  2:2962191/76‑1 (MQ=255)
                                            gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACAcc                     >  1:635957/1‑76 (MQ=255)
                                             tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCa                    >  1:3130878/1‑76 (MQ=255)
                                              gAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                   <  1:214128/76‑1 (MQ=255)
                                              gAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGCCCACCGCTAAGGTCACCATTGCTCACAccca                   >  2:913999/1‑74 (MQ=255)
                                              gAAATAGCGCCCAGCCCGCCAACGTGGCGGGCCGCGGCAGAACGTCCACCGGTAAGTTCACCATTGCTCACACCAt                   <  2:1637227/76‑1 (MQ=255)
                                               aaaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATa                  >  2:3054910/1‑76 (MQ=255)
                                                aaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATaa                 <  1:2881375/76‑1 (MQ=255)
                                                 atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                <  2:3226731/76‑1 (MQ=255)
                                                 atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                <  1:2335424/76‑1 (MQ=255)
                                                 atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                <  1:3544151/76‑1 (MQ=255)
                                                   atCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa              <  1:1192668/76‑1 (MQ=255)
                                                    tCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAt             <  1:2632750/76‑1 (MQ=255)
                                                     cctCCACCCCGCCCACCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            <  1:375264/74‑1 (MQ=255)
                                                     cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            <  1:4451753/76‑1 (MQ=255)
                                                     cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            >  1:4222503/1‑76 (MQ=255)
                                                     cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            >  2:3643504/1‑76 (MQ=255)
                                                      gTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa           >  1:4676587/1‑76 (MQ=255)
                                                      gTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa           >  1:2616767/1‑76 (MQ=255)
                                                        ccAGCTCGCCAACGTGGCGGGCCACGGCATAACCGCCACCGGGAAGTTCACCATTGCTCACACCATAACGATTacc         >  1:1551835/1‑74 (MQ=255)
                                                         cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCg        <  1:581127/76‑1 (MQ=255)
                                                          aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg       <  2:4545758/76‑1 (MQ=255)
                                                          aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg       >  1:1558195/1‑76 (MQ=255)
                                                          aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg       <  2:2723278/76‑1 (MQ=255)
                                                            cTCGCCAGGGTGGCGGGCCACGGCATGACGTCCACCGGTAAGTTCACCATTGCTCGCACCATAACGATTATCGGGa     <  2:1723638/76‑1 (MQ=255)
                                                            cTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGCTTATCGGGa     >  2:1242343/1‑76 (MQ=255)
                                                             tCGCCAACGTGGCGGGGCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa    >  1:2256771/1‑76 (MQ=255)
                                                             tCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa    <  1:3623359/76‑1 (MQ=255)
                                                                    ||                                                                     
GCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATC  >  NC_000913/3166867‑3167005

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