Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 522,463 T→C V683A (GTG→GCG)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913522,4630TC98.4% 215.7 / ‑5.6 61V683A (GTG→GCG) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/31);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCAC  >  NC_000913/522330‑522595
                                                                                                                                     |                                                                                                                                    
ttCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                 >  2:1162809/1‑139 (MQ=255)
  ccGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtggtg                                                                                                                               <  1:846613/139‑1 (MQ=255)
             aTGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCAc                                                                                                                    >  2:449537/1‑139 (MQ=255)
               ggCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCg                                                                                                                  <  2:298898/139‑1 (MQ=255)
                    cGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTAc                                                                                                                           <  2:148329/125‑1 (MQ=255)
                    cGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTAc                                                                                                                           >  1:148329/1‑125 (MQ=255)
                    cGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGc                                                                                                             >  1:737664/1‑139 (MQ=255)
                     ggCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCg                                                                                                            <  1:449537/139‑1 (MQ=255)
                          gTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       >  1:719469/1‑139 (MQ=255)
                                    cTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtgttgctgttgc                                                                                             <  1:658098/139‑1 (MQ=255)
                                            ccAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcag                                                                                     >  1:14564/1‑139 (MQ=255)
                                            ccAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcag                                                                                     >  2:1419418/1‑139 (MQ=255)
                                              aGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggt                                                                                   <  1:1162809/139‑1 (MQ=255)
                                                      aCCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       >  1:307398/1‑111 (MQ=255)
                                                      aCCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       <  2:307398/111‑1 (MQ=255)
                                                       ccATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGCCGGTCAGGTGCTGGAG‑AGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGa                                                                         <  2:838254/139‑1 (MQ=255)
                                                          ttAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAAt                                                                       >  1:385033/1‑139 (MQ=255)
                                                            aGCCTGTTAGATATTGGCACGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgc                                                                     >  2:837257/1‑139 (MQ=255)
                                                                  ttAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagc                                                               <  1:837257/139‑1 (MQ=255)
                                                                   tAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       >  2:1154237/1‑98 (MQ=255)
                                                                   tAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       <  1:1154237/98‑1 (MQ=255)
                                                                   tAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcg                                                              >  2:462557/1‑139 (MQ=255)
                                                                       tatTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGtcat                                                          >  1:471219/1‑139 (MQ=255)
                                                                          tGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtgttg                                                                                                    <  1:1131871/94‑1 (MQ=255)
                                                                          tGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtgttg                                                                                                    >  2:1131871/1‑94 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                 <  1:632119/79‑1 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                 >  2:632119/1‑79 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgcag                                                                               <  1:612993/113‑1 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgcag                                                                               >  2:612993/1‑113 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggt                                                                                   <  2:330211/109‑1 (MQ=255)
                                                                            gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggt                                                                                   >  1:330211/1‑109 (MQ=255)
                                                                                      aaaaCAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                           <  2:916121/139‑1 (MQ=255)
                                                                                           aGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGgcg                                      <  1:1422653/139‑1 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGCGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTa                                                                                                                            >  2:864039/1‑52 (MQ=38)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTa                                                                                                                            <  1:864039/52‑1 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGcac                                                                                       >  2:963979/1‑89 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGcac                                                                                       <  1:963979/89‑1 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                        >  1:14442/1‑88 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                        <  2:14442/88‑1 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                           >  2:51423/1‑133 (MQ=255)
                                                                                            ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                           <  1:51423/133‑1 (MQ=255)
                                                                                             gtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       >  1:350986/1‑72 (MQ=255)
                                                                                             gtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       <  2:350986/72‑1 (MQ=255)
                                                                                                  tcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTc                               <  1:462557/139‑1 (MQ=255)
                                                                                                       tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGgtgt                                          >  2:804792/1‑123 (MQ=255)
                                                                                                       tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGgtgt                                          <  1:804792/123‑1 (MQ=255)
                                                                                                                aGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGa                     <  1:156562/135‑1 (MQ=255)
                                                                                                                aGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGa                     >  2:156562/1‑135 (MQ=255)
                                                                                                                 ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgc                                                                                 <  1:1069945/74‑1 (MQ=255)
                                                                                                                 ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgc                                                                                 >  2:1069945/1‑74 (MQ=255)
                                                                                                                 ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                  <  2:1369972/73‑1 (MQ=255)
                                                                                                                 ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                  >  1:1369972/1‑73 (MQ=255)
                                                                                                                    aaGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcg                                                                    >  1:480219/1‑84 (MQ=255)
                                                                                                                    aaGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcg                                                                    <  2:480219/84‑1 (MQ=255)
                                                                                                                      gTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActgct           <  2:471219/139‑1 (MQ=255)
                                                                                                                            cgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                  >  2:1409194/1‑62 (MQ=255)
                                                                                                                            cgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                  <  1:1409194/62‑1 (MQ=255)
                                                                                                                            cgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTAc     >  2:1408044/1‑139 (MQ=255)
                                                                                                                               tGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCAc  <  1:1408044/139‑1 (MQ=255)
                                                                                                                                 gAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       >  1:1126989/1‑36 (MQ=255)
                                                                                                                                 gAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtg                                                                                                       <  2:1126989/36‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TTCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCAC  >  NC_000913/522330‑522595

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