Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 778,365 A→G 100% N472S (AAC→AGC)  PP_0669 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947778,3650AG94.6% 109.4 / ‑1.7 37N472S (AAC→AGC) PP_0669ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (12/23);  total (13/24)
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.21e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCGGGTGGCGGCCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAACAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCACTCAGGAGCAAAA  >  NC_002947/778222‑778509
                                                                                                                                               |                                                                                                                                                   
cccGGGTGGCGGCCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCaa                                                                                                                                                  <  2:146700/147‑1 (MQ=255)
 ccGGGTGGCGGCCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGa                                                                                                                                               <  8:24850/149‑1 (MQ=255)
  cGGGTGGCGGCCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGAc                                                                                                                                              <  6:10321/149‑1 (MQ=255)
        gcggcCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGTCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑Tc                                                                                                                                     <  7:42641/149‑1 (MQ=255)
          ggcCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGa                                                                                                                                   <  2:253467/149‑1 (MQ=255)
            ccccGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACt                                                                                                                                 <  6:364422/149‑1 (MQ=255)
                   cTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACgg                                                                                                                          <  5:79160/149‑1 (MQ=255)
                   cTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACgg                                                                                                                          <  4:40180/149‑1 (MQ=255)
                                   cgctCCACTTTCACTTAGATGACACCCAGACCAGCCTCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGTGCCGCCAGCATGAGCTGATCTTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCAcc                                                                                                          <  1:210529/149‑1 (MQ=255)
                                   cgcgACACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCAcc                                                                                                          <  1:215888/149‑1 (MQ=255)
                                        cACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCg                                                                                                     <  2:164998/149‑1 (MQ=255)
                                        cACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCg                                                                                                     <  3:152293/149‑1 (MQ=255)
                                        cACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGAACGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCg                                                                                                     >  7:243718/1‑149 (MQ=255)
                                            ttCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTa                                                                                                  >  2:343015/1‑148 (MQ=255)
                                             tCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTa                                                                                                  <  1:348605/147‑1 (MQ=255)
                                                    gATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGCTT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGa                                                                                         >  4:41469/1‑149 (MQ=255)
                                                       gACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACa                                                                                       <  5:292895/148‑1 (MQ=255)
                                                         cacCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGa                                                                                    <  1:28133/149‑1 (MQ=255)
                                                          acCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGaa                                                                                   >  2:115664/1‑149 (MQ=255)
                                                               gACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAAcc                                                                              >  3:213903/1‑149 (MQ=255)
                                                                   aGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACAtcgtc                                                                                           >  7:39412/1‑132 (MQ=255)
                                                                   aGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACAtcgtc                                                                                           <  8:39412/132‑1 (MQ=255)
                                                                       tCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCg                                                                      >  3:159536/1‑149 (MQ=255)
                                                                            ggCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGa                                                                 <  8:243718/149‑1 (MQ=255)
                                                                                                 cAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGAt                                                                >  3:52013/1‑129 (MQ=255)
                                                                                                 cAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGAt                                                                <  4:52013/129‑1 (MQ=255)
                                                                                                    ctgctgGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCaa                                         >  3:173337/1‑149 (MQ=255)
                                                                                                     tgctgGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAAc                                        <  4:213903/149‑1 (MQ=255)
                                                                                                       ctgGGCCGCCAGCATTAGCTGATCTTGTGTGGCTAAATGAGCAATCACTATTTTCTCGACTATCAGCGGGGCCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCaa                                         <  1:115664/149‑1 (MQ=255)
                                                                                                       ctgGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGt                                      >  1:145895/1‑149 (MQ=255)
                                                                                                                cAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGTCAGCCCTCAACGTCAACCAGTg                             <  1:343015/149‑1 (MQ=255)
                                                                                                                cAGCATGAGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTg                             >  2:74193/1‑149 (MQ=255)
                                                                                                                       aGCTGATCGTGGGTGGCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATAc                      <  4:173337/149‑1 (MQ=255)
                                                                                                                               gtgggtggCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCa               <  1:74193/148‑1 (MQ=255)
                                                                                                                                tgggtggCAATTGGAACAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCACt             >  6:18581/1‑149 (MQ=255)
                                                                                                                                      gCAATTGGAGCAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCACTCAGGAg       >  6:266818/1‑149 (MQ=255)
                                                                                                                                            ggAACAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCACTCAGGAGCaaaa  <  5:18581/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCCGGGTGGCGGCCCCGACCTGTTCACCCTCAACCCGCGCCACTTTCACTTCGATGACACCCAGACCAGCATCGACGGCTATTTCACCGGCCCGTTCCAGCTGCTGGGCCGCCAGCATGAGCTGATCGTGGGTGGCAATTGGAACAAGGACGATT‑‑‑TCGACTATCACGGGGGGCGCGATGCCACCTACCGCTACATCGTCGACATGAACAACCTGGCAGCGTTCGATCCGCCCAAGCCGACAGCCCTCAACGTCAACCAGTGGCAATACAACCGCACTCAGGAGCAAAA  >  NC_002947/778222‑778509

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