Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 427,752 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,7520AC10.2% 171.4 / 11.2 59L16R (CTG→CGG) ‡PP_0350ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base C (4/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATT  >  NC_002947/427607‑427892
                                                                                                                                                 |                                                                                                                                            
gATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAggg                                                                                                                                           <  2:217244/149‑1 (MQ=255)
           cAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTg                                                                                                                                <  8:159518/149‑1 (MQ=255)
           cAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTg                                                                                                                                <  6:75415/149‑1 (MQ=255)
            aaTGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCGGGCTgg                                                                                                                               <  6:176141/149‑1 (MQ=255)
             aTGCCGGCCGGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTggg                                                                                                                              <  8:226907/149‑1 (MQ=255)
              tGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCAAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTgggg                                                                                                                             >  1:192649/1‑149 (MQ=255)
               gCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGgc                                                                                                                            >  4:158197/1‑149 (MQ=255)
                  ggCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGgcg                                                                                                                           >  5:332864/1‑147 (MQ=255)
                  ggCCGGGCTGGCGAGCGGGTGGAGGGTGGTGGCCAGGGGCGCGCGGGGCAGGGGGGAGGCGCGCGACGTGGATTGCGCGGCGGCGACCACGGCCGGGCTGGGCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCGGGGGCGCa                                                                                                                         <  3:210830/149‑1 (MQ=255)
                   gCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAg                                                                                                                        >  3:176456/1‑149 (MQ=255)
                        gctggcGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGata                                                                                                                    <  8:294005/148‑1 (MQ=255)
                                    gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCAc                                                                                                       >  6:226531/1‑149 (MQ=255)
                                     ttGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTggg                                                                                                                              >  4:141454/1‑125 (MQ=255)
                                     ttGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTggg                                                                                                                              <  3:141454/125‑1 (MQ=255)
                                     ttGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACg                                                                                                      >  2:343622/1‑149 (MQ=255)
                                          ggtggtGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAg                                                                                                                        >  7:153689/1‑126 (MQ=255)
                                          ggtggtGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAg                                                                                                                        <  8:153689/126‑1 (MQ=255)
                                          ggtggtGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGAc                                                                                                  <  2:160673/148‑1 (MQ=255)
                                                  ccAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGcggcgg                                                                                         >  1:123187/1‑149 (MQ=255)
                                                          gcgcTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGa                                                                                 >  4:115792/1‑149 (MQ=255)
                                                                     gTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgc                                                                      >  3:286588/1‑149 (MQ=255)
                                                                           ggCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCt                                                                                    <  5:349004/129‑1 (MQ=255)
                                                                           ggCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCt                                                                                    >  6:349004/1‑129 (MQ=255)
                                                                           ggCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTtgctg                                                                >  1:244314/1‑149 (MQ=255)
                                                                               cgcgACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTc                                                            <  8:93013/149‑1 (MQ=255)
                                                                               cgcgACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTc                                                            <  3:156977/149‑1 (MQ=255)
                                                                               cgcgACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTc                                                            <  3:115792/149‑1 (MQ=255)
                                                                                      tGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcg                                                     <  5:226531/149‑1 (MQ=255)
                                                                                        gATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgag                                                   >  3:42423/1‑149 (MQ=255)
                                                                                          ttGCTCGGCGGCGAACAAGGCGGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                                  >  2:181505/1‑148 (MQ=255)
                                                                                              tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATc                                                                                                               <  3:229314/83‑1 (MQ=255)
                                                                                              tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATc                                                                                                               >  4:229314/1‑83 (MQ=255)
                                                                                                 gcggcgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                         >  1:157070/1‑134 (MQ=255)
                                                                                                 gcggcgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                         <  2:157070/134‑1 (MQ=255)
                                                                                                         cAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGAACAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGa                                  >  1:308359/1‑149 (MQ=255)
                                                                                                             gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttgc                                                                  >  2:205798/1‑113 (MQ=255)
                                                                                                             gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttgc                                                                  <  1:205798/113‑1 (MQ=255)
                                                                                                             gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAgg                              <  2:192649/149‑1 (MQ=255)
                                                                                                             gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAgg                              <  2:244314/149‑1 (MQ=255)
                                                                                                               aGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGt                            <  1:343622/149‑1 (MQ=255)
                                                                                                                 ggCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGc                                                       <  8:130295/120‑1 (MQ=255)
                                                                                                                 ggCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGc                                                       >  7:130295/1‑120 (MQ=255)
                                                                                                                 ggCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGa                          <  4:286588/149‑1 (MQ=255)
                                                                                                                   cTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAAt                        >  1:170031/1‑149 (MQ=255)
                                                                                                                       tCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTg                     >  2:284717/1‑148 (MQ=255)
                                                                                                                         aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACgg                                                                              <  2:185511/89‑1 (MQ=255)
                                                                                                                         aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACgg                                                                              >  1:185511/1‑89 (MQ=255)
                                                                                                                         aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGaa                   <  2:170031/148‑1 (MQ=255)
                                                                                                                                   ccAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTgag         >  1:25362/1‑148 (MQ=255)
                                                                                                                                       gCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAatta    <  3:158197/149‑1 (MQ=255)
                                                                                                                                        ctgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttgc                                                                  >  7:362735/3‑86 (MQ=255)
                                                                                                                                        ctgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                                  <  6:34720/100‑1 (MQ=255)
                                                                                                                                        ctgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                                  >  5:34720/3‑102 (MQ=255)
                                                                                                                                        ccGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattat   >  5:304682/1‑149 (MQ=255)
                                                                                                                                         tgATGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattatt  >  5:186489/2‑149 (MQ=255)
                                                                                                                                         tgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagcg                                                 >  1:248126/2‑102 (MQ=255)
                                                                                                                                         tgATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagcg                                                 <  2:248126/101‑1 (MQ=255)
                                                                                                                                         cGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattatt  <  2:123187/149‑1 (MQ=255)
                                                                                                                                         cGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattatt  <  4:176456/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATT  >  NC_002947/427607‑427892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑