Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,658 C→T 7.0% H247H (CAC→CAT)  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,6580CT 7.0% 162.4 / 2.1 57H247H (CAC→CAT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (23/30);  new base T (2/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

TGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGC  >  NC_002947/4316514‑4316806
                                                                                                                                                |                                                                                                                                                    
tGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgt                                                                                                                                                   <  6:350015/148‑1 (MQ=255)
 gTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgtg                                                                                                                                                  >  6:338513/1‑148 (MQ=255)
    aCGGCCGTGAAATCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGTCGCCGTGTTCGTGCCTGACGACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgt                                                                                                                                                   <  4:204531/144‑1 (MQ=255)
    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:223375/149‑1 (MQ=255)
                 cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                                 <  5:338513/149‑1 (MQ=255)
                 cccAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                                 <  5:249935/149‑1 (MQ=255)
                           gTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                           <  3:360292/129‑1 (MQ=255)
                           gTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                           >  4:360292/1‑129 (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                                                                           <  1:291067/98‑1 (MQ=255)
                                                                                                                          tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAg                                                                           >  2:291067/1‑98 (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)
                                                                                                                                                |                                                                                                                                                    
TGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGCGCTTGC  >  NC_002947/4316514‑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: ‑