Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,686 C→T 7.5% A242V (GCC→GTC)  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,6860CT 7.5% 146.6 / 2.6 53A242V (GCC→GTC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base T (2/2);  total (23/30)
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.26e-01

TTACGCGCTCGGCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGT  >  NC_002947/4316537‑4316823
                                                                                                                                                     |                                                                                                                                         
ttACGCGCTCGGCAGCGAAGATGTCACGGCCCTGAAAACCCCGGCCTCATTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGc                                                                                                                                           >  2:156816/1‑150 (MQ=255)
ttACGCGCTCGGCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGc                                                                                                                                           <  1:675896/150‑1 (MQ=255)
          ggCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCt                                                                                                                                 >  1:1054764/1‑150 (MQ=255)
          ggCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCt                                                                                                                                 <  2:724923/150‑1 (MQ=255)
          ggCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCt                                                                                                                                 <  2:726587/150‑1 (MQ=255)
             aGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCt                                                                                                                                         <  1:955600/139‑1 (MQ=255)
             aGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCt                                                                                                                                         >  2:955600/1‑139 (MQ=255)
                gAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCAc                                                                                                                           >  2:412996/1‑150 (MQ=255)
                caaGACGTCACAGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCTTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCAc                                                                                                                           <  2:95696/149‑1 (MQ=255)
                 aaGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTACTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                          <  2:660088/150‑1 (MQ=255)
                   gATGTCACGGCCCTGAAAACCCAGTCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTTAGCACGgt                                                                                                                        <  2:1077357/150‑1 (MQ=255)
                   gATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgt                                                                                                                        <  2:723615/150‑1 (MQ=255)
                   gATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgt                                                                                                                        <  1:649170/150‑1 (MQ=255)
                      gTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTc                                                                                                                     <  2:1054764/150‑1 (MQ=255)
                       tCACGGCCCTGAAAACCCAGGCCTCGTTCGGCCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                    <  2:32900/150‑1 (MQ=255)
                        cACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCAc                                                                                                                   <  2:54391/150‑1 (MQ=255)
                         aCGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCAcc                                                                                                                  >  1:718674/1‑150 (MQ=255)
                         aCGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCAcc                                                                                                                  >  1:162657/1‑150 (MQ=255)
                          cGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCACCa                                                                                                                 <  2:718674/150‑1 (MQ=255)
                          cGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCACCa                                                                                                                 <  2:162657/150‑1 (MQ=255)
                           ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGGCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGGCCGGTGTCACCAt                                                                                                                >  2:603620/1‑150 (MQ=255)
                              ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                             >  1:970815/1‑150 (MQ=255)
                              ccTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCt                                                                                                             <  2:161385/150‑1 (MQ=255)
                                tGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGcc                                                                                                                                          >  2:948557/1‑119 (MQ=255)
                                tGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGcc                                                                                                                                          <  1:948557/119‑1 (MQ=255)
                                tGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCtca                                                                                                           >  1:897500/1‑150 (MQ=255)
                                  aaaaCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                         <  2:1040988/150‑3 (MQ=255)
                                           gCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTAccc                                                                                                >  1:380552/1‑150 (MQ=255)
                                                         cTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTACCAATTACCTGGCGGCTTGCTTCTCTGATTACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagca                                                                                  >  2:678692/1‑150 (MQ=255)
                                                          tCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCAt                                                                                 <  2:895004/150‑1 (MQ=255)
                                                                   gaggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGcc                                                                        >  2:190929/1‑150 (MQ=255)
                                                                     ggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAg                                                                      >  2:1100356/1‑150 (MQ=255)
                                                                          tgctgcCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCgg                                                                 <  1:163228/150‑1 (MQ=255)
                                                                               ccTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                            <  1:412996/150‑1 (MQ=255)
                                                                               ccTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                            >  1:208692/1‑150 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCa                                                                                           <  1:35221/117‑1 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCa                                                                                           >  2:35221/1‑117 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaa                                                          <  2:208692/150‑1 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaa                                                          >  1:764685/1‑150 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTACATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagc                                                                                   >  2:850110/1‑125 (MQ=255)
                                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTACATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcagc                                                                                   <  1:850110/125‑1 (MQ=255)
                                                                                      tttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                     <  2:1069737/150‑1 (MQ=255)
                                                                                        tCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcctgc                                                   <  1:156816/150‑1 (MQ=255)
                                                                                          cgccgGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCa                                                 <  1:1100356/150‑1 (MQ=255)
                                                                                                 gcCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                         <  2:497748/135‑1 (MQ=255)
                                                                                                 gcCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatcaat                                                         >  1:497748/1‑135 (MQ=255)
                                                                                                                 aCCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgct                          >  2:234665/1‑150 (MQ=255)
                                                                                                                   cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                        >  1:707523/1‑150 (MQ=255)
                                                                                                                   cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                        <  2:897500/150‑1 (MQ=255)
                                                                                                                   cacaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgctgc                        >  1:551549/1‑150 (MQ=255)
                                                                                                                    gcaACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                       >  1:588301/2‑150 (MQ=255)
                                                                                                                             tCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt              >  2:251863/1‑150 (MQ=255)
                                                                                                                              aatcTTTCCCAAGCACCTGTCTGCCTACACCGCCGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCATAATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTg             <  1:76937/147‑1 (MQ=255)
                                                                                                                                  tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGga         <  1:678692/150‑1 (MQ=255)
                                                                                                                                     cccAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCCCCCTCAACATCCTACTTCGCCACCAGCAGCATCTTCTCGCCAGCCCGTCAATTGATGCGTGCCATCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagtt      <  2:970815/150‑1 (MQ=255)
                                                                                                                                       cAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCATGCATATCGTTTGCTGCGGAGAACCTTGGGgagttga    >  1:396778/1‑150 (MQ=255)
                                                                                                                                         aTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgagt  <  2:764685/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
TTACGCGCTCGGCAGCGAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCGAGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGT  >  NC_002947/4316537‑4316823

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