Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 500,548 C→T 10.2% A366V (GCG→GTG)  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,5480CT10.2% 133.2 / 4.7 49A366V (GCG→GTG) PP_0411polyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (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 = 9.95e-01

ATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAC  >  NC_002947/500403‑500688
                                                                                                                                                 |                                                                                                                                              
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                                                                                       >  1:156110/1‑135 (MQ=255)
                                                                  cagcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACtgt                                                                                       <  2:156110/135‑1 (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                                                                            >  2:242154/1‑134 (MQ=255)
                                                                              gAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataa                                                                            <  1:242154/134‑1 (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                                      >  3:43470/4‑121 (MQ=255)
                                                                                                                                 atcgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                      <  4:43470/118‑1 (MQ=255)
                                                                                                                                   cgAGCGCGCGGGCGTGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGACCAGGCCATGGcc        >  1:119069/2‑149 (MQ=255)
                                                                                                                                     aGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAtggcctggccag      <  7:178392/149‑3 (MQ=255)
                                                                                                                                      gCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCg       >  2:169941/1‑147 (MQ=255)
                                                                                                                                         cgcgggcgtgcTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAc  <  2:119069/144‑1 (MQ=255)
                                                                                                                                         cGCCCGCGCGCTCGATCCGATTGCCGATAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGa   >  1:193045/1‑149 (MQ=255)
                                                                                                                                         cGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAc  <  1:169941/149‑1 (MQ=255)
                                                                                                                                             cgcgcgCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCAcc                                                                                                    <  4:151293/47‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCGATTGCCGAAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAC  >  NC_002947/500403‑500688

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: