Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,277,598 T→C R433R (CGT→CGC yejF → putative oligopeptide ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,277,5980TC98.4% 210.3 / ‑5.0 62R433R (CGT→CGCyejFputative oligopeptide ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGCTTACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGTATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGCCT  >  NC_000913/2277461‑2277734
                                                                                                                                         |                                                                                                                                        
ggCTTACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCAt                                                                                                                                        >  2:139889/1‑140 (MQ=255)
  cTTACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTg                                                                                                                                      >  2:552336/1‑140 (MQ=255)
  cTTACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTg                                                                                                                                      >  2:567700/1‑140 (MQ=255)
   ttACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGc                                                                                                                                     <  1:103630/140‑1 (MQ=255)
              ccAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAggg                                                                                                                          <  1:69617/140‑1 (MQ=255)
                 gCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCAt                                                                                                                       >  2:611973/1‑140 (MQ=255)
                 gCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCAt                                                                                                                       <  2:953827/140‑1 (MQ=255)
                  ccGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCAtt                                                                                                                      <  2:1150962/140‑1 (MQ=255)
                   cGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCg                                                                                                                                    >  2:382046/1‑125 (MQ=255)
                   cGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCg                                                                                                                                    <  1:382046/125‑1 (MQ=255)
                   cGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGcc                                                                                                                              <  1:772948/131‑1 (MQ=255)
                   cGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGcc                                                                                                                              >  2:772948/1‑131 (MQ=255)
                    gACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTaa                                                                                                                    <  1:106347/140‑1 (MQ=255)
                            cTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAg                                                                                                            <  1:156830/140‑1 (MQ=255)
                             tGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGc                                                                                                           <  2:151465/140‑1 (MQ=255)
                              gCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGcc                                                                                                          <  1:929255/140‑1 (MQ=255)
                                  cacaACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCg                                                                                                      <  1:846635/140‑1 (MQ=255)
                                  cacaACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCg                                                                                                      <  1:1067383/140‑1 (MQ=255)
                                      aCGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGa                                                                                                  <  2:194245/140‑1 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTaa                                                                                                             <  1:369305/128‑1 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTaa                                                                                                             >  2:369305/1‑128 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGccc                                                                                                         <  1:30733/132‑1 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGccc                                                                                                         >  2:30733/1‑132 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGAt                                                                                                 <  2:651398/140‑1 (MQ=255)
                                       cgcgAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGAt                                                                                                 >  2:968248/1‑140 (MQ=255)
                                                 aaGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTtgat                                                                                       >  1:27276/1‑140 (MQ=255)
                                                 aaGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTtgat                                                                                       <  1:1055554/140‑1 (MQ=255)
                                                    tGATAGCCGTGATGCATGAAGTGGGATTGGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGaa                                                                                    >  1:1050536/1‑140 (MQ=255)
                                                    tGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGaa                                                                                    <  2:27276/140‑1 (MQ=255)
                                                          ccGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGAca                                                                              >  1:522083/1‑140 (MQ=255)
                                                             tGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGAcatca                                                                           >  1:1063392/1‑140 (MQ=255)
                                                               aTGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACatcatc                                                                         >  2:1004371/1‑140 (MQ=255)
                                                                     gAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGa                                                                   <  1:1071867/140‑1 (MQ=255)
                                                                      aaGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGAc                                                                  >  2:779476/1‑140 (MQ=255)
                                                                           ggATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAAc                                                             >  1:1092631/1‑140 (MQ=255)
                                                                            gATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGTATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACa                                                                 <  1:419845/135‑1 (MQ=255)
                                                                            gATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGAATTGCGATTGCCAGGGCATTAATTCTTAAGCCTTCGCTGATCATACTTGATGAACCGACATCATCACTCGACa                                                                 >  2:419845/1‑135 (MQ=255)
                                                                               tAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTa                                                         <  1:567700/140‑1 (MQ=255)
                                                                                aGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCGTCACTCGACAAAACGGTAc                                                        >  2:457758/1‑140 (MQ=255)
                                                                                 gATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACa                                                       <  1:611973/140‑1 (MQ=255)
                                                                                  aTCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAg                                                      >  1:248661/1‑140 (MQ=255)
                                                                                  aTCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAg                                                      <  2:1063392/140‑1 (MQ=255)
                                                                                  aTCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAg                                                      >  2:625547/1‑140 (MQ=255)
                                                                                   tCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAgg                                                     >  1:328298/1‑140 (MQ=255)
                                                                                           cacGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAata                                             >  1:615694/1‑140 (MQ=255)
                                                                                             cGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATAtt                                           >  2:326838/1‑140 (MQ=255)
                                                                                             cGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATAtt                                           >  1:153136/1‑140 (MQ=255)
                                                                                             cGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATAtt                                           >  2:102361/1‑140 (MQ=255)
                                                                                               ccACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGa                                         <  2:328298/140‑1 (MQ=255)
                                                                                               ccACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGa                                         >  1:795573/1‑140 (MQ=255)
                                                                                                 aCCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACg                                       <  1:326838/140‑1 (MQ=255)
                                                                                                        tCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGa                                >  1:370467/1‑140 (MQ=255)
                                                                                                            gcggAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATc                            >  2:651650/1‑140 (MQ=255)
                                                                                                              ggAGTTCTCTGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCAt                          <  1:625547/140‑1 (MQ=255)
                                                                                                                     tctGGTGGTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGcaaca                   <  2:1092631/140‑1 (MQ=255)
                                                                                                                       tggtggTCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACaaa                 >  1:635380/1‑140 (MQ=255)
                                                                                                                            gtCAGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCAt            <  2:635380/140‑1 (MQ=255)
                                                                                                                               aGCGACAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCaa         >  1:139286/1‑140 (MQ=255)
                                                                                                                                    cAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGc    >  2:1114067/1‑140 (MQ=255)
                                                                                                                                    cAACGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGc    <  1:1004371/140‑1 (MQ=255)
                                                                                                                                      aCGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGCCt  <  1:139889/140‑1 (MQ=255)
                                                                                                                                      aCGCATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGCCt  <  1:651650/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GGCTTACGGGTTCACCAGCCGACGCTTTCTGCCGCACAACGCGAACAACAAGTGATAGCCGTGATGCATGAAGTGGGATTAGATCCTGAAACACGCCACCGTTATCCGGCGGAGTTCTCTGGTGGTCAGCGACAACGTATTGCGATTGCCAGGGCATTAATTCTTAAGCCCTCGCTGATCATACTTGATGAACCGACATCATCACTCGACAAAACGGTACAGGCGCAAATATTGACGCTATTGAAATCATTGCAACAAAAGCATCAACTGGCCT  >  NC_000913/2277461‑2277734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: