Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,823,340 C→G 8.2% G2617A (GGC→GCC) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,823,3400CG 8.2% 147.7 / 4.1 49G2617A (GGC→GCC) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base G (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

TCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAC  >  NC_002947/4823196‑4823487
                                                                                                                                                |                                                                                                                                                   
tCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGccc                                                                                                                                                   <  6:332666/147‑1 (MQ=255)
tCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCgg                                                                                                                                                 >  2:165558/1‑149 (MQ=255)
  gAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGc                                                                                                                                                >  7:356122/1‑148 (MQ=255)
  gAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGc                                                                                                                                                <  8:356122/148‑1 (MQ=255)
    gcagcaCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGGGCGGCAAg                                                                                                                                             <  6:7645/149‑1 (MQ=255)
      agcaCGAGATGCGCAGGCAGCATGTAGTGGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAACTGCTTGCCGCCCGGCAAGTc                                                                                                                                           <  2:312574/149‑1 (MQ=255)
         aCCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAg                                                                                                                                        <  2:107568/149‑1 (MQ=255)
             gATGCGCAGGCACCATGTAGTCGGGCACGTGCAGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCACTTGTTGCACCTCGTTACCGGCGGTGGCATGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGAAAGTCCAgcgc                                                                                                                                     >  1:44942/1‑148 (MQ=255)
               tGCGCAGGCACCGTGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAg                                                                                                                                  <  5:59979/149‑1 (MQ=255)
               tGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAg                                                                                                                                  <  3:177093/149‑1 (MQ=255)
                            tGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAg                                                                                                                                             <  2:207877/125‑1 (MQ=255)
                                 tCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGAc                                                                                                                <  7:36502/149‑1 (MQ=255)
                                    ggCAGGTGCTGGCGGGCATCGGCTTTCCCCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATc                                                                                                             >  7:197827/1‑149 (MQ=255)
                                                  ggCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCa                                                                                               <  1:228314/149‑1 (MQ=255)
                                                     aTCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggc                                                                                            >  6:48192/1‑149 (MQ=255)
                                                        ggCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAggcggc                                                                                         >  3:97176/1‑149 (MQ=255)
                                                         gCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGct                                                                                        >  8:272399/1‑149 (MQ=255)
                                                           tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctct                                                                                      >  2:76494/1‑149 (MQ=255)
                                                           tttCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCATGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGctct                                                                                      >  1:258733/1‑149 (MQ=255)
                                                                        cgcAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAg                                                                         >  3:365675/1‑149 (MQ=255)
                                                                         gcagcaGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGt                                                                        >  7:242373/1‑149 (MQ=255)
                                                                             cagTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGa                                                                    >  2:211590/1‑149 (MQ=255)
                                                                              agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGGCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCTGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                   <  2:44942/149‑1 (MQ=255)
                                                                              agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGCCATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                   >  7:82528/1‑149 (MQ=255)
                                                                              agTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGAt                                                                   <  8:197827/149‑1 (MQ=255)
                                                                                ttgttgCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGc                                                                 <  1:76494/149‑1 (MQ=255)
                                                                                         gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                        <  8:242373/149‑1 (MQ=255)
                                                                                         gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                        <  3:204810/149‑1 (MQ=255)
                                                                                         gtcgtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCAGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCcgc                                                        <  2:209087/149‑1 (MQ=255)
                                                                                         gtagtTACGGGCGGTGGCCTGCGTGCAAAGCTGATAGCCCACCAGCTGGTTGCGGCCCGGCAAGTCGAGCGCCAGCACGAGGTCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCGCAGTTCGATGGGGAAGCcg                                                         <  4:365675/145‑1 (MQ=255)
                                                                                            gtTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGa                                                     >  8:268616/1‑149 (MQ=255)
                                                                                               aCCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATc                                                   <  7:268616/148‑1 (MQ=255)
                                                                                                 cggcggTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTc                                                <  1:165558/149‑1 (MQ=255)
                                                                                                   gcggTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCAc                                              <  5:48192/149‑1 (MQ=255)
                                                                                                           ctgcctgcACACCAGATAGCCCACCATCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACCGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCTATATCTCCCAGTTCGATGAGGAAGCCGCGGATCTTCACCTGCTGGt                                      >  3:62778/1‑149 (MQ=255)
                                                                                                            tgcctgcACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTc                                     >  5:334420/1‑149 (MQ=255)
                                                                                                                   acacCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACgg                              >  6:267548/1‑149 (MQ=255)
                                                                                                                       cAGATAGACCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCCTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATGTCGCCCAGTACGATGCGGAAGCCGCGGATCTACACCTGCTGGACGGCCCGGCCg                           >  5:260431/1‑148 (MQ=255)
                                                                                                                         gATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACa                        <  4:97176/149‑1 (MQ=255)
                                                                                                                            aGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATa                      >  2:322429/1‑148 (MQ=255)
                                                                                                                                caccaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTcca                 >  6:46348/1‑149 (MQ=255)
                                                                                                                                 accaGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTccac                >  1:148450/1‑149 (MQ=255)
                                                                                                                                     gctggatGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACTGCCGACATACTCCACCAgt            <  8:82528/143‑2 (MQ=255)
                                                                                                                                       tgctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGc           >  3:171228/1‑148 (MQ=255)
                                                                                                                                       gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAg             <  6:88484/144‑1 (MQ=255)
                                                                                                                                       gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAg             >  5:88484/3‑146 (MQ=255)
                                                                                                                                       gactTGCCGGGCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGc           >  7:57127/3‑148 (MQ=255)
                                                                                                                                        gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATAc                     <  7:241760/137‑1 (MQ=255)
                                                                                                                                        gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATAc                     >  8:241760/1‑137 (MQ=255)
                                                                                                                                        gctTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCa         >  3:27871/1‑149 (MQ=255)
                                                                                                                                               gccCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAc  >  1:241984/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCGAGCAGCACCAGATGCGCAGGCACCATGTAGTCGGGCAGGTGCTGGCGGGCATCGGCTTTCACCGCTTCGCGCAGCAGTTGTTGCACGTCGTTACCGGCGGTGGCCTGCCTGCACACCAGATAGCCCACCAGCTGCTTGCCGCCCGGCAAGTCCAGCGCCAGCACGACGGCCTCGTCGACATCCTCGTGCGCTAGCAGGCGGCTCTCGATTTCGCCCAGTTCGATGCGGAAGCCGCGGATCTTCACCTGCTGGTCGGCACGGCCGACATACTCCACCAGGCCATCGGTAC  >  NC_002947/4823196‑4823487

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: