Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 856,822 (GC)5→4 coding (860‑861/1593 nt) ybiT → ABC transporter ATP‑binding protein YbiT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913856,8130C.93.3% 111.2 / 2.2 30coding (851/1593 nt)ybiTABC transporter ATP‑binding protein YbiT
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (10/18);  total (11/19)
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.58e-01
*NC_000913856,8140G.93.3% 111.2 / 2.1 30coding (852/1593 nt)ybiTABC transporter ATP‑binding protein YbiT
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (10/18);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCAAAGG  >  NC_000913/856679‑856945
                                                                                                                                      ||                                                                                                                                   
gTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcg                                                                                                                                 >  1:133458/1‑140 (MQ=255)
  aCATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcgcg                                                                                                                               <  2:462293/140‑1 (MQ=255)
  aCATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcgcg                                                                                                                               <  2:621739/140‑1 (MQ=255)
  aCATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTcgcgcgcg                                                                                                                               >  1:431248/1‑140 (MQ=255)
         ggcggcGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGACAGCGCAGATTGCTGAGGTGCAATCTTCGTTAGCCGCTTTAGCGCCAACGTCTCGAAATCTCGCGAGGCACCTTagcccgcgccatattg                                                                                                                       >  2:469045/1‑133 (MQ=255)
            ggcgACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGAt                                                                                                                   <  1:402739/140‑1 (MQ=255)
              cgACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCGGATTGATaa                                                                                                                 >  2:379727/1‑140 (MQ=255)
              cgACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATaa                                                                                                                 >  2:387899/1‑140 (MQ=255)
                    aGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCaa                                                                                                           <  2:463005/140‑1 (MQ=255)
                         cgTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCTGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTgg                                                                                                      <  2:387305/140‑1 (MQ=255)
                          gTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGa                                                                                                     <  1:144737/140‑1 (MQ=255)
                           tGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGaa                                                                                                    >  2:121650/1‑140 (MQ=255)
                                           gATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAg                                                                                   <  1:469045/140‑1 (MQ=255)
                                            aTAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAg                                                                                   >  1:92820/1‑140 (MQ=255)
                                                  ccAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACAT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCa                                                                             >  2:641251/1‑140 (MQ=255)
                                                   cAAGAAGAAAGCGCACATTGCTCAGCTGCCATCTTTCTTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAg                                                                            <  1:181373/140‑1 (MQ=255)
                                                           aaGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTGAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGtt                                                                    >  2:420359/1‑140 (MQ=255)
                                                                    ttGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGttt                                                           >  1:158515/1‑140 (MQ=255)
                                                                    ttGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGttt                                                           >  2:32604/1‑140 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  2:413705/140‑1 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  2:431248/140‑1 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  2:133458/140‑1 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  1:47009/140‑1 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  1:420359/140‑1 (MQ=255)
                                                                        tGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGaa                                                       <  2:586617/140‑1 (MQ=255)
                                                                         gAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAAc                                                      <  2:92820/140‑1 (MQ=255)
                                                                          aGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACa                                                     <  1:121650/140‑1 (MQ=255)
                                                                             tGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGa                                                  >  2:131915/1‑140 (MQ=255)
                                                                               cAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATa                                                <  1:131915/140‑1 (MQ=255)
                                                                                      tCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACt                                         >  1:517985/1‑140 (MQ=255)
                                                                                            gCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTc                                     >  1:142598/1‑140 (MQ=255)
                                                                                                        ccAACGCCTCGAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTATTGAACAAGATGAGGAACTGTGCCGTAACGCGCtgga                       >  1:168169/1‑140 (MQ=255)
                                                                                                                  gAAATCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGt             <  2:40027/140‑1 (MQ=255)
                                                                                                                   aaaTCTCGCCAGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTc            <  1:43927/140‑1 (MQ=255)
                                                                                                                             aGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCaaa    <  2:142598/140‑1 (MQ=255)
                                                                                                                             aGGCAACTT‑‑CGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCAAAgg  <  2:139718/140‑1 (MQ=255)
                                                                                                                                      ||                                                                                                                                   
GTACATGACGGCGGCGACCCAGGCGCGTGAACGTCTGCTGGCCGATAACGCCAAGAAGAAAGCGCAGATTGCTGAGTTGCAATCTTTCGTTAGCCGCTTTAGCGCCAACGCCTCGAAATCTCGCCAGGCAACTTCGCGCGCGCGCCAGATTGATAAAATCAAACTGGAAGAGGTGAAAGCCTCCAGCCGTCAGAACCCGTTCATCCGTTTTGAACAGGATAAGAAACTGTTCCGTAACGCGCTGGAAGTGGAAGGTCTGACCAAAGG  >  NC_000913/856679‑856945

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