Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,065,864 G→A 6.1% intergenic (‑100/+204) cobU ← / ← yeeH cobinamide kinase and cobinamide phosphate guanylyltransferase/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,065,8640GA 6.1% 212.3 / 3.9 65intergenic (‑100/+204)cobU/yeeHcobinamide kinase and cobinamide phosphate guanylyltransferase/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

GGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAT  >  NC_000913/2065719‑2065985
                                                                                                                                                 |                                                                                                                         
ggCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTa                                                                                                                         <  2:544987/148‑1 (MQ=255)
       ctTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    >  3:343912/1‑146 (MQ=255)
       ctTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    <  4:343912/146‑1 (MQ=255)
                         ccccGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATtct                                                                                                 >  3:1042673/1‑147 (MQ=255)
                         ccccGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATtct                                                                                                 <  4:1042673/147‑1 (MQ=255)
                                   aaaaTCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    <  3:1019870/118‑1 (MQ=255)
                                   aaaaTCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    >  4:1019870/1‑118 (MQ=255)
                                           catACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    <  1:537919/110‑1 (MQ=255)
                                           catACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    >  2:537919/1‑110 (MQ=255)
                                              aCCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATAGTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATg                                                                           >  1:304444/1‑148 (MQ=255)
                                                   cTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATAGTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCAtt                                                                      <  2:304444/148‑1 (MQ=255)
                                                           tCAAATGTATGATGGTATTCATTGAAGAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAATATGGCATTACGGAATa                                                              >  1:212859/1‑148 (MQ=255)
                                                           tCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAAt                                                               <  4:187803/147‑1 (MQ=255)
                                                           tCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAAt                                                               >  3:187803/1‑147 (MQ=255)
                                                                   atgatgGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACt                                                      >  2:1011687/1‑148 (MQ=255)
                                                                   atgatgGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCAAGACt                                                      <  1:1011687/148‑1 (MQ=255)
                                                                        ggTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATcc                                                                                                         >  3:809582/1‑92 (MQ=255)
                                                                        ggTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATcc                                                                                                         <  4:809582/92‑1 (MQ=255)
                                                                         gTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    >  2:401505/1‑80 (MQ=255)
                                                                         gTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    <  1:401505/80‑1 (MQ=255)
                                                                             tCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    >  1:899162/1‑76 (MQ=255)
                                                                             tCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAAtata                                                                                                                    <  2:899162/76‑1 (MQ=255)
                                                                             tCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTa                                                                                                                <  4:344130/80‑1 (MQ=255)
                                                                             tCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTa                                                                                                                >  3:344130/1‑80 (MQ=255)
                                                                              cATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTaaaa                                                                                                             >  4:205971/1‑82 (MQ=255)
                                                                              cATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTaaaa                                                                                                             <  3:205971/82‑1 (MQ=255)
                                                                                ttGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATaa                                                                                                                   <  4:392828/74‑1 (MQ=255)
                                                                                ttGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATaa                                                                                                                   >  3:392828/1‑74 (MQ=255)
                                                                                 tGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCGAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGc                                                                         <  3:585143/115‑1 (MQ=255)
                                                                                 tGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCGAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGc                                                                         >  4:585143/1‑115 (MQ=255)
                                                                                 tGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAata                                                                                                                      >  4:1125032/1‑70 (MQ=255)
                                                                                 tGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAata                                                                                                                      <  3:1125032/70‑1 (MQ=255)
                                                                                        aacaATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATg                                                                           <  3:499405/106‑1 (MQ=255)
                                                                                        aacaATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATg                                                                           >  4:499405/1‑106 (MQ=255)
                                                                                             tGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGAt                                                                            >  2:694084/1‑100 (MQ=255)
                                                                                             tGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGAt                                                                            <  1:694084/100‑1 (MQ=255)
                                                                                             tGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTa                                                                     >  4:911583/1‑107 (MQ=255)
                                                                                             tGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTa                                                                     <  3:911583/107‑1 (MQ=255)
                                                                                             tGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATg                            >  1:304549/1‑148 (MQ=255)
                                                                                                   ccccACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCAt                      <  2:304549/148‑1 (MQ=255)
                                                                                                          gTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGc               >  3:872155/1‑148 (MQ=255)
                                                                                                                  ccGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTc                                      <  3:120563/117‑1 (MQ=255)
                                                                                                                  ccGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTc                                      >  4:120563/1‑117 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCAt                                                                                         <  4:645608/61‑1 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCAt                                                                                         >  3:645608/1‑61 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATg                            >  3:773135/1‑122 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATg                            <  4:773135/122‑1 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGtt           <  3:618326/139‑1 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGtt           >  4:618326/1‑139 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  <  2:212859/148‑1 (MQ=255)
                                                                                                                       aaaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  <  4:872155/148‑1 (MQ=255)
                                                                                                                        aaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGtt                                            <  4:330182/105‑1 (MQ=255)
                                                                                                                        aaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGtt                                            >  3:330182/1‑105 (MQ=255)
                                                                                                                        aaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCa                       <  4:1140798/126‑1 (MQ=255)
                                                                                                                        aaCATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCa                       >  3:1140798/1‑126 (MQ=255)
                                                                                                                                cgcgATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGc                           >  4:931330/1‑114 (MQ=255)
                                                                                                                                cgcgATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGc                           <  3:931330/114‑1 (MQ=255)
                                                                                                                                           cTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  <  3:589564/128‑1 (MQ=255)
                                                                                                                                           cTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  >  4:589564/1‑128 (MQ=255)
                                                                                                                                            tGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATg                            >  4:989657/1‑101 (MQ=255)
                                                                                                                                            tGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATg                            <  3:989657/101‑1 (MQ=255)
                                                                                                                                            tGCGCATAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  >  4:745308/1‑127 (MQ=255)
                                                                                                                                            tGCGCATAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  <  3:303894/127‑1 (MQ=255)
                                                                                                                                            tGCGCATAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  <  3:745308/127‑1 (MQ=255)
                                                                                                                                            tGCGCATAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAt  >  4:303894/1‑127 (MQ=255)
                                                                                                                                                cgTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGtt           <  3:326963/114‑1 (MQ=255)
                                                                                                                                                cgTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGtt           >  4:326963/1‑114 (MQ=255)
                                                                                                                                                 |                                                                                                                         
GGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAACATTTACGCGATTAATACTGCGCGTAATATAATTAAAAATCCAGCATTCTCAATCCATTACGAAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGAT  >  NC_000913/2065719‑2065985

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