Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,663 T→A 8.3% V249D (GTC→GAC)  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,6630TA 8.3% 161.1 / 3.4 57V249D (GTC→GAC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (23/29);  new base A (3/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGC  >  NC_002947/4316518‑4316806
                                                                                                                                                 |                                                                                                                                               
aCGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCCGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCAc                                                                                                                                              >  4:221963/1‑149 (MQ=255)
    cccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTg                                                                                                                                          <  3:377560/149‑1 (MQ=255)
     ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGc                                                                                                                                         <  7:122369/149‑1 (MQ=255)
      cTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                                                        >  3:160379/1‑149 (MQ=255)
      cTGAAAACCCAGGCCTCGTTCGACCCCTCGATACAGGAGGCCGTGCTGCCTTTTTTTTCCGCCTCCGCCGTGTTCGTGCCTGACAACAACCGAGTCAACTTTCCCAATTAACTGTCGGCCTGAATAGCTGAGCACGGTTTCACCATGc                                                                                                                                         >  3:191323/1‑148 (MQ=255)
            aCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTcatcat                                                                                                                                  >  8:168571/1‑149 (MQ=255)
             cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                                 <  5:338513/149‑1 (MQ=255)
             cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                                 <  5:223375/149‑1 (MQ=255)
             cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                                 <  5:249935/149‑1 (MQ=255)
                       gTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                           >  4:360292/1‑129 (MQ=255)
                       gTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                           <  3:360292/129‑1 (MQ=255)
                           gACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCt                                                                                                                   >  4:356068/1‑149 (MQ=255)
                           gACCCCTCGATCCAGGACTCCGTGCTTCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGATCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCt                                                                                                                   <  4:191323/149‑1 (MQ=255)
                             ccccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgc                                                                                                                 <  2:225752/149‑1 (MQ=255)
                              cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGc                                                                                                                        >  7:226813/1‑141 (MQ=255)
                              cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGc                                                                                                                        <  8:226813/141‑1 (MQ=255)
                               ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                               <  3:221963/149‑1 (MQ=255)
                                   gATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCAtc                                                                                                           <  7:168571/149‑1 (MQ=255)
                                        aggaggCCGTGCGTCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTtctc                                                                                                      <  7:51674/149‑1 (MQ=255)
                                           aggCCGTGCTGCC‑TTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCa                                                                                                  <  4:361249/149‑1 (MQ=255)
                                              ccGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGc                                                                                                <  5:249687/149‑1 (MQ=255)
                                                 tgctgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTAc                                                                                                                             <  1:316717/117‑1 (MQ=255)
                                                 tgctgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTAc                                                                                                                             >  2:316717/1‑117 (MQ=255)
                                                        ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                       >  6:197759/1‑148 (MQ=255)
                                                           tttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                                                    <  3:356068/148‑1 (MQ=255)
                                                                  gccgGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcctgc                                                                              >  3:308897/1‑147 (MQ=255)
                                                                  gccgGCGCCGTGTTCGTGCCAGACCACAACCGAGTCAAATTACACAAATACCTCTACGCCTGCATAGCTCAGCACCGAGACAACATGCTAATCATGGAACCCATCCACCTAAACCATATTCACCACACCCCGTCAAACAAtcactgcc                                                                             >  6:123043/1‑148 (MQ=255)
                                                                    cgGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGc                                                                          <  4:160379/149‑1 (MQ=255)
                                                                              ttCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGc                                                                 <  4:308897/148‑1 (MQ=255)
                                                                                 gTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCatat                                                             >  6:117550/1‑149 (MQ=255)
                                                                                       gACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtg                                                       >  3:61387/1‑149 (MQ=255)
                                                                                        aCCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGt                                                          >  6:89952/1‑145 (MQ=255)
                                                                                        aCCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGt                                                          <  5:89952/145‑1 (MQ=255)
                                                                                         ccACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgct                                                     >  8:326954/1‑149 (MQ=255)
                                                                                             aaCCGAGTCAACTTTCACAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCgg                                                 <  5:123043/149‑1 (MQ=255)
                                                                                                      aaCTTTCCCAATTACCTGTCGGCCTTCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTg                                        <  5:117550/149‑1 (MQ=255)
                                                                                                      aaCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                         <  3:102446/148‑1 (MQ=255)
                                                                                                           tCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgag                                   >  1:190649/1‑149 (MQ=255)
                                                                                                               aaTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttga                               <  5:197759/149‑1 (MQ=255)
                                                                                                                      tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAg                                                                           >  2:291067/1‑98 (MQ=255)
                                                                                                                      tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAg                                                                           <  1:291067/98‑1 (MQ=255)
                                                                                                                      tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGAt                         >  1:210575/1‑148 (MQ=255)
                                                                                                                           gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt                     >  8:84829/1‑147 (MQ=255)
                                                                                                                             cTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTgcgc                 <  2:210575/149‑1 (MQ=255)
                                                                                                                              tGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGAAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCt                >  3:203011/1‑149 (MQ=255)
                                                                                                                                cATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgc                                                                                  <  2:346038/81‑1 (MQ=255)
                                                                                                                                cATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgc                                                                                  >  1:346038/1‑81 (MQ=255)
                                                                                                                                 aTCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtc              <  4:61387/148‑1 (MQ=255)
                                                                                                                                  tCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcat            >  5:313779/1‑149 (MQ=255)
                                                                                                                                    gatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaat                                                                                        >  2:275697/3‑71 (MQ=37)
                                                                                                                                    gatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaat                                                                                        <  1:275697/69‑1 (MQ=37)
                                                                                                                                    gatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTtcat            >  3:13558/3‑147 (MQ=255)
                                                                                                                                        aGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgcgc      <  2:190649/149‑1 (MQ=255)
                                                                                                                                        aGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgcg       <  7:84829/148‑1 (MQ=255)
                                                                                                                                            cGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGc  >  7:242978/1‑149 (MQ=255)
                                                                                                                                             ggTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                      >  7:350635/1‑112 (MQ=255)
                                                                                                                                             ggTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                      <  8:350635/112‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGC  >  NC_002947/4316518‑4316806

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