Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,827,099 C→A 11.6% S1379I (AGT→ATT) ‡ pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,827,0990CA11.6% 120.8 / 6.6 43S1379I (AGT→ATT) ‡pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base A (3/2);  total (22/21)
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

TGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAG  >  NC_002947/4826962‑4827247
                                                                                                                                         |                                                                                                                                                    
tGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCt                                                                                                                                          >  2:259420/1‑150 (MQ=255)
    gTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCATCACGCGGTCACGCCGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCGCTTGCCTGCGGCTCCTg                                                                                                                                      <  2:94270/150‑1 (MQ=255)
    gTTGATGGCGGTTTCAGTCGGGCCATAACGGATGTGCAGTAGCACCTGAGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTg                                                                                                                                      >  2:434469/1‑150 (MQ=255)
      tGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGa                                                                                                                                    <  2:253880/150‑1 (MQ=255)
       gATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTccctgat                                                                                                                                   <  2:429389/150‑6 (MQ=255)
           gCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAac                                                                                                                               <  1:290795/150‑1 (MQ=255)
             ggTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCt                                                                                                                                       >  2:272354/1‑140 (MQ=255)
             ggTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCt                                                                                                                                       <  1:272354/140‑1 (MQ=255)
             ggTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAacac                                                                                                                             <  1:30689/150‑1 (MQ=255)
                    gTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAAc                                                                                                                      >  1:549119/1‑150 (MQ=255)
                    gTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAAc                                                                                                                      >  1:303859/1‑150 (MQ=255)
                       gggCCACACCGGTTGCGCAGTTGCACCTGCGGTTGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAgc                                                                                                                   <  2:497629/150‑1 (MQ=255)
                            aTACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACa                                                                                                                     >  1:142472/1‑143 (MQ=255)
                            aTACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACa                                                                                                                     <  2:142472/143‑1 (MQ=255)
                                   ttGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCTAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAacac                                                                                                                             <  1:110751/128‑1 (MQ=255)
                                   ttGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAacac                                                                                                                             >  2:110751/1‑128 (MQ=255)
                                          gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCa                                                                                                   <  2:592086/147‑1 (MQ=255)
                                          gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCa                                                                                                   >  1:592086/1‑147 (MQ=255)
                                          gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcagta                                                                                                                                          >  2:572696/1‑104 (MQ=255)
                                          gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCAAGTGCCTGcagta                                                                                                                                          <  1:572696/108‑5 (MQ=255)
                                                   gCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCAccc                                                                                       <  2:591485/150‑1 (MQ=255)
                                                          caCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGc                                                                                >  1:237847/1‑150 (MQ=255)
                                                                        ggTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGc                                                                    <  2:256917/148‑1 (MQ=255)
                                                                        ggTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGc                                                                    >  1:256917/1‑148 (MQ=255)
                                                                                            aaCGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGa                                              <  2:601239/150‑1 (MQ=255)
                                                                                                 ttCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACACCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTc                                         >  1:453590/1‑150 (MQ=255)
                                                                                                    gCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAcc                                      >  1:9923/1‑150 (MQ=255)
                                                                                                           gagaACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCGACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGcc                                                                               >  1:407458/1‑102 (MQ=255)
                                                                                                           gagaACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCGACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGcc                                                                               <  2:407458/102‑1 (MQ=255)
                                                                                                             gaACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCTTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGc                             <  1:12449/150‑1 (MQ=255)
                                                                                                                   ccAGCGCAAGCTACTGCAGTCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAACGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAACCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCAc                       >  2:1645/1‑150 (MQ=255)
                                                                                                                    cAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTcgc                                                >  1:170359/1‑124 (MQ=255)
                                                                                                                    cAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTcgc                                                <  2:170359/124‑1 (MQ=255)
                                                                                                                     aGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCgg                                                     >  2:111163/1‑118 (MQ=255)
                                                                                                                     aGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCgg                                                     <  1:111163/118‑1 (MQ=255)
                                                                                                                             gagccgCAGGCAAGTGCCTGCAGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAg             >  1:475046/6‑150 (MQ=255)
                                                                                                                                ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAg                                                                                                  >  2:450399/3‑62 (MQ=38)
                                                                                                                                ccgCAGGCAAGTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAg                                                                                                  <  1:450399/60‑1 (MQ=38)
                                                                                                                                cTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAc                                       >  2:398902/1‑121 (MQ=255)
                                                                                                                                cTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAc                                       <  1:398902/121‑1 (MQ=255)
                                                                                                                                 tGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGGACCAGCTTGCAGCCgg         >  2:568897/1‑150 (MQ=255)
                                                                                                                                  gCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGt        <  2:549119/150‑1 (MQ=255)
                                                                                                                                        aCTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAg  >  2:235305/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TGACGTTGATGGCGGTTTCAGTCGGGCCATACCGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAG  >  NC_002947/4826962‑4827247

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: