Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,378,669 A→G 7.1% M21T (ATG→ACG)  PP_2979 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,378,6690AG 7.1% 159.2 / 3.0 57M21T (ATG→ACG) PP_2979hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

AAGATCACATCGTCCAGGTAGCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCATCTTGGTGTAGC  >  NC_002947/3378528‑3378816
                                                                                                                                             |                                                                                                                                                   
aaGATCACATCGTCCAGGTAGCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAAc                                                                                                                                              <  7:209107/149‑1 (MQ=255)
  gATCACATCGTCCAGGTAGCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACtt                                                                                                                                            >  6:8845/1‑149 (MQ=255)
   aTCACATCCGACATGTATCCAAGGTCCGGTATGTACGCTGGAATCAGATCGTTTTGGCGTACGGCATAAGCTATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTg                                                                                                                                           <  2:295451/149‑1 (MQ=255)
            tCCAGGTAGCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAg                                                                                                                                  >  2:55962/1‑149 (MQ=255)
                    gCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTgg                                                                                                                          >  2:110811/1‑149 (MQ=255)
                     ccAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGc                                                                                                                         >  7:108744/1‑149 (MQ=255)
                               gAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCcgcg                                                                                                                >  3:53512/1‑148 (MQ=255)
                                    aaGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCg                                                                                                          >  2:74046/1‑149 (MQ=255)
                                      gTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGtt                                                                                                        >  6:23653/1‑149 (MQ=255)
                                      gTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGtt                                                                                                        >  8:55328/1‑149 (MQ=255)
                                      gTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTATGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGCGACTTGGCCCAGCCGCGAACTCGtt                                                                                                        >  6:19837/1‑149 (MQ=255)
                                       tCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTg                                                                                                       >  8:150023/1‑149 (MQ=255)
                                       tCTGGAATCAGATCGATTGGGATTAGGGCATTAGCTATTACGATGATCGCTATAACCTTCGGCATCCATGGCGTTTGCGGATTCCGGCAACTAAACCACAGCATCTTAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCTAACTCGTTg                                                                                                       <  8:226768/149‑1 (MQ=255)
                                         tGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGc                                                                                                      >  4:316506/1‑148 (MQ=255)
                                         tGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGc                                                                                                      >  1:188171/1‑148 (MQ=255)
                                               cAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCGTCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGt                                                                                               >  2:172343/1‑149 (MQ=255)
                                               cAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCGTCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGt                                                                                               >  2:167197/1‑149 (MQ=255)
                                                aGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGt                                                                                               <  3:316506/148‑1 (MQ=255)
                                                aGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGt                                                                                               <  5:13782/148‑1 (MQ=255)
                                                 gATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTa                                                                                             <  4:53512/149‑1 (MQ=255)
                                                 gATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTa                                                                                             <  5:187336/149‑1 (MQ=255)
                                                  aTCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTAc                                                                                            >  5:209692/1‑149 (MQ=255)
                                                  aTCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTAc                                                                                            >  6:228713/1‑149 (MQ=255)
                                                     gATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCTGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTc                                                                                          <  5:8845/148‑1 (MQ=255)
                                                      aTTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCAc                                                                                        >  8:65093/1‑149 (MQ=255)
                                                       ttGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCCAACTCGTTGCAATACGTTACTCACt                                                                                       <  6:88488/149‑1 (MQ=255)
                                                             ttAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACt                                                                                           >  7:96634/1‑139 (MQ=255)
                                                             ttAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACt                                                                                           <  8:96634/139‑1 (MQ=255)
                                                                      aaGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGa                                                                                >  2:222621/1‑141 (MQ=255)
                                                                      aaGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGa                                                                                <  1:222621/141‑1 (MQ=255)
                                                                               cGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAgc                                                               <  1:110811/149‑1 (MQ=255)
                                                                                gATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAgcg                                                              >  5:89217/1‑149 (MQ=255)
                                                                                gATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTAACTTTCTAGAATCCGTAGCCAACCAgcg                                                              >  5:147008/1‑149 (MQ=255)
                                                                                gATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTAACTTTCTAGAATACGTAGACAACCAgcg                                                              >  5:142229/1‑149 (MQ=255)
                                                                                           taACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCtt                                                    <  7:65093/148‑1 (MQ=255)
                                                                                             aCCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCg                                                 <  7:55328/149‑1 (MQ=255)
                                                                                              ccTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTTCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCgg                                                <  5:228713/149‑1 (MQ=255)
                                                                                               cTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATAATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGt                                               >  4:231241/1‑149 (MQ=255)
                                                                                                   ggCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCGTCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGc                                            <  1:167197/148‑1 (MQ=255)
                                                                                                   ggCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCGTCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGc                                            <  1:172343/148‑1 (MQ=255)
                                                                                                   ggCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGcc                                           <  5:23653/149‑1 (MQ=255)
                                                                                                    gCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTg                                                   >  6:327776/1‑140 (MQ=255)
                                                                                                    gCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTg                                                   <  5:327776/140‑1 (MQ=255)
                                                                                                    gCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGcc                                           >  5:297750/1‑148 (MQ=255)
                                                                                                       aacATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGt                                       <  8:108744/147‑1 (MQ=255)
                                                                                                       aCCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGt                                       <  6:209692/149‑1 (MQ=255)
                                                                                                         cATGGCGTTTGCGGATGCCGGCAACTAAACCAGAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGttt                                     <  7:138939/149‑1 (MQ=255)
                                                                                                         cATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAgcgc                                                             >  5:143002/1‑125 (MQ=255)
                                                                                                         cATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAgcgc                                                             <  6:143002/125‑1 (MQ=255)
                                                                                                                          ccGGCAACTAAACCACAGCCTCATAACTTGCTGTTGCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCg                     <  8:260027/148‑1 (MQ=255)
                                                                                                                          ccGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCgg                    <  7:150023/149‑1 (MQ=255)
                                                                                                                             gCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAgg                 >  8:180814/1‑149 (MQ=255)
                                                                                                                               aaCTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGcc               <  6:297750/149‑1 (MQ=255)
                                                                                                                               aaCTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCACCGAACTCGTTGCAATACGATACACACTTTCTAGAATCCGTATCCAACCAGCGCCAATATCTTGCTGTGGCCGAGTTTGTACGAGACTTAGCCGGAGGc                >  1:323182/1‑148 (MQ=255)
                                                                                                                                aCTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCa              <  1:55962/149‑1 (MQ=255)
                                                                                                                                    aaCCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCGACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCATCtt          >  4:229476/1‑149 (MQ=255)
                                                                                                                                    aaCCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCATCtt          >  8:20117/1‑149 (MQ=255)
                                                                                                                                            catcatAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCATCTTGGTGTAGc  <  2:188171/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AAGATCACATCGTCCAGGTAGCCAAGGACCGGAATGAAGTCTGGAATCAGATCGATTGGGCTTAGGGCATAAGCGATTACGATCATCGCTATAACCTTCGGCAACCATGGCGTTTGCGGATGCCGGCAACTAAACCACAGCATCATAACTTGCTGTTTCAGGGCCTTGGCCCAGCCGCGAACTCGTTGCAATACGTTACTCACTTTCTAGAATCCGTAGCCAACCAGCGCCAATAGCTTGCGGTGGCCGAGTTTGTCGGCGACGTAGCCGGAGGCCATCTTGGTGTAGC  >  NC_002947/3378528‑3378816

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