Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,471,227 C→T 16.0% E450K (GAA→AAA)  PP_3084 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,471,2270CT16.0% 51.4 / 4.4 25E450K (GAA→AAA) PP_3084ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base T (2/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GTTGCTGTCATAACGCTTGCTGGCCGCCAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCGGCATACAGCATCCAGGCGTCGTTCAGGTCGTACTCGGCGTTGA  >  NC_002947/3471078‑3471376
                                                                                                                                                     |                                                                                                                                                     
gTTGCTGTCATAACGCTTGCTGGCCGCCGGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTtc                                                                                                                                                       <  1:222946/150‑1 (MQ=255)
               cTTGCTGGCCGCCAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATa                                                                                                                                        >  1:436522/1‑150 (MQ=255)
                      gccgccAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAAc                                                                                                                                 >  2:653904/1‑150 (MQ=255)
                      gccgccAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAAc                                                                                                                                 <  2:508572/150‑1 (MQ=255)
                        cgccAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAACATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGCCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAAccc                                                                                                                               <  1:811978/150‑1 (MQ=255)
                             gcgcCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCCGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCcgccgc                                                                                                                          >  1:300012/1‑150 (MQ=255)
                              cgcCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCGAGACCGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCATCGGGACCGCAGAGGAACCCTCcgcg                                                                                                                         >  1:613012/1‑150 (MQ=255)
                               gcCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGt                                                                                                                        >  2:146183/1‑150 (MQ=255)
                                    gTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGt                                                                                                                   <  2:623963/150‑1 (MQ=255)
                                    gTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGt                                                                                                                   >  1:34414/1‑150 (MQ=255)
                                          cGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCGGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCAtt                                                                                                             <  2:300012/150‑1 (MQ=255)
                                              cgcTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCg                                                                                                         <  1:322010/150‑1 (MQ=255)
                                                   ctCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTg                                                                                                                  >  2:9422/1‑136 (MQ=255)
                                                   ctCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTg                                                                                                                  <  1:9422/136‑1 (MQ=255)
                                                        cccACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACAt                                                                                               >  1:82654/1‑150 (MQ=255)
                                                                         ccAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCg                                                                              >  2:781646/1‑150 (MQ=255)
                                                                          cAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGc                                                                             >  1:533051/1‑150 (MQ=255)
                                                                               ttGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTtc                                                                        <  2:436522/150‑1 (MQ=255)
                                                                                     cGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATCCTCGCCGTTCTCACgg                                                                  <  2:613012/150‑1 (MQ=255)
                                                                                                           cAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCgg                                            <  2:533051/150‑1 (MQ=255)
                                                                                                                    cATTGACGCCGGCGACGGCGCTCTTGTTGTCTTTGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATag                                                                                             <  1:756523/92‑1 (MQ=255)
                                                                                                                    cATTGACGCCGGCGACGGCGCTCTTGTTGTCTTTGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATag                                                                                             >  2:756354/1‑92 (MQ=255)
                                                                                                                    cATTGACGCCGGCGACGGCGCTCTTGTTGTCTTTGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATag                                                                                             >  2:756523/1‑92 (MQ=255)
                                                                                                                    cATTGACGCCGGCGACGGCGCTCTTGTTGTCTTTGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATag                                                                                             <  1:756354/92‑1 (MQ=255)
                                                                                                                      ttGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCGGCATACAGCATc                                 <  2:82654/150‑1 (MQ=255)
                                                                                                                                             gttgtCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCGGCATACAGCATCCAGGCGTCGTTCAGGTCGTACTc          >  1:21572/1‑150 (MQ=255)
                                                                                                                                                     cgtcgtGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCGGCATACAGCATCCAGGCGTCGTTCAGGTCGTACTCGGCGTTGa  <  2:21572/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGCTGTCATAACGCTTGCTGGCCGCCAGCGCCTCGTAGGCCGAGCGCTGCTCACCCCACAGGCCGGACAGGCCAAAATTGACCCGGTGGCTGACCGGGCCGGTGGCAAAGCGGCCATTGACGCCGGCGACGGCGCTCTTGTTGTCTTCGTCGTGGGGCGCATAGAGGAACCCGCCGCGTGCAGTGCCATTCCCGTCACTGACATAGAGCGAGGAATACTCGCCGTTCTCACGGGTATGCTTGGCCCCGGCGCCGGCATACAGCATCCAGGCGTCGTTCAGGTCGTACTCGGCGTTGA  >  NC_002947/3471078‑3471376

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