Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,147 C→T 7.2% intergenic (‑130/+41) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1470CT 7.2% 182.5 / 3.9 57intergenic (‑130/+41)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base T (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATA  >  NC_000913/2027999‑2028293
                                                                                                                                                    |                                                                                                                                                  
tttttGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATACTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCacac                                                                                                                                                    >  2:239994/1‑149 (MQ=255)
 ttttGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACAcc                                                                                                                                                   <  4:395341/149‑1 (MQ=255)
     gTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGc                                                                                                                                               <  2:296088/149‑1 (MQ=255)
          gTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAt                                                                                                                                          <  1:575065/149‑1 (MQ=255)
                   ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                 >  3:150631/1‑149 (MQ=255)
                    cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGa                                                                                                                                           <  3:675393/138‑1 (MQ=255)
                    cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGa                                                                                                                                           >  4:675393/1‑138 (MQ=255)
                         gacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                           <  4:150631/149‑1 (MQ=255)
                           cTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACAGCAGGCCGGATGTGATCCGGCCAGATGt                                                                                                                         >  4:516125/1‑149 (MQ=255)
                           cTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                                                         <  2:601516/149‑1 (MQ=255)
                              ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCa                                                                                                                                                  <  2:103578/121‑1 (MQ=255)
                              ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCa                                                                                                                                                  >  1:103578/1‑121 (MQ=255)
                              ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGTCCTGATGTTaa                                                                                                                      <  3:516125/149‑1 (MQ=255)
                               cAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTTGCATGATAAGACGCGACTATCATTCATAGCAGGCATTAGCACACCAGGCCGAATGTGATCCGGCCAGATGTTAAt                                                                                                                     >  3:676924/1‑149 (MQ=255)
                                aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                    >  2:26533/1‑149 (MQ=255)
                                 aCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGa                                                                                                                   <  1:26533/149‑1 (MQ=255)
                                 aCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                               >  2:107277/1‑137 (MQ=255)
                                 aCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                               <  1:107277/137‑1 (MQ=255)
                                     gCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                      <  1:634856/126‑1 (MQ=255)
                                     gCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                      >  2:634856/1‑126 (MQ=255)
                                     gCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc                                                                                                               >  2:246503/1‑149 (MQ=255)
                                       ttttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCt                                                                                                             >  1:465305/1‑149 (MQ=255)
                                              tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                                                      >  2:198650/1‑149 (MQ=255)
                                               gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTggt                                                                                                                           >  4:662710/1‑125 (MQ=255)
                                               gATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTggt                                                                                                                           <  3:662710/127‑3 (MQ=255)
                                                aTTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                                                    >  1:730188/1‑149 (MQ=255)
                                                 ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                              >  2:281067/1‑122 (MQ=255)
                                                 ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                              <  1:281067/122‑1 (MQ=255)
                                                 ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                                                     <  4:511457/131‑1 (MQ=255)
                                                 ttCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                                                     >  3:511457/1‑131 (MQ=255)
                                                    gggCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                                                                <  3:48481/148‑1 (MQ=255)
                                                      gCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgt                                                                                                                         >  4:336925/1‑122 (MQ=38)
                                                      gCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgt                                                                                                                         <  3:336925/122‑1 (MQ=38)
                                                             cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                                                       >  2:378729/1‑149 (MQ=255)
                                                              tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct                                                                                                  >  1:308547/1‑137 (MQ=255)
                                                              tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct                                                                                                  <  2:308547/137‑1 (MQ=255)
                                                               ggCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAc                                                                                     <  1:246503/149‑1 (MQ=255)
                                                                 cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                                                                   <  1:293312/149‑1 (MQ=255)
                                                                  gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                  >  1:461241/1‑149 (MQ=255)
                                                                  gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                  >  4:447857/1‑149 (MQ=255)
                                                                   aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActccc                                                                                 <  1:239994/149‑1 (MQ=255)
                                                                                     gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTc                                                               <  2:461241/149‑1 (MQ=255)
                                                                                       gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                                                             >  4:793866/1‑149 (MQ=255)
                                                                                            ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                        >  2:175491/1‑149 (MQ=255)
                                                                                            ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                        <  2:406417/149‑1 (MQ=255)
                                                                                            ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                                                        <  4:409101/149‑1 (MQ=255)
                                                                                                  aaTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                  <  4:792505/149‑1 (MQ=255)
                                                                                                   aTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg                                                 >  3:351166/1‑149 (MQ=255)
                                                                                                              aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg                                      >  3:428207/1‑149 (MQ=255)
                                                                                                                       gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTaa                             <  2:730188/149‑1 (MQ=255)
                                                                                                                        aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAc                            <  2:465305/149‑1 (MQ=255)
                                                                                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCAt                      >  1:166874/1‑149 (MQ=255)
                                                                                                                               ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                     >  2:797016/1‑149 (MQ=255)
                                                                                                                                       aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGt                         >  4:186943/1‑137 (MQ=255)
                                                                                                                                       aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGt                         <  3:186943/137‑1 (MQ=255)
                                                                                                                                           aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc         <  2:809097/149‑1 (MQ=255)
                                                                                                                                           aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc         >  1:809097/1‑149 (MQ=255)
                                                                                                                                           aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc         >  1:486081/1‑149 (MQ=255)
                                                                                                                                           aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc         >  1:327689/1‑149 (MQ=255)
                                                                                                                                                  cacCAGGCCGGATGTGAGCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACCCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTCACCCCGTAACCGTCATACCAACGTCACTCAACCATa  <  4:676924/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATA  >  NC_000913/2027999‑2028293

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