Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,130,169 G→C 8.8% A107G (GCC→GGC)  PP_5378 ← metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,130,1690GC 8.8% 171.6 / 6.1 56A107G (GCC→GGC) PP_5378metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base C (3/2);  total (31/25)
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.73e-01

GCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGA  >  NC_002947/6130026‑6130313
                                                                                                                                               |                                                                                                                                                
gCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGAt                                                                                                                                            >  1:190753/1‑150 (MQ=255)
           tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGcc                                                                                                                                 >  2:700352/1‑150 (MQ=255)
           tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGcc                                                                                                                                 >  2:514945/1‑150 (MQ=255)
             aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTggg                                                                                                                                        >  2:87818/1‑141 (MQ=255)
             aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTggg                                                                                                                                        <  1:87818/141‑1 (MQ=255)
             aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCgg                                                                                                                               >  2:353738/1‑150 (MQ=255)
              aGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGa                                                                                                                              >  1:297781/1‑150 (MQ=255)
                gAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATa                                                                                                                            >  2:217928/1‑150 (MQ=255)
                   aGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGATAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAg                                                                                                                         <  1:704480/150‑1 (MQ=255)
                      gACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAgcgc                                                                                                                      <  2:719296/150‑1 (MQ=255)
                         aaTACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTg                                                                                                                   <  2:297781/150‑1 (MQ=255)
                         aaTACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTg                                                                                                                   >  1:771066/1‑150 (MQ=255)
                              cccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                              >  2:556380/1‑150 (MQ=255)
                              aacAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                              >  1:353151/3‑150 (MQ=255)
                               ccAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAg                                                                                                             >  2:13127/1‑150 (MQ=255)
                                         tGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTa                                                                                                   <  1:353738/150‑1 (MQ=255)
                                              aCCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAc                                                                                              >  2:204904/1‑150 (MQ=255)
                                                ccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTggg                                                                                                                                        >  2:628583/1‑106 (MQ=255)
                                                ccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTggg                                                                                                                                        <  1:628583/106‑1 (MQ=255)
                                                ccATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCAccc                                                                                            >  2:197278/1‑150 (MQ=255)
                                                     ctctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTggg                                                                                       <  1:13127/150‑1 (MQ=255)
                                                       ctTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGc                                                                                     <  1:217928/150‑1 (MQ=255)
                                                            cccAAGCGGGCATACCGGTAGCCTTGATACCGTGCTTCATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGt                                                                                >  2:178477/1‑150 (MQ=255)
                                                             ccAAGCGGACATACCGGTAGCCTTGATACCGTGGTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTg                                                                               >  2:428522/1‑150 (MQ=255)
                                                                           cGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                      >  1:367252/1‑113 (MQ=255)
                                                                           cGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                      <  2:367252/113‑1 (MQ=255)
                                                                             gTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATa                                                               >  2:336823/1‑150 (MQ=255)
                                                                                  cTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCAtt                                                          <  2:353151/150‑1 (MQ=255)
                                                                                  cTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCAtt                                                          >  2:535407/1‑150 (MQ=255)
                                                                                  cTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCAtc                                                          <  2:547065/150‑2 (MQ=255)
                                                                                     gATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACACAGCATTGTa                                                       <  1:197278/150‑1 (MQ=255)
                                                                                      aTACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAg                                                      <  1:556380/150‑1 (MQ=255)
                                                                                          cGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                      >  1:435350/1‑98 (MQ=255)
                                                                                          cGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGtt                                                                                                      <  2:435350/98‑1 (MQ=255)
                                                                                          cGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTac                                                  >  2:276523/1‑150 (MQ=255)
                                                                                                        cAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGGTTGGGCGGTGACGGACAAAGCGCTTTGCTGAGCTCGGATATATCCAGCCCCGGGGCGAGGTAGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACgg                                    >  2:10202/1‑150 (MQ=255)
                                                                                                              gCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAgg                              <  2:190753/150‑1 (MQ=255)
                                                                                                              gCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAgg                              <  1:700352/150‑1 (MQ=255)
                                                                                                                    gCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAAt                                                                      <  1:554459/104‑1 (MQ=255)
                                                                                                                    gCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAAt                                                                      >  2:554459/1‑104 (MQ=255)
                                                                                                                          tCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTgg                                                                                        >  1:697730/1‑80 (MQ=255)
                                                                                                                          tCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTgg                                                                                        <  2:697730/80‑1 (MQ=255)
                                                                                                                          tCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCttt                  <  1:276523/150‑1 (MQ=255)
                                                                                                                            aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                       >  1:781223/1‑111 (MQ=255)
                                                                                                                            aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                       <  2:781223/111‑1 (MQ=255)
                                                                                                                                cTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                              >  1:219612/1‑52 (MQ=255)
                                                                                                                                cTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGa                                                                                                              <  2:219612/52‑1 (MQ=255)
                                                                                                                                    gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGtt             <  1:767763/143‑1 (MQ=255)
                                                                                                                                    gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGtt             >  2:767763/3‑145 (MQ=255)
                                                                                                                                    gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGaa        >  1:771650/3‑150 (MQ=255)
                                                                                                                                    gcccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGaa        >  2:497242/3‑150 (MQ=255)
                                                                                                                                     aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGt                            <  2:642711/129‑1 (MQ=255)
                                                                                                                                     aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGt                            >  1:642711/1‑129 (MQ=255)
                                                                                                                                       cAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCt     >  2:675343/1‑150 (MQ=255)
                                                                                                                                        aaTCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTc    <  1:428522/150‑1 (MQ=255)
                                                                                                                                          tCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGa  <  1:497242/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTCGA  >  NC_002947/6130026‑6130313

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: ‑