Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,296,545 G→A 7.7% Q125* (CAG→TAG) ‡ PP_2898 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,296,5450GA 7.7% 145.4 / 2.8 52Q125* (CAG→TAG) ‡PP_2898hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

TGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGT  >  NC_002947/3296404‑3296693
                                                                                                                                                  |                                                                                                                                                    
tgtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACAAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTgc                                                                                                                                              >  2:985577/1‑150 (MQ=255)
 gtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCCAGCCGTTgca                                                                                                                                             <  2:380460/150‑2 (MQ=255)
    cAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTg                                                                                                                                          >  1:937690/1‑150 (MQ=255)
    cAGGGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTg                                                                                                                                          <  2:444741/150‑1 (MQ=255)
      gTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                        >  2:823658/1‑150 (MQ=255)
          gTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGAt                                                                                                                                    <  2:36454/150‑1 (MQ=255)
                   cccTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTc                                                                                                                           >  2:300900/1‑150 (MQ=255)
                    ccTCAATCGCACGGGTAGCTGCTCAGCAGACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCg                                                                                                                          >  2:346614/1‑150 (MQ=255)
                        aaTCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCaa                                                                                                                      >  2:480513/1‑150 (MQ=255)
                         aTCGCACGGGTAGCTGCTCAGCAAACTCGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCaaa                                                                                                                     >  1:759072/1‑150 (MQ=255)
                            gCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGCGCAACGGCCAGCCGTTgca                                                                                                                                             >  2:840371/1‑122 (MQ=255)
                            ccACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGCGCAACGGCCAGCCGTTgca                                                                                                                                             <  1:840371/122‑2 (MQ=255)
                             cACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTgg                                                                                                                 <  1:480513/150‑1 (MQ=255)
                                gggTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTggcgg                                                                                                              <  2:759072/150‑1 (MQ=255)
                                       tgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCGACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCGGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACa                                                                                                       >  2:627716/1‑150 (MQ=255)
                                       tgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACa                                                                                                       <  2:414451/150‑1 (MQ=255)
                                       tgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACa                                                                                                       >  2:84000/1‑150 (MQ=255)
                                            aGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTAcc                                                                                                  <  1:84000/150‑1 (MQ=255)
                                               aaaCTGGCGCTATCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTgg                                                                                               >  1:538462/1‑150 (MQ=255)
                                                aaCTGGCGCTACCCAGATAGGCTTGAGGGGTGGAGGGCATGCGCTGCTGCAACTTGGGCCGATCGACCAGTTCGTGGCATCGGCGGTGTTCACC‑‑ATGAACG‑‑GCG‑TGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTCCAGTACCTggg                                                                                              >  1:253012/1‑150 (MQ=255)
                                                     gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGGAGGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                         <  2:511374/150‑1 (MQ=255)
                                                       gcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGAtt                                                                                       <  1:533811/150‑1 (MQ=255)
                                                       gcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGAtt                                                                                       <  1:678244/150‑1 (MQ=255)
                                                          aCCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAg                                                                                    <  2:208872/150‑1 (MQ=255)
                                                           cccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                                   <  2:736110/150‑1 (MQ=255)
                                                           cccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGt                                                                                   <  1:985577/150‑1 (MQ=255)
                                                            ccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTg                                                                                  >  1:822837/1‑150 (MQ=255)
                                                            ccAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTg                                                                                  >  2:563271/1‑150 (MQ=255)
                                                              aGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGAAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGc                                                                                >  2:595074/1‑150 (MQ=255)
                                                                    gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTg                                                                                                                              <  1:712149/98‑1 (MQ=255)
                                                                    gCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTg                                                                                                                              >  2:712149/1‑98 (MQ=255)
                                                                                gcgcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTGGGCAAATTTTCGCCGTAACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                              <  1:300900/150‑1 (MQ=255)
                                                                                gcgcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                              <  1:419009/150‑1 (MQ=255)
                                                                                          tCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGttt                                                    >  2:862246/1‑150 (MQ=255)
                                                                                               ggcAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCgg                                               <  1:823658/150‑1 (MQ=255)
                                                                                               ggcAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCgg                                               >  1:422324/1‑150 (MQ=255)
                                                                                                             aTCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTgg                                                                                 <  1:329383/102‑1 (MQ=255)
                                                                                                             aTCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTgg                                                                                 >  2:329383/1‑102 (MQ=255)
                                                                                                                 aCCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCgg                             >  2:1001321/1‑150 (MQ=255)
                                                                                                                  ccAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCggg                            <  1:602840/150‑1 (MQ=255)
                                                                                                                   cag‑‑c‑agGGA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                           >  1:375706/1‑150 (MQ=255)
                                                                                                                    ag‑‑c‑agGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTTATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCt                          <  2:538462/150‑1 (MQ=255)
                                                                                                                     g‑‑c‑agGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGGtg                         >  1:680161/1‑150 (MQ=255)
                                                                                                                          agGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCGGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCtgtg                       <  2:287057/150‑1 (MQ=255)
                                                                                                                            gCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCNGCCTGCGCCAGCGTTAGTCGAGCCAGGTCGGGCAAGCGGGCTGTggg                     <  1:145949/150‑1 (MQ=255)
                                                                                                                              a‑caGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                          <  2:1057674/127‑1 (MQ=255)
                                                                                                                              a‑caGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAg                                          >  1:1057674/1‑127 (MQ=255)
                                                                                                                                     cagcgCAACGGCCAGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGATACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTg              >  2:730183/5‑150 (MQ=255)
                                                                                                                                       gTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGaa            >  1:1032631/1‑150 (MQ=255)
                                                                                                                                        tGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTCTCTTCCGCTCTGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGaaa           <  1:346614/150‑1 (MQ=255)
                                                                                                                                            aCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTa       >  1:954383/1‑150 (MQ=255)
                                                                                                                                                 tGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGt  <  2:937690/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAG‑‑C‑AGGCA‑CAGC‑GCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAATGTACCGGT  >  NC_002947/3296404‑3296693

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