Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,702 C→T 10.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,7020CT10.0% 127.4 / 5.4 51H247H (CAC→CAT) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base T (2/3);  total (23/28)
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.31e-01

GAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGC  >  NC_002947/4316553‑4316844
                                                                                                                                                      |                                                                                                                                              
gAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTTCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCAc                                                                                                                                                >  2:692708/1‑150 (MQ=255)
 aaGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                                               <  1:419897/150‑1 (MQ=255)
 aaGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                                               >  1:355107/1‑150 (MQ=255)
 aaGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                                               <  2:42358/150‑1 (MQ=255)
   gATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGgt                                                                                                                                             <  2:1205321/150‑1 (MQ=255)
       tCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                                         <  2:355107/150‑1 (MQ=255)
       tCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCa                                                                                                                                         >  1:997232/1‑150 (MQ=255)
           ggCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAAACACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCAt                                                                                                                                     <  1:473458/150‑1 (MQ=255)
             cccTGAAACCCCAGGCCTCTTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGc                                                                                                                                   <  2:616512/150‑1 (MQ=255)
               cTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACGTGGCGCCCTGCATCGCTGATCACGGTGTCACCATGCtc                                                                                                                                 >  1:454004/1‑150 (MQ=255)
               cTGAAAACCAAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCATAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                                               >  2:1297367/1‑136 (MQ=255)
               cTGAAAACCAAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCATAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACg                                                                                                                                               <  1:1297367/136‑1 (MQ=255)
                     aCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                           <  1:148617/150‑1 (MQ=255)
                     aCCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATg                                                                                                                           <  1:33260/150‑1 (MQ=255)
                         aGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTAc                                                                                                                       >  1:408774/1‑150 (MQ=255)
                            ccTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCa                                                                                                                    <  1:1324011/150‑1 (MQ=255)
                                 ttCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               <  2:997232/150‑1 (MQ=255)
                                       cccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTgcag                                                                                                         >  2:428027/1‑150 (MQ=255)
                                        ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                              <  1:1370203/128‑3 (MQ=255)
                                        ccTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCATGGTGACACCGTGCtcagc                                                                                                                              >  2:1370203/1‑126 (MQ=255)
                                                  ggaggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTTGTGCCTGACCACAACTG‑AGTCAGCTTTTCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGTTACCCAGCCCCCTGCAGCATCTTCTCGc                                                                                              >  2:746205/1‑150 (MQ=255)
                                                     ggCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTTAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAg                                                                                           >  1:263354/1‑150 (MQ=255)
                                                               ccTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGAACAAACTCGAACTCGA‑AAACCCAATTACCTGTCGGTCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCAACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 >  1:603047/1‑150 (MQ=255)
                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 <  2:191396/148‑1 (MQ=255)
                                                                 ttttttttCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGcaatc                                                                                 >  1:191396/1‑148 (MQ=255)
                                                                         ccgccgGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                                          <  2:48441/147‑1 (MQ=255)
                                                                         ccgccgGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAAtgcct                                                                          >  1:48441/1‑147 (MQ=255)
                                                                                         cGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCTGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCatat                                                       <  1:970863/150‑1 (MQ=255)
                                                                                          gTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATc                                                      >  1:438772/1‑150 (MQ=255)
                                                                                             ccTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGtt                                                   <  2:740742/150‑1 (MQ=255)
                                                                                                gACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTtgc                                                >  2:1060856/1‑150 (MQ=255)
                                                                                                   cacaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               <  1:1378904/84‑1 (MQ=255)
                                                                                                   cacaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCAc                                                                                                               >  2:1378904/1‑84 (MQ=255)
                                                                                                    acaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                            <  2:408774/150‑1 (MQ=255)
                                                                                                    acaACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCg                                            <  1:428027/150‑1 (MQ=255)
                                                                                                         cG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGaa                                       <  1:796929/150‑1 (MQ=255)
                                                                                                              tCAACTTGCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCCGGCATATCGTTTGCTGCGGAGAACCtt                                   <  2:229406/150‑1 (MQ=255)
                                                                                                                aaCTTTCCCAATTACCTGTCGGCCTGCATAGCTGCGCAGAGTGGCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 >  1:983624/1‑150 (MQ=255)
                                                                                                                aaCTGTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 <  2:603047/150‑1 (MQ=255)
                                                                                                                 aCTTTCCCAATTGCCTGCCGGCCTGCACCGCCGAGCACGGTGTCCCCATGCTCATCATGGTACCCAGCCACCTGCAGCCTCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                <  1:622231/150‑1 (MQ=255)
                                                                                                                 aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                <  1:1060856/150‑1 (MQ=255)
                                                                                                                 aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                >  1:532344/1‑150 (MQ=255)
                                                                                                                 aCTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTggg                                >  2:787452/1‑150 (MQ=255)
                                                                                                                   tttCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGga                              >  2:1425098/1‑150 (MQ=255)
                                                                                                                         aaTTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGAAATCCAAGCATATCGTTTGCTGCGGAGAACCTTGGGgagttgag                        >  1:1085439/1‑150 (MQ=255)
                                                                                                                                tGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCa                 <  2:983624/150‑1 (MQ=255)
                                                                                                                                  tCGGCCTGCACCGCTGACAACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGt               <  1:17159/150‑1 (MQ=255)
                                                                                                                                     gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                   >  2:290161/1‑127 (MQ=255)
                                                                                                                                     gCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCtt                                   <  1:290161/127‑1 (MQ=255)
                                                                                                                                           catgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 <  2:585104/118‑1 (MQ=255)
                                                                                                                                           catgatGAGCATGGTGACACCGTGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTgg                                 >  1:585104/6‑123 (MQ=255)
                                                                                                                                               cTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAgc  >  2:22577/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GAAGATGTCACGGCCCTGAAAACCCAGGCCTCGTTCGACCCCTCGATCCAGGAGGCCGTGCTGCCTTTTTTTTCCGCCGGCGCCGTGTTCGTGCCTGACCACAACCG‑AGTCAACTTTCCCAATTACCTGTCGGCCTGCATCGCTGAGCACGGTGTCACCATGCTCATCATGGTACCCAGCCACCTGCAGCATCTTCTCGCCAGCCCGGCAATCAATGCCTGCCAGCACCTCAGGCATATCGTTTGCTGCGGAGAACCTTGGGGAGTTGAGTTGATCAGTGCGCTTCATCAGC  >  NC_002947/4316553‑4316844

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