Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,808,682 C→G 10.5% intergenic (+26/+15) yniA → / ← yniB putative kinase YniA/uncharacterized protein YniB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,808,6820CG10.5% 165.8 / 11.3 57intergenic (+26/+15)yniA/yniBputative kinase YniA/uncharacterized protein YniB
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base G (4/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

ATTTCCTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGT  >  NC_000913/1808538‑1808831
                                                                                                                                                |                                                                                                                                                     
aTTTCCTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCg                                                                                                                                                  <  1:561200/150‑1 (MQ=255)
   tCCTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGcc                                                                                                                                               <  1:293378/150‑1 (MQ=255)
    ccTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCACGGAGCGGcca                                                                                                                                              <  1:1113867/150‑2 (MQ=255)
     cTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCACGGAGCGGccat                                                                                                                                             <  1:1044945/150‑3 (MQ=255)
       tGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCtttt                                                                                                                                           >  2:1361423/1‑150 (MQ=255)
       tGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCtttt                                                                                                                                           >  1:986694/1‑150 (MQ=255)
        gAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTg                                                                                                                                          <  1:1125536/150‑1 (MQ=255)
          aCGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGAt                                                                                                                                        <  1:296234/150‑1 (MQ=255)
           cGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGAtt                                                                                                                                       <  1:1300033/150‑1 (MQ=255)
           cGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGAtt                                                                                                                                       <  2:406158/150‑1 (MQ=255)
            gTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATta                                                                                                                                      <  1:126388/150‑1 (MQ=255)
              cAACCGGTTTACAAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGAGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATtaaa                                                                                                                                    >  2:234243/1‑150 (MQ=255)
               aaCCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATtaaat                                                                                                                                   >  1:825837/1‑150 (MQ=255)
                       tACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCaa                                                                                                                           <  2:1116499/150‑1 (MQ=255)
                        aCCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAg                                                                                                                          >  1:1013203/1‑150 (MQ=255)
                         ccAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGc                                                                                                                         >  1:905287/1‑150 (MQ=255)
                                 acacGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATCTGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGaaaa                                                                                                                 <  2:1325797/150‑1 (MQ=255)
                                  cacGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAg                                                                                                                <  2:139250/150‑1 (MQ=255)
                                  cacGCTGATAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAg                                                                                                                >  2:1049749/1‑150 (MQ=255)
                                    cGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACACAAGCAGAGAAAAGaa                                                                                                              >  1:244437/1‑150 (MQ=255)
                                     gCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAac                                                                                                             <  2:242350/150‑1 (MQ=255)
                                                   cAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTgcagc                                                                                               >  2:770737/1‑150 (MQ=255)
                                                       gTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATa                                                                                           <  1:1044017/150‑1 (MQ=255)
                                                         tATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACa                                                                                         <  1:1361423/150‑1 (MQ=255)
                                                              ggCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGAGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATaa                                                                                    >  1:1257376/1‑150 (MQ=255)
                                                                 ggTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGTTAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCg                                                                                 >  2:970705/1‑150 (MQ=255)
                                                                    cAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCa                                                                              >  2:1065338/1‑150 (MQ=255)
                                                                         tttGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAATTAACCCCAGGAGGGGAAAGGAAACAATACCGGGAGGAAAAAAAATTACCGGCaaggatt                                                                        >  1:526762/1‑144 (MQ=255)
                                                                            ggTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATAc                                                                      <  1:896716/150‑1 (MQ=255)
                                                                            ggTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATAc                                                                      >  1:1388338/1‑150 (MQ=255)
                                                                             gTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACa                                                                     <  1:1049749/150‑1 (MQ=255)
                                                                                aTTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTg                                                                  <  2:791771/150‑1 (MQ=255)
                                                                                    cTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGaa                                                              >  1:251116/1‑150 (MQ=255)
                                                                                        gcagTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAAc                                                          <  2:196998/150‑1 (MQ=255)
                                                                                         cagTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACACAAGCAGAGAAAAGAACAAATAACCTGCAGCAATACAAATAACAGGCAGGATATACAGTGAGAAAAACt                                                         >  1:778485/1‑150 (MQ=255)
                                                                                            tCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCa                                                      <  2:972413/150‑1 (MQ=255)
                                                                                             cATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAg                                                     <  1:552430/150‑1 (MQ=255)
                                                                                             cATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAg                                                     >  1:40187/1‑150 (MQ=255)
                                                                                             cATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAg                                                     <  1:770737/150‑1 (MQ=255)
                                                                                             cATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAg                                                     <  2:825837/150‑1 (MQ=255)
                                                                                                  gATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGaaaa                                                <  1:823210/150‑1 (MQ=255)
                                                                                                  gATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGaaaa                                                >  2:1112806/1‑150 (MQ=255)
                                                                                                   aTAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAAt                                               <  1:796792/150‑1 (MQ=255)
                                                                                                         tattaGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATg                                         <  1:733904/150‑1 (MQ=255)
                                                                                                           ttaGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGt                                       <  2:986694/150‑1 (MQ=255)
                                                                                                                  caTGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGaac                                >  1:948717/1‑150 (MQ=255)
                                                                                                                     gATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGaacaac                             >  2:883369/1‑150 (MQ=255)
                                                                                                                               gAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGAc                   >  1:1355445/1‑150 (MQ=255)
                                                                                                                                aGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACt                  >  2:1430249/1‑150 (MQ=255)
                                                                                                                                 atcaaaaggCCGCACGGAGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTc                 >  1:716141/8‑150 (MQ=255)
                                                                                                                                  tcaaaaggCCGCACGGAGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCa                >  1:790748/7‑150 (MQ=255)
                                                                                                                                  tcaaaaggCCGCACGGAGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCa                >  2:173597/7‑150 (MQ=255)
                                                                                                                                  gATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCa                >  1:1014315/1‑150 (MQ=255)
                                                                                                                                  gATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCa                >  2:149348/1‑150 (MQ=255)
                                                                                                                                   caaaaggCCGCACGGAGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCaa               >  2:202518/6‑150 (MQ=255)
                                                                                                                                         gCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGtt         >  1:326133/1‑150 (MQ=255)
                                                                                                                                              tCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCAcgc    <  2:40187/150‑1 (MQ=255)
                                                                                                                                                cGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGt  >  1:1003744/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATTTCCTTGAACGTCAACCGGTTTACCAACTCTACACGCTGTTAAATCGTGCAAGGTTATTTGGCGGTCAGCATTTGGTTATTGCTCAGCAGTCATTGGATAGATTATTAGCAGCATGATATGGGTTGAGGATAATGGCCGCTCCGTGCGGCCTTTTGATTAAATAAACCCAAGCAGAGAAAAGAACACATAACCTGCAGCAATACAAATAACCGGCAGGATATACAGTGAGAAAAACTGCAGGAAAATAGTATGGTGCGGAACAACAATACGGGACTCAATTTGTTCACGCGT  >  NC_000913/1808538‑1808831

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