Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,688,589 G→A 8.7% P330P (CCC→CCT)  ydcT ← polyamine ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,688,5890GA 8.7% 125.0 / 3.4 46P330P (CCC→CCT) ydcTpolyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (2/2);  total (21/25)
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.16e-01

GGTAAATCGCGCCGAACCAGGTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTGCCTTGGTACTGGACATCGTGCA  >  NC_002947/1688443‑1688737
                                                                                                                                                  |                                                                                                                                                    
ggTAAATCGCGCCGAACCAGGTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGc                                                                                                                                                   >  1:177019/1‑150 (MQ=255)
          gcCGAACCAGGTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCt                                                                                                                                         <  1:1009372/150‑1 (MQ=255)
               aCCAGGTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAgg                                                                                                                                    <  1:1121388/150‑1 (MQ=255)
                   ggTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGtt                                                                                                                                <  1:841263/150‑1 (MQ=255)
                   ggTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGCCCTCTCGGGGCCAGTGCGCCTGCAGGCGtt                                                                                                                                <  1:477235/150‑1 (MQ=255)
                     tCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGc                                                                                                                              <  1:807175/150‑1 (MQ=255)
                         cggcggGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGa                                                                                                                          <  1:572484/150‑1 (MQ=255)
                         cggcggGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGa                                                                                                                          <  2:826035/150‑1 (MQ=255)
                                   aaCAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGc                                                                                                                <  1:554583/150‑1 (MQ=255)
                                    aCAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGcc                                                                                                               >  2:745786/1‑150 (MQ=255)
                                    aCAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGcc                                                                                                               >  1:540849/1‑150 (MQ=255)
                                     cagcagTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCCCGCGAGGCCatga                                                                                                                                              >  2:1038647/1‑115 (MQ=255)
                                     cagcagTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCCCGCGAGGCCatga                                                                                                                                              <  1:1038647/118‑4 (MQ=255)
                                       gcagTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCCCCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGc                                                                                                            >  1:242228/1‑150 (MQ=255)
                                       gcagTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGc                                                                                                            <  2:177019/150‑1 (MQ=255)
                                             ccAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGc                                                                                                      >  1:178567/1‑150 (MQ=255)
                                              cAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCt                                                                                                     <  2:964407/150‑1 (MQ=255)
                                               aGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTACTGCAGCCCCGTCCTGGCCTCGCGGGGCCAGGGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCtt                                                                                                    >  1:117201/1‑150 (MQ=255)
                                                       gTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGc                                                                                            >  1:606098/1‑150 (MQ=255)
                                                        ttGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCg                                                                                           <  2:540849/150‑1 (MQ=255)
                                                                  ggTACAGGAGGTTCGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGCCATGGCCTCGCGGGGCCAGTGCGACTGCAGGCGTTGCCCGAGCTGCCAGGCCCGCATCTGCTTCTGCCAGCGGTCATTGGcc                                                                                 <  2:367815/150‑1 (MQ=255)
                                                                           ggTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACt                                                                        >  2:600367/1‑150 (MQ=255)
                                                                               ggACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCca                                                                    <  2:882158/150‑1 (MQ=255)
                                                                                 aCAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCcagc                                                                  >  1:404997/1‑150 (MQ=255)
                                                                                   aGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCcagcag                                                                <  2:178567/150‑1 (MQ=255)
                                                                                             gtgtGCTCATGTCAGCGGCCCTCGGTTTCCTGCAACACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCACGGGGGCTCAGGCCCGCTTCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGCCGACCGCTg                                                      >  1:78861/1‑150 (MQ=255)
                                                                                                  cTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTc                                                                                        >  2:386132/1‑111 (MQ=255)
                                                                                                  cTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTc                                                                                        <  1:386132/111‑1 (MQ=255)
                                                                                                   tCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGc                                                                                                      <  1:576823/96‑1 (MQ=255)
                                                                                                   tCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGc                                                                                                      >  2:576823/1‑96 (MQ=255)
                                                                                                   tCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGc                                                >  2:637088/1‑150 (MQ=255)
                                                                                                    cATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCt                                               <  1:637088/150‑1 (MQ=255)
                                                                                                           gCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTc                                        >  2:887080/1‑150 (MQ=255)
                                                                                                           gCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTc                                        <  1:887080/150‑1 (MQ=255)
                                                                                                                 cTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCg                                  >  1:899240/1‑150 (MQ=255)
                                                                                                                  tCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCgg                                 >  1:377122/1‑150 (MQ=255)
                                                                                                                  tCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCgg                                 <  1:858180/150‑1 (MQ=255)
                                                                                                                  tCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCgg                                 <  2:606098/150‑1 (MQ=255)
                                                                                                                    ggTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGt                               <  2:404997/150‑1 (MQ=255)
                                                                                                                      tttCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTgg                             <  1:745786/150‑1 (MQ=255)
                                                                                                                           tGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTg                        >  1:927950/1‑150 (MQ=255)
                                                                                                                           tGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTg                        >  1:997072/1‑150 (MQ=255)
                                                                                                                           tGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTg                        >  2:89648/1‑150 (MQ=255)
                                                                                                                                          gccccgcgAGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTg                                                                       <  1:558506/84‑1 (MQ=255)
                                                                                                                                          gCCCCGCGAGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTg                                                                       >  2:558506/1‑88 (MQ=255)
                                                                                                                                            cTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTGCCTTGGTACTGGACATc       <  2:899240/150‑1 (MQ=255)
                                                                                                                                                 ggggCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGCAGCGGGTGGCGCTGCCTTGGTACTGGACATCGTgca  <  2:1042626/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTAAATCGCGCCGAACCAGGTCAGCGGCGGGATCAACAGCAGTGCCAGGTACAGGTTGGGGCGCCGGTACAGCAGGTTGGACAGGGTGCGCAGTGTGCTCATGTCAGCGGCCCTCGGTTTCCTGCAGCACCGTCATGGCCTCGCGGGGCCAGTGCGCCTGCAGGCGTTGCCCGAGCTGCCAGGCCTGCATCTGCTTCTGCCAGCGGTCATTGGCCTGGCTTACTGCCAGCAGTTGACCGTTGTCCAGCTGCAGTTCGTAGCGGGTGGCGCTGCCTTGGTACTGGACATCGTGCA  >  NC_002947/1688443‑1688737

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