Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 427,755 A→C 10.2% L16R (CTG→CGG) ‡ PP_0350 ← ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947427,7550AC10.2% 117.2 / 7.7 49L16R (CTG→CGG) ‡PP_0350ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base C (3/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGAT  >  NC_002947/427610‑427904
                                                                                                                                                 |                                                                                                                                                     
gATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGc                                                                                                                                                   <  1:774123/150‑1 (MQ=255)
 aTGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGcc                                                                                                                                                  >  2:726778/1‑150 (MQ=255)
  tGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCa                                                                                                                                                 <  1:90325/150‑1 (MQ=255)
       agggCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAgg                                                                                                                                            >  2:1041671/3‑150 (MQ=255)
         ggCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAgcgc                                                                                                                                          <  1:44765/150‑4 (MQ=255)
          gCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGCCGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTg                                                                                                                                         <  1:726778/150‑1 (MQ=255)
           cAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGCCTGATGGCCCTGGCCATcggcctgg                                                                                                                                        <  1:317522/150‑8 (MQ=255)
              tGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGgc                                                                                                                                     <  1:857443/150‑1 (MQ=255)
                 tggCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCa                                                                                                                                  <  1:701899/149‑1 (MQ=255)
                   gggtggcgggcGAGCTGGGTGAGGGGGCGGGGCAGTGGAGCGCTGGGCAGGTTGTGGTGGCGCGGCGTGGATTGCTGGGCGGCGCACGAGGGAGGGTTGGTCAGCGGGGCGGCCAGGCCGAGGGCCAGGGCCATCAGGCTGGGGCGCAgg                                                                                                                                <  2:67095/147‑1 (MQ=255)
                                   ggTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCAc                                                                                                                <  1:870448/150‑1 (MQ=255)
                                     ttGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCCCGGCGGCGAACTAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGGTCAACGGCACgg                                                                                                              >  1:897585/1‑150 (MQ=255)
                                            tggtggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGa                                                                                                                          <  1:578541/131‑1 (MQ=255)
                                            tggtggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGa                                                                                                                          >  2:578541/1‑131 (MQ=255)
                                                ggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATgcggcg                                                                                                   >  1:1183501/1‑150 (MQ=255)
                                                 gCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGcggcgg                                                                                                  <  1:642197/150‑1 (MQ=255)
                                                    aGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTc                                                                                               <  2:1269304/150‑1 (MQ=255)
                                                                  cAGGGTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAgg                                                                                 >  1:418657/1‑149 (MQ=255)
                                                                     gTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgct                                                                              >  1:346474/1‑150 (MQ=255)
                                                                      ttGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgctt                                                                             <  2:346474/150‑1 (MQ=255)
                                                                                        gATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                                           >  1:326472/1‑150 (MQ=255)
                                                                                              tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACg                                                                                        >  1:91997/1‑115 (MQ=255)
                                                                                              tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACg                                                                                        <  2:91997/115‑1 (MQ=255)
                                                                                                ggcggcGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCt                                                    <  1:407493/149‑1 (MQ=255)
                                                                                                ggcggcGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGAGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCt                                                    >  2:407493/1‑149 (MQ=255)
                                                                                                   ggcgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaa                                                >  2:405367/1‑150 (MQ=255)
                                                                                                     cgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAg                                              <  2:930946/150‑1 (MQ=255)
                                                                                                      gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGt                                             >  1:85807/1‑150 (MQ=255)
                                                                                                      gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGt                                             >  1:180154/1‑150 (MQ=255)
                                                                                                        aCAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGa                                           >  2:429851/1‑150 (MQ=255)
                                                                                                           aGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGAGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAg                                        <  1:1113755/150‑1 (MQ=255)
                                                                                                            ggCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAgg                                       >  2:457501/1‑150 (MQ=255)
                                                                                                             gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAggg                                      >  2:271030/1‑150 (MQ=255)
                                                                                                                         aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGaaaa                          >  1:31853/1‑150 (MQ=255)
                                                                                                                         aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATAAGGACGGGGAGTCTGCTACCGTCGTCGTCACTGGACATGCCGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGaaaa                          >  1:24757/1‑150 (MQ=255)
                                                                                                                              ggCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                     <  2:85807/150‑1 (MQ=255)
                                                                                                                              ggCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                     <  2:31853/150‑1 (MQ=255)
                                                                                                                              ggCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                     <  2:326472/150‑1 (MQ=255)
                                                                                                                              ggCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                     <  1:429851/150‑1 (MQ=255)
                                                                                                                                   ccagcctgATGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGaa                >  2:1064382/8‑150 (MQ=255)
                                                                                                                                       cctgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAatat                      <  1:696095/137‑1 (MQ=255)
                                                                                                                                       cctgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAatat                      >  2:696095/4‑140 (MQ=255)
                                                                                                                                         tgATGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTAttt          >  1:535431/2‑150 (MQ=255)
                                                                                                                                          gATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTa         >  2:623709/1‑150 (MQ=255)
                                                                                                                                           aTGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAatta             >  2:4590/1‑145 (MQ=255)
                                                                                                                                           aTGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAatta             <  1:4590/145‑1 (MQ=255)
                                                                                                                                             ggCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGTCAGGGTGAATCTGAAAAATATTGAGAATTATTTAATg      <  1:1041671/150‑1 (MQ=255)
                                                                                                                                              gCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATgg     <  1:457501/150‑1 (MQ=255)
                                                                                                                                                cAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGa   >  2:921932/1‑150 (MQ=255)
                                                                                                                                                 aGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGAt  >  1:1137349/1‑150 (MQ=255)
                                                                                                                                                 aGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGAt  <  2:1137349/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGAT  >  NC_002947/427610‑427904

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