Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 500,539 A→G 10.2% H363R (CAC→CGC)  PP_0411 → polyamine ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947500,5390AG10.2% 127.0 / 5.4 49H363R (CAC→CGC) PP_0411polyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (3/2);  total (25/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 = 6.87e-01

CGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGC  >  NC_002947/500396‑500684
                                                                                                                                               |                                                                                                                                                  
cGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTCGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGccc                                                                                                                                               >  5:233118/1‑149 (MQ=255)
 gAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCcg                                                                                                                                              >  7:300306/1‑149 (MQ=255)
       aTTTACGTGGGCTAGCACGTGCGGGTGCGTCTGGAGGTTTGCGGCATAGGCGACTTCTTCGTGAAGCAGCGGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGAGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCt                                                                                                                                        <  8:199622/149‑1 (MQ=255)
            ccTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATc                                                                                                                                  >  2:134126/1‑149 (MQ=255)
                ggCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTTGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGa                                                                                                                               <  1:21315/148‑1 (MQ=255)
                       aCGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCg                                                                                                                         >  3:344279/1‑147 (MQ=255)
                                tGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGa                                                                                                              <  1:134126/149‑1 (MQ=255)
                                 gCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGAt                                                                                                             >  8:261084/1‑149 (MQ=255)
                                        gAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGc                                                                                                      >  4:96956/1‑149 (MQ=255)
                                          ggTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCtt                                                                                                    <  3:96956/149‑1 (MQ=255)
                                           gTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCtt                                                                                                    >  3:81487/1‑148 (MQ=255)
                                             ttGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGc                                                                                                                     <  5:99313/129‑1 (MQ=255)
                                             ttGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGc                                                                                                                     >  6:99313/1‑129 (MQ=255)
                                               gCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCAcc                                                                                                <  8:300306/148‑1 (MQ=255)
                                                   cAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAAccc                                                                                           >  8:22147/1‑149 (MQ=255)
                                                        gCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTAc                                                                                      <  7:261084/149‑1 (MQ=255)
                                                           aCTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACtg                                                                                    <  1:72790/149‑1 (MQ=255)
                                                                   gTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgagaga                                                                            <  6:101843/148‑1 (MQ=255)
                                                                   gTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaa                                                                           >  7:268923/1‑149 (MQ=255)
                                                                    tGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaat                                                                          <  8:268923/149‑1 (MQ=255)
                                                                         cagcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACtgt                                                                                   <  2:156110/135‑1 (MQ=255)
                                                                         cagcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACtgt                                                                                   >  1:156110/1‑135 (MQ=255)
                                                                           gcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATg                                                                   <  4:81487/149‑1 (MQ=255)
                                                                              gATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataataa                                                                     <  7:129616/144‑1 (MQ=255)
                                                                              gATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataataa                                                                     >  8:129616/1‑144 (MQ=255)
                                                                                ttGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATgcgc                                                                >  5:23649/1‑147 (MQ=255)
                                                                                     gAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataa                                                                        <  1:242154/134‑1 (MQ=255)
                                                                                     gAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataa                                                                        >  2:242154/1‑134 (MQ=255)
                                                                                              ccGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGa                                                                                                                               >  5:49035/1‑70 (MQ=255)
                                                                                              ccGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGa                                                                                                                               <  6:49035/70‑1 (MQ=255)
                                                                                               cGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGAc                                               >  6:144606/1‑149 (MQ=255)
                                                                                                       gCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTTGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTg                                        <  7:22147/148‑1 (MQ=255)
                                                                                                           tGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgag                                                                               <  8:34598/105‑1 (MQ=255)
                                                                                                           tGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgag                                                                               >  7:34598/1‑105 (MQ=255)
                                                                                                             ggCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                  <  6:23649/148‑1 (MQ=255)
                                                                                                              gCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTggg                                >  1:206569/1‑149 (MQ=255)
                                                                                                                    tagtaCCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGtt                          <  2:206569/149‑1 (MQ=255)
                                                                                                                       taCCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCAc                       >  2:52962/1‑149 (MQ=255)
                                                                                                                          cGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGc                    >  8:198282/1‑149 (MQ=255)
                                                                                                                               ggctgggAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTggc                <  6:233118/148‑1 (MQ=255)
                                                                                                                                  tgggAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTggccatg            >  8:178392/1‑147 (MQ=255)
                                                                                                                                        gTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATgg      <  4:344279/149‑1 (MQ=255)
                                                                                                                                        atcgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTgggg                               <  8:88294/121‑1 (MQ=255)
                                                                                                                                        atcgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTgggg                               >  7:88294/4‑124 (MQ=255)
                                                                                                                                        atcgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                  <  4:43470/118‑1 (MQ=255)
                                                                                                                                        atcgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                  >  3:43470/4‑121 (MQ=255)
                                                                                                                                          cgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGACCAGGCCATGGcc    >  1:119069/2‑149 (MQ=255)
                                                                                                                                            aGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAtggcctggccag  <  7:178392/149‑3 (MQ=255)
                                                                                                                                             gCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCg   >  2:169941/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGC  >  NC_002947/500396‑500684

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: