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,9350GC100.0% 369.1 / NA 139L105V (CTT→GTT) b3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (57/82);  total (57/82)
*NC_0009133,166,9360CG99.3% 373.7 / NA 139Q104H (CAG→CACb3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (56/82);  total (57/82)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

GCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATC  >  NC_000913/3166867‑3167005
                                                                    ||                                                                     
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 >  1:2862217/1‑76 (MQ=255)
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 <  2:126449/76‑1 (MQ=255)
gccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                 >  1:3144941/1‑76 (MQ=255)
gccgcTTGCAGGCAGCCAGTGGCGAATGAGGCCCCACAACCTGTGTGTAATATCGTCCAGCTCCTCAACCTGCCgg                                                                 >  2:696128/1‑76 (MQ=255)
 ccgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCggg                                                                <  2:1398164/76‑1 (MQ=255)
  cgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                               >  2:1759720/1‑76 (MQ=255)
   gcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGcc                                                              >  2:4889714/1‑76 (MQ=255)
    ctggcaggcagCCAGTGGCGAATGATGCCGCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCa                                                             <  2:3991426/73‑1 (MQ=255)
    cTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCa                                                             <  2:3031590/76‑1 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            <  1:362962/76‑1 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAc                                                            <  2:4150650/76‑1 (MQ=255)
     ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGGGCCGGGCCAc                                                            >  2:2061274/1‑76 (MQ=255)
      tGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                           <  2:1863995/76‑1 (MQ=255)
      tGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                           <  1:518595/76‑1 (MQ=255)
      tGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                           >  1:5203892/1‑76 (MQ=255)
       gcaggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACgg                                                          <  1:3357049/76‑1 (MQ=255)
        cgggcggcCAGTGTCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGc                                                         <  1:5403751/74‑1 (MQ=255)
         aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCCACGTGGCGGGCCACCGcc                                                        >  2:1440096/1‑75 (MQ=255)
           gtagCCAGTGGCGAATGAGGCTCCATCCCGTGTGTGAAATATCGTGCAGCTCGCCAACGTGGCGGGCCACGGCATa                                                      <  2:3116158/74‑1 (MQ=255)
           gcagCCAGTCGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATa                                                      >  1:5812056/1‑76 (MQ=255)
            cagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATaa                                                     <  1:5143339/76‑1 (MQ=255)
            cagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATaa                                                     >  2:2725195/1‑76 (MQ=255)
             agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAAc                                                    <  1:3258802/76‑1 (MQ=255)
             agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAAc                                                    <  1:5456243/76‑1 (MQ=255)
             agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAAc                                                    >  1:2326044/1‑76 (MQ=255)
             agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAAc                                                    >  1:4140447/1‑76 (MQ=255)
              gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCCCGGCATAACg                                                   >  2:530482/1‑76 (MQ=255)
              gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACg                                                   <  1:1600674/76‑1 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                               <  1:5397597/63‑1 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                               >  2:5397596/1‑63 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACgc                                                  >  2:3397775/1‑75 (MQ=255)
               ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGt                                                  >  1:3991009/1‑76 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATCACGTc                                                 >  2:5577950/1‑76 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                 >  1:2308010/1‑76 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTc                                                 <  2:3983630/76‑1 (MQ=255)
                cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGAACAACGTc                                                 >  2:5181815/1‑76 (MQ=255)
                 aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTcc                                                <  2:405461/76‑1 (MQ=255)
                 aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTcc                                                >  2:3615756/1‑76 (MQ=255)
                  gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCa                                               >  2:4146870/1‑76 (MQ=255)
                  gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCa                                               >  1:1251962/1‑76 (MQ=255)
                   tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                              >  1:3403466/1‑76 (MQ=255)
                   tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                              <  1:2145343/76‑1 (MQ=255)
                   gggCGAATGATGCCCCATACAGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                              <  1:4013561/75‑1 (MQ=255)
                     gCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCg                                            <  2:3276069/76‑1 (MQ=255)
                      cGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcccg                                           >  1:529020/1‑74 (MQ=255)
                       gAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGt                                          >  2:3985464/1‑76 (MQ=255)
                        aaTGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTa                                         <  1:4678034/76‑1 (MQ=255)
                        aaTGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTa                                         <  2:3189479/76‑1 (MQ=255)
                         atgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGAATAACGCCCCCCgggac                                        >  2:2205785/1‑73 (MQ=255)
                          tgatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                       >  1:5144093/1‑76 (MQ=255)
                           gatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCCGGCCACGGGATGAAGTACACCGGTaatg                                      >  1:605243/1‑74 (MQ=255)
                            atgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGtt                                     <  1:1440096/76‑1 (MQ=255)
                             tgCCCGATACCGTGTTTGGACTATCGTCCAGCACTCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTc                                    <  1:2496357/76‑1 (MQ=255)
                             tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTc                                    <  1:2966020/76‑1 (MQ=255)
                               ccccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAc                                  >  1:2249480/1‑76 (MQ=255)
                               ccccATACCGTGTGTGAAATAGAGTCCAGCGAGCCAACGGGGAGGGACACGGGAAAACATAGACCGGGAAGTTCAc                                  >  1:698645/1‑76 (MQ=255)
                                  cATACCGTTTGTGAAATATCGTGCAGATCGCCAACGGGTCGGTCCACGGCACGACGTCCACCGGGAAGTACAACAt                               <  2:2018538/76‑1 (MQ=255)
                                  cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                               <  1:3399046/76‑1 (MQ=255)
                                  cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                               <  2:1962273/76‑1 (MQ=255)
                                  cATACCGTGTGTGAAATATCGCCCAGCCCCCCAACGTGGCGGGCCACGGCATAACGCCCCCCGGTAAGTTCACCAt                               <  1:4949312/76‑1 (MQ=255)
                                   aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAtt                              >  2:3607976/1‑76 (MQ=255)
                                   aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAtt                              <  1:1194532/76‑1 (MQ=255)
                                    tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                             <  2:4786949/76‑1 (MQ=255)
                                    tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                             <  2:2994915/76‑1 (MQ=255)
                                    tACCGTGTGTGAAAGATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                             <  2:2734090/76‑1 (MQ=255)
                                       cGTGTGTGGAATATCGTCCAGCCCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTc                          >  1:5214020/1‑76 (MQ=255)
                                       cGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTc                          <  1:1844622/76‑1 (MQ=255)
                                       cGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTc                          <  1:3431569/76‑1 (MQ=255)
                                        gtgtgtGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTca                         <  2:2397567/76‑1 (MQ=255)
                                         tgtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTcac                        >  2:549681/1‑76 (MQ=255)
                                         tgtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGCAAGTTCAACATGTCTcac                        >  1:2522502/1‑76 (MQ=255)
                                          gtgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTcaca                       <  1:4392116/76‑1 (MQ=255)
                                           tgtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTcaccc                      >  2:3131360/1‑74 (MQ=255)
                                            gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCCCGGCATAACGCTCCCCCGTAAGGTCCACATTGCTTCCAcc                     >  1:2758127/1‑76 (MQ=255)
                                            gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACAcc                     >  2:5281517/1‑76 (MQ=255)
                                             tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCa                    >  1:1097052/1‑76 (MQ=255)
                                               aaaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAATGCCCACACCCTa                  >  2:2178483/1‑76 (MQ=255)
                                                aaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATaa                 <  2:716158/76‑1 (MQ=255)
                                                 atatAGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAAc                <  1:3080187/76‑1 (MQ=255)
                                                   atCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa              <  1:3467964/76‑1 (MQ=255)
                                                     ccgccAGCGCGCCAACGGGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            <  2:129520/73‑1 (MQ=255)
                                                     cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCGCGCCATAACGAtt            <  2:2224638/76‑1 (MQ=255)
                                                     cGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt            <  1:320179/76‑1 (MQ=255)
                                                      gTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa           <  2:3823882/76‑1 (MQ=255)
                                                       tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCCCACCATAACGAttct          >  2:5173933/1‑74 (MQ=255)
                                                       tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTAt          <  1:3225031/76‑1 (MQ=255)
                                                       tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTAt          <  1:737905/76‑1 (MQ=255)
                                                       tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTAt          <  1:4608158/76‑1 (MQ=255)
                                                        ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCGCCGGTAAGTTCACCATTGCTCACACCATAACGATTATc         <  1:2151257/76‑1 (MQ=255)
                                                        ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATc         >  1:4109049/1‑76 (MQ=255)
                                                          aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCgg       <  1:4870242/76‑1 (MQ=255)
                                                           gCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCCTAACGATTATCggg      <  2:3165596/76‑1 (MQ=255)
                                                            cTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGa     <  2:2668876/76‑1 (MQ=255)
                                                             tCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa    >  1:981195/1‑76 (MQ=255)
                                                             tCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa    <  2:2326044/76‑1 (MQ=255)
                                                              cGCCAACGTGTCGTGCCGCGGCATCTCGTCCACCTGTCAGTTCACCATTGCTCACACCATAACGATTATCgggcat   <  1:4482154/76‑4 (MQ=255)
                                                              cGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt   <  1:5462920/76‑1 (MQ=255)
                                                              cGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt   <  2:4483803/76‑1 (MQ=255)
                                                               gCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATc  <  2:5511040/76‑1 (MQ=255)
                                                               gCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATc  <  2:1240421/76‑1 (MQ=255)
                                                                    ||                                                                     
GCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATC  >  NC_000913/3166867‑3167005

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