Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,818,324 A→G 8.9% V4304A (GTC→GCC)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,818,3240AG 8.9% 116.1 / 2.5 45V4304A (GTC→GCC) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/2);  total (22/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 = 7.04e-01

CAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTGG  >  NC_002947/4818178‑4818462
                                                                                                                                                   |                                                                                                                                          
cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCt                                                                                                                                         >  2:559238/1‑150 (MQ=255)
    aaCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtc                                                                                                                                     <  2:284790/150‑1 (MQ=255)
     aCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtca                                                                                                                                    <  2:989516/150‑1 (MQ=255)
          cGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCatcatc                                                                                                                               >  1:144738/1‑150 (MQ=255)
           gACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATcc                                                                                                                              <  2:52365/150‑1 (MQ=255)
              tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                                                                           <  1:559238/150‑1 (MQ=255)
              tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                                                                           >  1:445172/1‑150 (MQ=255)
                ccTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCCCTTCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgcgc                                                                                                                         >  1:674073/1‑150 (MQ=255)
                   gattgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACt                                                                                                                      >  1:636515/1‑150 (MQ=255)
                       gatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGccc                                                                                                                  >  1:1100010/1‑150 (MQ=255)
                        atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtg                                                                                                                                     >  2:548311/1‑129 (MQ=255)
                        atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtc                                                                                                                                     <  1:548311/130‑1 (MQ=255)
                        atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACGGCCCt                                                                                                                 >  2:1056283/1‑150 (MQ=255)
                         tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTg                                                                                                                >  1:430728/1‑150 (MQ=255)
                           tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa                                                                                                              <  2:667676/150‑1 (MQ=255)
                               aaTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                                                                          <  2:1100010/150‑1 (MQ=255)
                                aTCCGGATGACCGCCTCAGCCGCGGTGTGCGCAGCCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCGCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGcccc                                                                                                         <  1:1245418/150‑1 (MQ=255)
                                 tCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCt                                                                                                        <  2:445172/150‑1 (MQ=255)
                                   cGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCATTGCCCTGCAGCCCCTGt                                                                                                      <  2:438399/150‑1 (MQ=255)
                                    ggATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGGGAAACTCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTa                                                                                                     >  1:91482/1‑150 (MQ=255)
                                           cGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTc                                                                                              <  2:723496/150‑1 (MQ=255)
                                                           tgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATTCACGGCCAACGCCTCCAcc                                                                              <  2:636515/150‑1 (MQ=255)
                                                                     gACCAACGAAGGGGTCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACt                                                                    <  2:91482/150‑1 (MQ=255)
                                                                            gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                                                                            >  2:9045/1‑87 (MQ=255)
                                                                            gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                                                                            <  1:9045/87‑1 (MQ=255)
                                                                            gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCCCCCGTTACGCACTCAGGTCACTGCGCCCGTCGACTTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACAt                                                             <  2:144738/150‑1 (MQ=255)
                                                                             aaaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATc                                                            >  1:431840/1‑150 (MQ=255)
                                                                              aaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCg                                                           >  1:847415/1‑150 (MQ=255)
                                                                                   gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGCCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc                                                      >  2:204421/1‑150 (MQ=255)
                                                                                      ccAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGc                                                   >  2:55693/1‑150 (MQ=255)
                                                                                                cTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCa                                         <  1:55693/150‑1 (MQ=255)
                                                                                                  gCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGc                                       <  2:491233/150‑1 (MQ=255)
                                                                                                     cATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                                                                       <  2:1188676/99‑1 (MQ=255)
                                                                                                     cATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                                                                       <  2:1188497/99‑1 (MQ=255)
                                                                                                     cATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                                                                       >  1:1188676/1‑99 (MQ=255)
                                                                                                     cATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGGCCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                                                                       >  1:1188497/1‑99 (MQ=255)
                                                                                                        cacaGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgc                                 <  2:430728/150‑1 (MQ=255)
                                                                                                        cacaGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCGCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgc                                 <  1:1056283/150‑1 (MQ=255)
                                                                                                         acaGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAgct                                >  1:1051171/1‑150 (MQ=255)
                                                                                                             ccccTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCtt                            >  1:853835/1‑150 (MQ=255)
                                                                                                              cccTCCCGTTTCCCATCCAGGGCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTc                           <  2:674073/150‑1 (MQ=255)
                                                                                                                                 aGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                                                                 <  2:331098/46‑1 (MQ=255)
                                                                                                                                 aGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                                                                 >  1:331098/1‑46 (MQ=255)
                                                                                                                                     cACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCCCCGTGCTGCACTCGAACATCGCCAGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAGGCGACCTTCg     >  1:177910/1‑150 (MQ=255)
                                                                                                                                        tGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTgg  >  1:1045947/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGGCACGTTCAGCTGTAGCTGCTTCTGCACCCGCGAAGCGATCTGCGTGG  >  NC_002947/4818178‑4818462

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: