Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,576 T→C 12.1% intergenic (+21/+220) ybiT → / ← ybiU ABC transporter ATP‑binding protein YbiT/DUF1479 domain‑containing protein YbiU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5760TC12.1% 119.7 / 8.6 54intergenic (+21/+220)ybiT/ybiUABC transporter ATP‑binding protein YbiT/DUF1479 domain‑containing protein YbiU
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base C (4/3);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

TTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAA  >  NC_000913/857434‑857718
                                                                                                                                              |                                                                                                                                              
tttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGACGATTACCTGCGTAGTAAAGGGATAGAGTAAGTGATTATCCGACACCCCCAt                                                                                                                                                >  2:155702/1‑143 (MQ=255)
 ttCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATa                                                                                                                                               >  2:575379/1‑143 (MQ=255)
      acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTCGAAATCACCCCGGAACGCGTGATCGACTTTAGCTGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                          >  2:987016/1‑143 (MQ=255)
           ggTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAGCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                     <  1:1063197/142‑1 (MQ=255)
           cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                     <  2:539456/143‑1 (MQ=255)
           cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                     <  1:992232/143‑1 (MQ=255)
               aGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGCGAcg                                                                                                                                 >  1:188552/1‑143 (MQ=255)
                gTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGCAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcgg                                                                                                                                >  2:893854/1‑143 (MQ=255)
                gTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGGAAAGGGATCGAGTAAGTGATTATCCGTCACCACATTGCTGTCGGGGTGCcgg                                                                                                                                >  2:502359/1‑143 (MQ=255)
                 ttCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcggc                                                                                                                               >  1:469961/1‑143 (MQ=255)
                   cGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCTGTCACCCCGACAGTATCGGGGTGAcggata                                                                                                                             <  2:601945/143‑4 (MQ=255)
                    gtaatttCCCTCGCTAGCCGCATTCTGGAAATCCCCCCCTAACGGGTGATCTCCTTTAGCGGTAATTACTAAGCTTCCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                                                            <  2:323842/137‑1 (MQ=255)
                    gTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACACTGTCGGGGTGACgacggc                                                                                                                            >  1:601605/1‑143 (MQ=255)
                    gTAAGCTCCCTGGCGACCCGCATTCCGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                                                                            <  1:254895/143‑1 (MQ=255)
                     tAATCTCCCGTGCCACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACGTTAGCGGTAATTACGAAGATTACCGGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCt                                                                                                                           <  1:145942/143‑1 (MQ=255)
                       aGCTCCCTGGCGACCCGCGTTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                                                                         <  1:290904/143‑1 (MQ=255)
                         cTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgatg                                                                                                                       >  1:165940/1‑143 (MQ=255)
                           cccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGACCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                                                                     >  1:302307/1‑143 (MQ=255)
                            ccTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACGTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTGACCCCGATACTGTCGGGGTGACGGCGGCTGATGACa                                                                                                                    >  2:256183/1‑143 (MQ=255)
                                        cATTCTGGAAATCACCCCGAAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTAACTGCGTAGTAAAGGGATCGAGTAAGAGATTATCCGTCACCCCGATACTGTCGGGGTGACGACGGCTGATGACAAAAGGAAAATaa                                                                                                        >  2:941215/1‑142 (MQ=255)
                                           tCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACACTGTCGGGGCGACGGCGGCTCATGACAAATGCAAAATAGCCt                                                                                                     >  2:145983/1‑143 (MQ=255)
                                             tGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                                                                   >  2:408306/1‑143 (MQ=255)
                                                  aTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATgcgc                                                                                              <  1:52579/143‑1 (MQ=255)
                                                            agcgcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCGCCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                                                                    <  2:960373/141‑1 (MQ=255)
                                                               gcgTGATCGACGTTAGCGGTAATTACGAAGATTACCTGCGTGGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCGTAGCAgg                                                                                 <  2:991499/143‑1 (MQ=255)
                                                                cgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGc                                                                                <  1:987016/143‑1 (MQ=255)
                                                                 gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATACAAACAAGCCTGATGCGCTACGCATATCAgcca                                                                               >  2:835278/1‑140 (MQ=255)
                                                                    aTCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                            <  1:959182/143‑1 (MQ=255)
                                                                          tttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                                                                      <  1:16040/143‑1 (MQ=255)
                                                                              gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                                                                  >  1:670846/1‑143 (MQ=255)
                                                                                  tAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                              >  1:680265/1‑143 (MQ=255)
                                                                                  tAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                              >  1:680273/1‑143 (MQ=255)
                                                                                    aTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGCCGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGATTATCCGGCCTACGCAGTACCTGCAatgt                                                            >  2:743842/1‑141 (MQ=255)
                                                                                      tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat                                                          <  1:274065/143‑1 (MQ=255)
                                                                                        tcatgATTACCTGCGTAGTAAAGGGATTTATTATGTGATTTTGCTTCACCCCGATATTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg                                                        <  2:117856/139‑1 (MQ=37)
                                                                                            gATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATc                                                    >  1:713194/1‑143 (MQ=255)
                                                                                                   tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGa                                             <  1:186735/143‑1 (MQ=255)
                                                                                                    gCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGAc                                            >  1:410959/1‑143 (MQ=255)
                                                                                                     cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACt                                           <  1:979827/143‑1 (MQ=255)
                                                                                                        agtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATTAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTg                                        <  1:756681/143‑1 (MQ=255)
                                                                                                         gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCCAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACGTTGt                                       >  2:48640/1‑143 (MQ=255)
                                                                                                          taAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTtata                                      >  1:287973/1‑143 (MQ=255)
                                                                                                                 aTCGAGCAAGTGATTATCCGTCACCCCGAGACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGa                               <  1:74362/143‑1 (MQ=255)
                                                                                                                            gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTac                    >  2:573171/1‑142 (MQ=255)
                                                                                                                            gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTaa                    >  2:381375/1‑143 (MQ=255)
                                                                                                                               tagccGTCACCCCGATAGTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAAcgc                 <  2:415886/140‑1 (MQ=255)
                                                                                                                                  ccGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTCTCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCCGGACTTTGTAGGCCGGATAAGCCGTTCAcgccgc              >  2:322541/1‑143 (MQ=255)
                                                                                                                                   cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTAATCAAGCCTAAGCAGCACCTGCAATGTATTGACTCTGCAGGACTCTGTAGACCGGATCAGACGTTAAggccgaa             >  2:93271/1‑141 (MQ=255)
                                                                                                                                    gTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCAt            <  2:274769/143‑1 (MQ=255)
                                                                                                                                    gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATTCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCAt            >  1:894515/1‑143 (MQ=25)
                                                                                                                                     tCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc           >  1:578349/1‑143 (MQ=255)
                                                                                                                                     tCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATc           >  2:431050/1‑143 (MQ=255)
                                                                                                                                        ccccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAgg        <  1:524947/143‑1 (MQ=255)
                                                                                                                                         cccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCCTCAGGc       >  1:567172/1‑143 (MQ=255)
                                                                                                                                          ccGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGCTAAGGCGTTAACGCCGCATCAGGCa      >  2:1021590/1‑143 (MQ=255)
                                                                                                                                              tACTGTCAGGGTGACGGCGGCTGAAAACAAAAAAAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa  <  2:145397/143‑1 (MQ=25)
                                                                                                                                              |                                                                                                                                              
TTTCACACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAA  >  NC_000913/857434‑857718

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