Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,712 A→G 9.1% M251V (ATG→GTG)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,316,7120AG 9.1% 143.0 / 4.9 55M251V (ATG→GTG) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGAACTGC  >  NC_002947/4316564‑4316858
                                                                                                                                                    |                                                                                                                                                  
ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                                   >  1:106278/1‑150 (MQ=255)
ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTCTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTTGCTTGGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                                   >  1:429814/1‑150 (MQ=255)
  cccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGc                                                                                                                                                 <  1:801328/150‑1 (MQ=255)
       aaaaCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCtcatc                                                                                                                                            <  2:106278/150‑1 (MQ=255)
         aaCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTcatcat                                                                                                                                          <  1:542161/150‑1 (MQ=255)
                   tCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGc                                                                                                                                <  2:810662/150‑1 (MQ=255)
                      ttCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                                             <  2:728890/150‑1 (MQ=255)
                           ccccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgca                                                                                                                        <  1:635689/150‑1 (MQ=255)
                            cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                                            >  2:106660/1‑127 (MQ=255)
                            cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                                            <  1:106660/129‑3 (MQ=255)
                                   tCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCTCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCAtcttc                                                                                                                >  1:117380/1‑150 (MQ=255)
                                      aggaggCCGTGCTGCCTTTTTTTTCCGCCTGCACCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCg                                                                                                             >  1:662948/1‑150 (MQ=255)
                                       ggaggCCGTGCTGCCTTTTTTTTCGGCCGTTGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGc                                                                                                            <  1:430326/150‑1 (MQ=255)
                                        gaggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTAATCATTGTACCCAGCCACCTGCAGCATCTTCTCGcc                                                                                                           >  2:707711/1‑150 (MQ=255)
                                                  tgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaa                                                                                                 >  2:794030/1‑150 (MQ=255)
                                                  tgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaa                                                                                                 >  2:794088/1‑150 (MQ=255)
                                                      ttttttttCAGCCGGCGCCGTGTTTGTGCCTGACCACAACCGAGTCAACTTTCCCAATTAACTGTCGGCCTGCATCGATGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCAAGCCCAGcaataaa                                                                                             >  2:134231/1‑150 (MQ=255)
                                                       tttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCTCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                            <  2:117380/150‑1 (MQ=255)
                                                         tttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgc                                                                                          <  2:851535/150‑1 (MQ=255)
                                                           tttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                                                        <  1:596399/150‑1 (MQ=255)
                                                             tCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTCCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCCACAGCATCTTCTCGCCAGCCCGGCAATCAAtgcctgc                                                                                      <  1:707711/150‑1 (MQ=255)
                                                                          tgtTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGcc                                                                                                       >  1:153376/1‑120 (MQ=255)
                                                                          tgtTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGcc                                                                                                       <  2:153376/120‑1 (MQ=255)
                                                                               gTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTCACCAGGCAGCGACAGCAACGGCGCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATc                                                                    >  1:101406/1‑150 (MQ=255)
                                                                                        cacaCCCAAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCAGATCGTTtgctgc                                                           >  2:425590/1‑150 (MQ=255)
                                                                                        cacaACCGAGTCAACTTTCCCAATCACCTGACGGCCTGCATCTCTGAGAACGCTGTAACAATGGGAATCATGGTACCCAGCCACCTGCCGCATCTTCTCGCCAGCCCGGCAAACACCCACTGCCTGCAGCGCAGGTATAGCATTtgctgc                                                           >  2:413332/1‑150 (MQ=255)
                                                                                         acaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                                          >  2:13457/1‑150 (MQ=255)
                                                                                           aaCCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGAAACTCAGGCATATCGTTTGCTGCGga                                                        >  1:391501/1‑150 (MQ=255)
                                                                                            aCCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGgag                                                       >  1:738183/1‑150 (MQ=255)
                                                                                               gAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAAc                                                    >  1:417011/1‑150 (MQ=255)
                                                                                                    aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAcc                                                                                                                            >  1:754901/1‑73 (MQ=255)
                                                                                                    aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAcc                                                                                                                            <  2:754901/73‑1 (MQ=255)
                                                                                                     aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGACTCGCGACGCTGGCGAGGGGACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              <  2:662948/150‑1 (MQ=255)
                                                                                                     aCTTTCCCAACTACCTGCCGGCCTGCATCGCTGATCCCGGTGTCACCATGCTCATCATGGGACCCCGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              <  1:304182/150‑1 (MQ=255)
                                                                                                       tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTTCGGAGAACCTTGGGga                                            <  2:738183/150‑1 (MQ=255)
                                                                                                       tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTATCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCgg                                                         <  2:676873/137‑1 (MQ=255)
                                                                                                       tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTATCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCgg                                                         >  1:676873/1‑137 (MQ=255)
                                                                                                                      tCGGCCTGCATTGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt                             <  2:103089/150‑1 (MQ=255)
                                                                                                                        ggCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCGCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTgc                           <  1:13457/150‑1 (MQ=255)
                                                                                                                            tGCATCGCTGAGCACGGTGTCACCATGCTCGTCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCtt                       <  2:429814/150‑1 (MQ=255)
                                                                                                                            tGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                                          >  1:264561/1‑115 (MQ=255)
                                                                                                                            tGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                                          <  2:264561/115‑1 (MQ=255)
                                                                                                                            tGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCtt                       <  2:417011/150‑1 (MQ=255)
                                                                                                                              cATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtca                     >  2:72115/1‑150 (MQ=255)
                                                                                                                                 tgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcatca                  >  1:798215/4‑150 (MQ=255)
                                                                                                                                 cGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATCTCGTTTTCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcatca                  <  1:794030/150‑1 (MQ=255)
                                                                                                                                 cGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTGTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcatca                  <  1:794088/150‑1 (MQ=255)
                                                                                                                                  gCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              >  2:587844/1‑121 (MQ=255)
                                                                                                                                  gCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              <  1:587844/121‑1 (MQ=255)
                                                                                                                                       gcatggTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              <  1:200524/112‑1 (MQ=255)
                                                                                                                                       gCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                              >  2:200524/1‑116 (MQ=255)
                                                                                                                                         ccGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGc          <  1:577062/149‑1 (MQ=255)
                                                                                                                                                 aCCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCATTGCGCTTCATCAGCGCTTGCCGAACTGc  >  1:628132/1‑150 (MQ=255)
                                                                                                                                                  ccATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGc                                                                                                            >  2:624676/1‑43 (MQ=255)
                                                                                                                                                  ccATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGc                                                                                                            <  1:624676/43‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGCCGAACTGC  >  NC_002947/4316564‑4316858

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