Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,994,718 G→A 7.8% L463L (TTG→TTA)  PP_5251 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,994,7180GA 7.8% 145.0 / 3.3 51L463L (TTG→TTA) PP_5251transporter
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GGTGCCGGTGTTGAGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCA  >  NC_002947/5994574‑5994861
                                                                                                                                                |                                                                                                                                               
ggtgccggtgTTGAGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTg                                                                                                                                            <  1:652096/150‑1 (MQ=255)
   gccggtgTTGAGGGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTgggg                                                                                                                                         <  2:277980/150‑1 (MQ=255)
       gtgtTGAGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGt                                                                                                                                       >  2:40783/1‑148 (MQ=255)
       gtgtTGAGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGt                                                                                                                                       <  1:40783/148‑1 (MQ=255)
             aGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCgg                                                                                                                               <  2:1135844/150‑1 (MQ=255)
             aGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCgg                                                                                                                               <  1:440429/150‑1 (MQ=255)
             aGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCgg                                                                                                                               <  2:839456/150‑1 (MQ=255)
                    aCATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGGCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTg                                                                                                                        <  2:765810/150‑1 (MQ=255)
                     cATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGc                                                                                                                       <  2:537896/150‑1 (MQ=255)
                     cATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCCCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATCGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGc                                                                                                                       <  2:1030525/150‑1 (MQ=255)
                      aTGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCt                                                                                                                      <  1:820465/150‑1 (MQ=255)
                      aTGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCACCATGGCGCACTGTCGTGGGGGTTTATGCGGCGCAGTGCt                                                                                                                      >  2:69630/1‑150 (MQ=255)
                          aCAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATc                                                                                                                  >  2:905940/1‑150 (MQ=255)
                              ggAAGCCGGACCGGTGATGGTCACCGCCGCATCCCATATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGc                                                                                                              <  2:96353/150‑1 (MQ=255)
                               gAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCg                                                                                                             >  2:516261/1‑150 (MQ=255)
                                   ccGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCCAATCCGGCGTtgc                                                                                                         <  2:923556/150‑2 (MQ=255)
                                       aCCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCt                                                                                                     <  2:636212/150‑1 (MQ=255)
                                                ggTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGttctt                                                                                            >  1:1170853/1‑150 (MQ=255)
                                                 gTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGttcttc                                                                                           >  1:1041291/1‑150 (MQ=255)
                                                                 acaTCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGc                                                                           <  2:1133629/150‑1 (MQ=255)
                                                                    tCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTgg                                                                        >  1:413470/1‑150 (MQ=255)
                                                                     ctcccccACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGa                                                                       <  1:905940/148‑1 (MQ=255)
                                                                       cccccACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGc                                                                     >  2:941807/1‑150 (MQ=255)
                                                                           caccaAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCt                                                                 >  1:490495/1‑150 (MQ=255)
                                                                                aGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGAc                                                            >  1:1116998/1‑150 (MQ=255)
                                                                                 ggCCGAGTCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTATCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGa                                                                       >  2:621597/1‑138 (MQ=255)
                                                                                 ggCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTATCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGa                                                                       <  1:621686/138‑1 (MQ=255)
                                                                                 ggCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTATCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGa                                                                       >  2:621686/1‑138 (MQ=255)
                                                                                 ggCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTATCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGa                                                                       <  1:621597/138‑1 (MQ=255)
                                                                                        gCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGGCGTGGGGGTTGATGCGGGACAGTGCTGATCCTGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATcaccac                                                    >  2:581903/1‑150 (MQ=255)
                                                                                           ttCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTg                                                                                                    <  1:1103555/99‑1 (MQ=255)
                                                                                           ttCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTg                                                                                                    >  2:1103555/1‑99 (MQ=255)
                                                                                           ttCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCgg                                                 <  2:413470/150‑1 (MQ=255)
                                                                                             ccAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGt                                               <  1:425819/150‑1 (MQ=255)
                                                                                               agcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGa                                             <  1:516261/150‑1 (MQ=255)
                                                                                                      ggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCgtgg                                      >  2:724308/1‑150 (MQ=255)
                                                                                                           gggAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAAc                                 <  1:724308/150‑1 (MQ=255)
                                                                                                              aGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGa                              <  2:1116998/150‑1 (MQ=255)
                                                                                                               gCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGaa                             >  2:211317/1‑150 (MQ=255)
                                                                                                                  ggAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGAt                          >  2:906164/1‑150 (MQ=255)
                                                                                                                  ggAAGCGATGCGCAAGCGCAATGGCGCATTGGCGAGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGAt                          >  2:741404/1‑150 (MQ=255)
                                                                                                                      gCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAg                      <  2:490495/150‑1 (MQ=255)
                                                                                                                            cgcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCgg                <  1:1112863/150‑1 (MQ=255)
                                                                                                                             gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGGCAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGCCGTCACCACCTGGTGACCCGTGGTGAACTGAAGATCGAGGCGCggg               >  2:264246/1‑150 (MQ=255)
                                                                                                                             gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCgtgg                                      >  2:352719/1‑127 (MQ=255)
                                                                                                                             gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCgtgg                                      <  1:352719/127‑1 (MQ=255)
                                                                                                                                agcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCa            <  1:941807/150‑1 (MQ=255)
                                                                                                                                  cgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAgg          <  1:211317/150‑1 (MQ=255)
                                                                                                                                   gcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTAGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGa         >  2:976994/1‑150 (MQ=255)
                                                                                                                                     aaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCGCCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtg       >  1:660017/1‑150 (MQ=255)
                                                                                                                                          cgcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGAATGGACGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCa  <  1:317355/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGTGCCGGTGTTGAGCGATGACATGGACAAGGAAGCCGGACCGGTGATGGTCACCGTCGAGTACCACATCGCCCCCACCAAGGCCGAGGCCTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCCGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCA  >  NC_002947/5994574‑5994861

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