Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 857,580 G→A 9.0% intergenic (+25/+216) 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,5800GA 9.0% 132.7 / 6.1 56intergenic (+25/+216)ybiT/ybiUABC transporter ATP‑binding protein YbiT/DUF1479 domain‑containing protein YbiU
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (2/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

ACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAA  >  NC_000913/857440‑857718
                                                                                                                                            |                                                                                                                                          
acGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTCGAAATCACCCCGGAACGCGTGATCGACTTTAGCTGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTc                                                                                                                                          >  2:987016/1‑143 (MQ=255)
     ggTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAGCGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                                                                     <  1:1063197/142‑1 (MQ=255)
     cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                     <  1:992232/143‑1 (MQ=255)
     cGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGt                                                                                                                                     <  2:539456/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:680273/1‑143 (MQ=255)
                                                                            tAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                                                              >  1:680265/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)
                                                                                                                                     cgtttcTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGAGGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGTACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCAt     <  2:7767/138‑1 (MQ=255)
                                                                                                                                        tccTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa  <  1:603028/141‑1 (MQ=255)
                                                                                                                                        tACTGTCAGGGTGACGGCGGCTGAAAACAAAAAAAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTaa  <  2:145397/143‑1 (MQ=25)
                                                                                                                                            |                                                                                                                                          
ACGACCGTGAGTTCGTAAGCTCCCTGGCGACCCGCATTCTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAA  >  NC_000913/857440‑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: